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I.   Introdnc t ion 

Free  space  is  a  commonly  used  transmission  medium  for  com- 
municating information.  Information  communicated  as  such  can  not 
be  kept  confidential  because  of  the  nature  of  the  medium.  Thus 
the  task  of  surveillance  is  established;  there  will  always  be 
someone  that  wishes  to  receive  confidential  information  intended 
for  someone  else.  Many  times  a  surveillance  receiver  is  not  in 
the  optimal  location  for  reception  which  results  in  low  signal- 
to-noise  ratios.  Also,  the  surveillance  receiver  can  not  expect 
to  know  a  priori  any  parameters  of  the  transmission.  Further- 
more, those  transmitting  the  information  will  try  to  keep  it  from 
getting  to  the  wrong  people.  For  these  three  reasons  alone, 
surveillance  receivers  are  among  the  most  challenging  receivers 
to  design. 

This  research  effort  is  focused  at  determining  the  baud  rate 
of  a  message  transmitted  using  frequency  shift  keying  (FSK) 
modulation.  The  received  signal  is  assumed  to  be  FSK.  This 
assumption  has  been  made  either  by  a  modulation  recognition 
algorithm  or  by  some  knowledgeable  insight.  Instead  of  relying 
on  conventional  methods  which  use  a  zero-crossing  algorithm,  the 
estimate  is  determined  using  adaptive  noise  cancelling  tech- 
niques. Figure  1  shows  the  two  methods  of  a  baud  rate  analysis 
system.  The  crux  of  the  problem  is  to  evaluate  the  performance 
of  the  adaptive  noise  cancelling  technique  in  estimating  the  baud 
rate  in  a  low  s igna 1-to-soise  ratio  environment. 
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Obviously  some  restrictions  must  be  made  since  one  receiver 
can  not  be  expected  to  work  for  all  signals.  The  first  restric- 
tion is  that  the  FSK  signal  bandwidth  be  no  greater  than  6000  Hz. 
This  fignre  was  arbitrarily  chosen  bnt  can  be  changed  to  meet  the 
needs  of  a  specific  application.  The  other  restriction  applies 
to  the  band  rate.  The  FSK  signal  bandwidth,  the  frequency  devia- 
tion of  the  FSK  signal,  and  the  band  rate  are  ronghly  related  by 

BFSK  &   2(AfFSK  +  R)  (1.1) 

In  order  to  demodulate  the  signal  with  some  degree  of  accuracy 
the  frequency  shift  must  be  as  large  as  the  baud  rate.  This 
means  that  the  maximum  baud  rate  for  a  given  signal  bandwidth  is 
one-fourth  the  bandwidth.  For  this  case  the  main  lobes  of  the 
message  spectrum  are  above  and  below  the  center  frequency.  To 
lower  the  probability  of  error  during  demodulation,  the  frequency 
deviation  is  generally  larger  than  the  band  rate.  In  order  that 
the  main  lobe  and  the  first  side  lobe  of  the  message  spectrum  be 
above  and  below  the  center  frequency  upon  modulation,  the  fre- 
quency deviation  must  be  at  least  twice  the  baud  rate.  Based  on 
this  premise,  the  restriction  for  the  maximum  baud  rate  is  1000 
Hz  given  the  6000  Hz  signal  bandwidth. 

Since  the  investigation  will  be  done  by  comparing  the  per- 
formance of  the  proposed  method  against  that  of  the  conventional 
method,  any  benefits  or  harms  that  may  result  from  the  arbi- 
trarily chosen  restictions  will  affect  both  methods.  Therefore, 
final  judgment  will  not  be  influenced  by  these  restrictions. 


II.   FSK  Band  Ra.te  Ana  lysis  System 

The  system  nsed  to  determine  the  band  rate  of  an  FSK  signal 
consists  of  two  major  devices  as  shown  in  Figure  1.  The  first  is 
an  FSK  demodulator,  the  second  is  the  band  rate  estimator.  The 
band  rate  estimate  is  determined  either  with  adaptive  noise 
cancelling  techniques,  the  method  nnder  investigation,  or  with 
zero-crossing  techniques,  the  conventional  method.  The  FSK  demo- 
dulator, illustrated  in  Figure  2,  has  a  delay-line  discriminator 
(OLD)  as  its  principal  component.  Before  the  DLD  is  a  bandpass 
filter  to  limit  the  input  noise  power,  and  following  the  DLD  is  a 
low-pass  filter  for  the  same  purpose.  The  limiter  is  to  stabil- 
ize the  DLD  response  which  is  dependent  on  the  input  signal 
ampl i tude . 
Bandpa  ss  Filter 

The  pre-de t e c t i on  filter  is  a  4-pole  Butterworth  bandpass 
filter.  Its  purpose  is  twofold.  First,  it  limits  the  noise 
power  that  enters  the  system.  Second,  it  allows  only  signals  for 
which  the  receiver  is  designed  to  operate  to  enter  the  system. 
For  this  reason  the  bandwidth  is  set  to  6000  Hz.  Also,  the 
filter  should  be  centered  on  the  carrier  frequency  of  the  signal 
which  has  been  determined  a  priori. 
Limi t e  r 

FSK  is  a  special  case  of  frequency  modulation  (FH)  and  the 
FSK  demodulator  can  actually  be  considered  as  an  FM  demodulator. 
It  makes  sense  then  to  eliminate  any  amplitude  modulation 
resulting  from  the  noise  perturbations  which  might  have  entered 
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the  signal.  Furthermore,  as  will  be  discussed  shortly,  the 
output  of  the  DLD  for  a  continuous-wave  input  is  dependent  upon 
the  squared  amplitude  of  the  input.  Thus,  if  the  input  amplitude 
is    left   to  vary,    the  DLD  reponse  will    also  vary. 

The    limiter    is    modeled    as    a    device    which   will    have    an    output 
of    constant    amplitude    but    will    have    identical    phase    properties    as 
the    input.       That    is,     the    limiter   will    block    any    amplitude    modula- 
tion  while    frequency    modulation    is    allowed    to    pass. 
De lav-Line    Discriminator 

The    DLD     is     a    device    whose    output     is    proportional     to    the 

instantaneous    frequency   deviation   of    the    input    signal     from    its 

center    frequency.       A    block    diagram    of     the    DLD     is     given     in    Figure 

3.       One    use    of    the    DLD    is    for    FSK    demodulation.       If    the    input    is 

/ 

an  FSK  signal  the  output  will  be  the  bit  sequence  of  the  binary 
me  s  sage . 

The  effect  of  the  DLD  is  apparent  when  the  input  consists  of 
a  noise-free  unmodulated  carrier,  namely, 

Uj(t)  »  A  cos[(2nfc  +  2nAf)t]  (2.1) 

The  frequency  offset  from  the  carrier  frequency  of  the  signal, 
f  ,  is  denoted  by  Af.  It  can  be  shown  that  the  output  of  the 
DLD  is  given  by 

m(t)  =  -4A2sin(2nAfT  +  2nfcT)  (2.2) 

where  T  is  the  time  of  delay  between  oj(t)  and  u2ft)-  Note  this 
output  is  only  for  a  continuous-wave  input.  This  result  plus 
the  output  of  the  DLD  for  a  general  input  are  d  e  v  e  1  op  e  d  w  i  t  h  low- 
pass  modeling  techniques  presented  in  the  following  section. 
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In  order  for  the  discriminator  to  work  properly,  (i.e.,  pro- 
dace  an  output  proportional  to  in s t a n t a ne on s  input  frequency 
deviation)  it  is  necessary  that  that  the  output  be  zero  when  the 
deviation  is  zero.   For  this  to  be,  it  is  necessary  that 

2nfcT  =  nit      for   n-1,2,...  (2.3) 

For  this  choice,  the  output  is  described  by 

m(t)  =  +4A2sin(2nAfT)  (2.4) 

The  +  cones  from  the  addition  of  nn,  +  for  n  odd  and  -  for  n 
even.  A  plot  of  one  cycle  of  (2.4)  is  shown  in  Figure  4  when  the 
delay  T  is  chosen  to  be  125  usee  and  n  is  even.  Observe  that  the 
curve  is  nearly  linear  when  the  frequency  deviation  is  small.  T 
should  be  chosen  so  that  the  discriminator  operates  in  this 
linear  range.  The  obvious  reason  for  the  choice  is  that  it  is 
easier  to  work  with  linear  or  nearly  linear  devices  as  opposed  to 
non-linear  devices.  Another  reason  is  that  the  response  of  the 
DLD  in  the  non-linear  region  for  general  inputs  is  not  understood 
well.  Since  the  pre-filter  has  limited  the  range  of  the  incoming 
deviations,  the  bandwidth  of  the  pre-filter  dictates  a  specific 
delay.   The  procedure  in  setting  T  is  now  discussed. 

The  linear  region  is  roughly  that  portion  of  the  spectrum 
centered  on  f  such  that  the  argument  of  the  sinusoid  is  less 
that  it/6  in  magnitude.  Then  the  relation  between  the  maximum 
delay  and  the  maximum  deviation  is 


(2.5) 
t 

Recall     that     the     pre-filter    passes     frequencies    from    ~Af_a       to 
Afmai.        Substituting    the    pre-filter    bandwidth,     Bgp,     for    twice    the 
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maximum  deviation  and  then  solving  for  the  maximum  delay  results 
in 


6B 


RF 


(2.6) 


Since    the    product    2fcT    must    be    an     integer    by    Eqn.     (2.3),     T    should 

be    chosen    as    the     largest    value    which    satasfies    Eqns.     (2.3)    and 

(2.6). 

Low-Pass  Filter 

The  purpose  of  the  post-filter  is  to  reduce  the  noise  in  the 
output  signal  of  the  DLD.  Since  the  output  is  at  baseband,  the 
filter  is  low-pass.  The  filter  implemented  is  a  4-pole  Butter- 
worth  filter.  The  bandwidth,  to  eliminate  as  much  noise  as 
possible,  should  be  only  slightly  larger  than  the  bandwidth  of 
the  output  signal.  However,  selecting  the  bandwidth  is  difficult 
since  the  signal  bandwidth  is  proportional  to  the  baud  rate.  As 
stated  in  the  introduction,  the  system  should  work  for  baud  rates 
up  to  1000  Hz.  Thus  the  post-filter  bandwidth  is  set  to  1000  Hz. 
Adaptive  Noise  Cance  Her 

The  method  of  baud  rate  analysis  under  study  consists  of  an 
adaptive  noise  canceller  (hereafter  referred  to  as  an  ANC)  and 
the  ANC  baud  rate  estimator.  The  ANC  baud  rate  estimator  is  an 
algorithm  which  evaluates  the  parameters  of  the  ANC  to  determine 
in  estimate  of  the  baud  rate.  The  ANC  implemented  is  given  in 
Figure  5  and  was  introduced  by  Yidrow  [1],  The  primary  input  to 
the  ANC  is  some  signal  plus  noise,  s+n~.  The  reference  input, 
specified  as  n j ,   is  also  noise.    It  is  uncorrelated  with  the 
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signal  s  but  is  correlated  with  uq  in  some  (possibly)  unknown 
way.  The  reference  input  nj  is  filtered  by  the  adaptive  filter 
to  produce  y  which  is  as  close  a  replica  of  ng  as  possible.  The 
output     z     is     the     difference    of    the    primary    input    s+n*    and    y. 

The  adaptive  filter  of  the  ANC  differs  from  a  conventional 
filter  in  that  its  impulse  response  is  automatically  adjusted  to 
produce  the  'best  result'  which  in  this  case  is  noise  cancella- 
tion. The  adjustment,  or  adaptation,  is  accomplished  through  an 
algorithm  that  responds  to  the  ANC  output.  Thus,  the  ANC  will 
operate  under  changing  conditions  or  in  situations  in  which  the 
noise  can  not  be  characterized  (something  that  a  conventional 
noise  cancelling  filter  depends  on).  Widrow  [1]  provides  a 
simple    but    convincing    argument    to     support     this     claim. 

The  ANC  output,  z=s+nQ-y,  is  referred  to  as  the  error  of  the 
system.  The  object  of  this  system  is  to  produce  an  output  z  that 
is  the  best  fit  in  the  least  squares  sense  to  'the  signal  s. 
Therefore,  the  object  is  to  minimize  the  mean  squared  error.  It 
can  be  shown  that  doing  so  does  not  reduce  the  signal  power  of 
the  output  but  rather  it  reduces  the  average  power  of  the  quan- 
tity nQ-y.  This  means  that  minimizing  the  error  must  maximize 
the  adaptive  filter's  objective  of  transforming  n«  into  y. 
Furthermore,  since  the  signal  power  remains  constant,  minimizing 
the  total  output  power  maximizes  the  output  signal-to-noise 
power     ratio. 

These  arguments  can  be  extended  to  cases  where  the  primary 
input    and    the    reference     input     do    not     consist     of    just     signal     plus 
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noise  and  noise,  respectively.  As  long  as  the  reference  input  is 
correlated  with  part  of  the  primary  input  and  uncorrelated  with 
the  rest,  the  two  parts  of  the  primary  input  can  be  separated. 

Based  on  this  premise,  a  scheme  was  devised  to  use  the  ANC 
to  separate  the  noise  from  the  message  in  the  DLD  output.  The 
derivation  of  the  scheme  required  information  about  the  correla- 
tion properties  of  the  noise  and  message  at  the  output  of  the 
DLD.  Ratcliffe  [2]  treats  this  subject  extensively.  Another 
approach  at  obtaining  this  information  is  to  use  a  simulation 
model  as  described  in  the  next  two  sections.  This  was  the  method 
used  to  obtain  the  information. 

First  a  noise-free  signal  consisting  of  an  FSK  modulated, 
pseudo-random  sequence  was  generated.  Next  this  signal  was 
passed  through  the  FSK  demodulator.  The  input  to  the  ANC  was 
then  recorded.  This  input  is  either  the  direct  output  of  the  DLD 
or  the  filtered  output  of  the  DLD  (depending  on  whether  post- 
filtering  was  performed  or  not).  In  either  case,  the  input  is 
some  variation  of  the  binary  message  bit  sequence.  To  avoid 
ambiguities,  the  result  will  be  referred  to  as  the  ANC  input. 
The  message  autocorrelation  estimate  was  obtained  from  an  ANC 
input  which  resulted  from  a  noise-free  FSK  signal  with  the 
estimation  being  done  with  FFT  techniques.  For  the  noise  auto- 
correlation estimate,  some  definition  of  the  noise  had  to  be 
made.  The  noise  at  the  ANC  input  was  arbitrarily  defined  as  the 
difference  of  the  input  resulting  from  a  noisy  FSK  signal  and  the 
input  resulting  from  a  noise-free  FSK  signal.   The  final  noise 
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autocorrelation  estimate  was  produced  by  averaging  ten  individual 
estimates.  This  procedure  was  completed  for  several  predetection 
signal-to-noise  ratios  and  for  the  cases  of  post-filtering  and  no 
post-filtering. 

Figure  6  shows  the  message  autocorrelation  plotted  with  the 
noise  autocorrelation  for  a  SNR  of  10  dB  and  with  a  post-filter 
bandwidth  of  1000  Hz.  Figure  6  also  shows  the  same  results  when 
the  post-filter  was  omitted.  Since  the  functions  are  even,  only 
the  positive  components  are  shown.  Furthermore,  the  estimates 
are  plotted  for  lag  components  0  to  40  since  the  rest  of  the 
function  contains  no  additional  information.  For  different 
SNR's,  the  autocorrelation  functions  differed  only  in  a  scaling 
of  the  amplitude  of  the  noise  autocorrelation.  Note  that  in  both 
the  filtered  and  non-filtered  cases  the  message  is  correlated 
until  component  numbe  r  3  2,  which  is  the  numbe  r  of  samples  per  bit 
for  this  case  (the  hump  past  component  32  in  the  filtered  case  is 
caused  by  the  correlation  introduced  by  the  post-filter).  This 
result  is  intuitively  obvious  since  the  sequence  is  pseudo- 
random. The  noise  autocorrelation  for  the  filtered  case  shows 
the  noise  is  practically  uncorrelated  past  component  20.  For  the 
non-filtered  case  the  correlation  extends  to  only  the  sixth 
component.  Note  that  the  post-filter  has  eliminated  some  noise 
power  as  evidenced  by  the  difference  in  peak  amplitudes  of  the 
noise  autocorrelation  for  the  two  cases. 

Now  suppose  the  primary  input  of  the  ANC  is  the  output  of 
the  FSK  demodulator,   i.e.,  both  message  and  noise.   The  ANC  will 
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t,   -   9600  Hz 


FIGURE  6:  Effect  of  post-filter  and  limiter  on  the 

correlation  properties  of  the  DLD  output. 
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cancel  one  of  these  two  if  the  reference  input  is  correlated  with 
one  and  is  nncorrelated  with  the  other.  The  portion  being  can- 
celled will  be  the  one  that  is  correlated  with  the  reference 
inpnt.  Notice  that  when  the  DLD  output  without  post-filtering  is 
delayed  by  8  samples  the  noise  portion  of  delayed  seqnence  is 
nncorrelated  with  the  noise  of  the  non-delayed  seqnence. 
Furthermore  the  message  portion  of  the  delayed  and  non-delayed 
sequences  still  remains  correlated.  Hence  the  ANC  will  attempt 
to  cancel  the  message  in  the  DLD  output,  i.e.  upon  adapting,  the 
adaptive  filter's  output  will  be  the  message  and  the  output  error 
will  be  the  noise.  The  weights  of  the  adaptive  filter  will  have 
converged  once  the  noise  of  the  error  becomes  stationary.  In 
this  mode,  the  ANC  is  being  used  as  a  self-tuning  filter.  This 
is  illustrated  in  Figure   7. 

Suppose  the  filtered  output  is  considered.  The  delay  would 
have  to  be  20  samples  instead  of  6  for  the  noise  to  become 
nncorrelated.  But  note  that  a  delay  this  long  weakens  the  corre- 
lation between  the  signal  portion  of  the  two  inputs  making  the 
job  of  cancelling  more  difficult.  Thus,  omitting  the  post-filter 
for  this  method  seems  like  a  good  idea  especially  since  deter- 
mining the  correct  bandwidth  is  very  difficult  as  mentioned 
earlier.  One  argument  for  keeping  the  post  filter  is  that  the 
magnitude  of  correlation  of  the  noise  is  reduced,  which  would 
improve  the  results  of  the  cancelling  tasks.  The  initial 
approach  taken  is  to  omit  the  post-filter  when  making  baud  rate 
estimates  with  the  ANC  and  to  use  a  delay  long  enough  so  the 


17 


e 


I 

3 

4-1 

I 


<e 


03 

in 


(0 


c- 

0) 


03 

>  c 

•d  QJ 

JJ  JJ 

a  r-l 

(0  -rt 

TT  <*- 

<o 


0) 

o 

c 

(0 

u 

m 
tn 
•rt 

a 
c 

o 

> 


a 

10 
■D 
ID 

a 

.c 
H 


S 
3 

CD 


18 


noise  is  uncorre 1  a  ted. 
Zero-Crossing  Est  ima tor 

The  proposed  method  of  band  rate  analysis  is  to  be  measured 
against  the  traditional  method  of  determining  the  FSK  bend  rate, 
that  being  the  observation  of  the  zero  crossings  of  a  demodulated 
FSK  signal.  The  concept  is  that  a  zero  crossing  in  the  message 
should  only  occur  at  intervals  equal  to  the  baud  rate  or  integral 
multiples  of  the  baud  rate.  Obviously  noise  will  perturb  the 
message  and  thus  the  zero  crossing  might  not  be  exactly  in  the 
proper  place.  Thus  the  algorithm  should  make  decisions  based  on 
some  criteria  in  order  to  obtain  a  baud  rate  estimate.  This 
algorithm  is  described  in  the  next  section. 


Ill .   Deve  lopment  £.f  a    System  Mode  1 

A  mathematical  model  of  the  system  must  be  developed  in 
order  to  simulate  the  system.  Since  the  simulation  is  performed 
on  a  digital  computer,  the  implications  of  Shannon's  sampling 
theorem  must  be  kept  in  mind.  Obviously,  since  the  received  FSK 
signal  is  a  bandpass  signal,  sampling  in  order  to  prevent  alias- 
ing would  require  vast  amounts  of  memory  for  any  practical  length 
of  received  signal.  Besides  requiring  an  enormous  memory  capa- 
city, the  simulation  would  also  demand  much  time  to  complete. 
Obviously  the  modeling  cannot  be  done  efficiently  using  bandpass 
representations.  The  solution  is  to  use  low-pass  modeling  tech- 
niques. The  technique  will  be  applied  to  the  pre-filter,  the 
limiter,  and  the  DLD  since  they  are  the  only  bandpass  devices  in 
the  system.  All  other  devices  are  baseband  and  can  be  modeled  as 
they  are.  Thus  the  FSK  demodulator  illustrated  again  in  Figure 
8a  is  modeled  by  its  low-pass  equivalent  representation  shown  in 
Figure  8b. 
Low  P.8.S.S  Technique  s 

Low-pass  modeling  [3]  relies  on  the  property  that  the  per- 
tinent information  of  a  band-pass  signal  is  contained  in  its 
magnitude  and  in  its  phase.  Both  of  these  are  slowly  varying  in 
relation  to  the  carrier  frequency  and  thus  are  easier  to  simulate 
with  a  digital  computer.  The  relation  between  a  bandpass  signal 
s(t)  and  its  equivalent  low-pass  counterpart  or  complex  envelope 

*(t)  is,  by  definition, 

j2nfct 
s(t)  -  Re{s(t)e       }  (3.1) 
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If  the  complex  envelope  is  expressed  in  terms  of  its  magnitude 

and  argument  then  it  follows  that 

s(t)  =   ?(t)  cos[2nfct  +  e(t>]  (3.2) 

where   s(t)   is  the  magnitude  of  s(t)  and  i  (  t )  is  the  argument. 

Evidently  s(t)   is  a  complex  function  with  the  property  that  its 

magnitude  is  the  envelope  of  s(t)  and  its  argument  is  the  phase 

of  s(t).   All  bandpass  functions,  both  deterministic  and  random, 

can  be  put  in  this  format. 

Bandpass  Pre-Fil tor 

A  bandpass   filter  with  an  impulse  response  h(t)  has  an 

equivalent  low-pass  impulse  response  E(t)  which  by  definition 

satisfies 

j2nfct 
h(t)  =  2  Re{n(t)e       }  (3.3) 

where  f  is  the  center  frequency  of  the  filter.  An  equivalent 
low-pass  transfer  function,  H(f),  can  be  derived  from  this  defi- 
nition [3].   It  is 

H(f)  =  H(f-fc)  +  I*(-f-fc)  (3.4) 

where  H(f)  and  H(f)  are  the  Fourier  transforms  of  h(t)  and  E(t), 
respectively.  Also,  represents  conjugation.  This  relationship 
can  be  shifted  in  frequency  to  obtain  the  equivalent  low-pass 
transfer  fnnction  in  terms  of  the  actual  bandpass  transfer  func- 
tion 

i(f)  =  IH(f+fc)]low_pass  term  (3.5) 

The  low-pass  model  of  the  filter  will  have  a  bandwidth  equal  to 


one-half  the  RF  bandwidth.   i.e.   Bj   =  BRF/2.    All   transform 
relations  for  linear  systems  are  also  valid,  for  low-pass  models. 
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e.g.,   for  an  input  of  X(f)  and  a  transfer  function  H(f)  the 

output  T(f)  is 

Y(f)  =  i(f)I(f)  ,  (3.6) 

Limi ter 

The  limiter  model  with  input  u(t)  has  an  output  Ujft)  which 

is  related  to  the  input  by 

j  arg[u(t)] 
Ij(t)  =  A  e  (3.7) 

where  A  is  some  constant  amplitude. 

Delay-Line  Discriminator 

The  block  diagram  of  the  low-pass  equivalent  model  of  a 

delay-line  discriminator  is  shown  in  Figure  9.   The  output  of  the 

delay  block  can  be  characterized  in  terms  of  its  input  by 

u2(t)  =  UjU-TJe      c  (3.8) 

This  can  be  verified  by  expanding  the  definition  of  the  complex 
envelope  of  Uj(t)  =  Uj(t-T). 

Each  output  of  the  90  hybrid  is  equal  to  the  corresponding 
input  minus  a  90  phase  shifted  copy  of  the  remaining  input.  In 
equations  this  can  be  stated  as 

'jltl  ■  Ijltl  -  J  Jjltl  (3.9a) 

z2(t)  -  u2(t)  -  j  ni(t)  (3.9b) 

The  action  of  the   square   law  detector  is  to  produce  an 

output  proportional   to  the  magnitude  squared  of  the  complex 

envelope  of  the  input,  i.e. 

(3.10a) 
(3.10b) 


v,(t)  -ME,  (t) 
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Note  that  the  output  is  no  longer  a  bandpass  signal.  Thns  the 
low-pass  model  is  the  sane  as  the  physical  system  from  this  point 
on. 

Finally     the     output     of     the     DLD     is     simply     the     difference 
between    the    two    square    law    detectors 

n(t)    =    Vjft)    -    v2(t)  (3.11) 

The     discriminator    equation     is     derived    as     follows.       The     out- 
puts   of    the    90       hybrid    in    terms    of    the    input    Ujft)    are 

-j2*fcT 


ZjU)     =    Ujft)     -    j     u1(t-T) 

-j2nfcT 
z2(t)    =    Ujtt-T)    e  -    j    SjU) 


(3.12a) 


(3.12b) 

From  these,  the  output  of  the  first  square  law  envelope  detector 
is  determined  to  be 


Vj(t)  =  |  z^t)! 


«j(t)»J(t) 


[u1(t)-jS1(t-T)e       ][5£(t)+jS£(t-T)« 


"l(t,l  2  +  ■>  *i<t)»J(t-T)« 


j2nfcT 


+       ux(t-T)     z    -    j    Uj(t)51(t-T)( 
Similarly 
v2(t)    =   |   ujft)!  2    -    j    UjCOuJU-TJe 

+  |u1(t-T)|2    +    j    uJUJtijU-Tje 
The     comparator    output     is    then 
m(t)    =    v^t)    -    v2(t) 


-j2nfcT 


(3.13a) 

f 
]       (3.13b) 

(3.13c) 


J2nf    T 


2ju, (t)u. (t-T)e 


j2nfcT 


-    2JU,  (t)u..  (t-T)e 


"j2nfcT 


(3.14) 


(3.15a) 


(3.15b) 
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2[u1<t)S1(t-T) 


j2nfcT  ir  .    _.      i2nfcT  Jr 

i       e  *+  n1(t)u1(t-T)e       e  *]      (3.15c) 


_♦       j2nfcT  j-  J2nfcT  j- 

2{«1(t)«j{t-T)l       e  2  +  [n1(t)n1(t-T)e       e  2]*)   (3.15d) 


4Re{n1(t)SJ(t-T)e 


J<2nfcT  +  T> 

2  }  (3.15e) 

To    observe    the    discriminator    effect    for    a    c on t i nuo u s -w a v e 

signal,     let    the     input    be 

Uj(t)    =    A    cos[2n(fc    +    Af)t]  (3.16) 

then 

j2rtAft 
Ujft)  =  A  e  (3.17) 

Substituting  into  the  output  expression  for  a  general   input, 

(3.15e),   gives 

j2«Aft     -j2nAf(t-T)   j(2nfcT  +  -) 
m(t)  =  4Re{[Ae       ] [Ae  ]e  2  )   (3.18) 

Simplifying  yields 

m(t)     =    4A2cos[2it(Af    +    f0)T    +    it/2] 


-4A2sin[2n(Af    +    fc)T] 


(3.19a) 
(3.19b) 

This  is  the  result  given  and  discnssed  in  the  DLD  description  of 

the  previous  section. 

It  is  important  to  note  that  the  input  u.(t)  will  consist  of 

two  terms  in  the  general  case,  a  signal  envelope  and  a  noise 

envelope.    Ratcliffe   [2]  discusses  the  statistics  and  the  power 

spectra  associated  with  the  general  case. 

Low-pass  Pos t-Pll ter 

The  post-filter,  being  lowpass,  can  be  implemented  directly 

into  the  system  model.   No  equivalent  model  is  necessary. 
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Adapt  ive  Noise  Cance  Her 

Tie  adaptive  filter  of  the  canceller  is  so  named  because  the 
filter  weights  adjnst  to  changing  conditions.  The  implementation 
of  the  adaptive  filter  and  the  algorithm  for  the  weight  adjust- 
ment will  now  be  discussed.  The  discusion  is  an  abbreviated 
version  of  the  one  in  Widrow  [1], 

The  adaptive  filter  is  modeled  with  the  adaptive  linear 
combiner  shown  in  Figure  10.  Note  that  X,  and  W  are  vectors, 
while  y.  and  d.  are  scalars.   The  vectors  are  defined  as 

Xj  =  [x0j  i1;j  .  .  .  xnj]T  (3.20) 

W  =  twQ  Wl  .  .  .  wn]T  (3.21) 

Generally  xq,  is  set  to  1  so  that  wg  can  be  considered  the 
biasing  weight.  The  output  y,  is  the  inner  product  of  X,  and  W. 
The  error  e  .  is  the  difference  between  the  desired  response  d, 
and  the  actual  response  y,.   That  is 

7j  =  XJW  =  WTXj  (3.22) 

6j  =  dj  "  Xj  "  dj  "  *TXJ  (3-23) 

Note  that  in  accordance  with  the  system  description  of  the  ANC, 

primary  input  and  the  reference  input. 

Since  the  desire  is  to  minimize  the  mean  squared  error,  the 
mean  squared  error  is  derived  first.   It  is 

Ete^]  =  E[dJ]  -  2E[djx]']W  +  WTE[Xjx]']W  (3.24) 

where  E[*]  denotes  expectation.  The  gradient  of  the  error 
function,  denoted  by  V,  is  obtained  by  differentiating  this 
result  with  respect  to  f.   The  gradient  is 
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V  =  V{E[e2^,  =  -2EtdjlJ]  +  2E[X.jXT]W  (3.25) 

The  optimal  weight  vector  W  is  obtained  by  setting  the  gradient 
to  ze  ro. 

W*  =  {E[Xjx]']}":lErdjxT]  (3.26) 

The  idea  behind  using  the  gradient  method  is  that  the  mean 
squared  function  is  a  quadratic  function  of  the  weights.  It  can 
be  pictured  as  a  concave  hype rpa r abo 1 o ida 1  surface  that  never 
goes  negative,  i.e.  a  bowl  resting  on  some  point  in  the  hyper- 
space.  The  object  is  to  get  to  the  bottom  of  the  bowl,  at  which 
point  the  gradient  is  the  zero  vector. 

The  LHS  algorithm  implements  the  method  of  steepest  descent. 
According  to  this  method,  the  next  weight  vector  estimate  W,  +  »  is 
equal  to  the  present  weight  vector  estimate  W.  pins  a  change 
proportional  to  the  negative  of  the  gradient;  i.e., 

Wj+1  =  Wj  "  *Tj  (3.27) 

The  parameter  u  controls  stability  and  convergence  rate.  The 
gradient  V.  is  the  true  gradient  evaluated  with  the  weight  vector 
of  the  jth  iteration. 

The  LHS  algorithm  assumes  i?  as  an  estimate  of  the  mean 
squared  error  E[e?].  The  gradient  is  then  taken  of  the  estimate. 
This  gradient,  denoted  by  V  .  ,  is  an  estimate  of  the  true  gradient 
and  can  be  determined  to  be 

*J  =  'fej]w=Wj  "  -2«jXj  <3.3«) 

Using  this  estimate  in  place  of  the  true  gradient  in  (3.27)  gives 

Wj  +  1  =  Wj  +  2uejXj  (3.29) 

The  adaptive  filter  of  Figure  7  is  implemented  with  a  tapped 
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delay  line  as  shown  in  Figure  11.  The  input  signal  vector  is 

XJ  =  [lj   Xj-1  •  •  •  xj-n+llT  <3.30) 

The  components  of  this  vector  are  the  delayed  components  of  z.+. 

which  is  a  sample  of  the  time  series  ....  xj+a-i«  x)+A'  x1+A+1' 


When  implementing  this  algorithm  the  parameter  u  must  be 
specified.  If  made  too  large,  the  algorithm  becomes  unstable. 
If  u  is  made  too  small,  the  weight  vector  takes  too  much  time  to 
converge  to  the  optimal  weight  vector.  Although  there  are  ways 
to  determine  the  proper  u  to  insure  stability,  it  is  generally 
easier  to  use  a  tria 1-and-error  procedure  in  this  problem.  This 
was  the  manner  used  to  set  u  for  fast  convergence  in  the 
s  imul a  t  ion. 

Besides  convergence  rate  and  stability,  another  concern  in 
picking  u  is  imbedded  in  the  LMS  equation,  Eqn.  (3.27).  There 
are  two  sources  of  information  that  contribute  to  update  of  the 
weight  vector.  One  is  the  current  weight  vector  which  represents 
a  history  of  the  weight  vectors  due  to  the  regressive  structure 
of  Eqn.  (3.27).  The  second  is  the  gradient  which  is  estimated  by 
the  the  prodnct  of  the  current  filter  input  vector  and  the  error. 
The  choice  of  u  decides  what  emphasis  should  be  given  to  each. 
If  (i  is  made  small  then  the  emphasis  is  placed  on  the  current 
weight  vector.  In  this  instance  a  weight  vector  with  many 
weights  is  useful  since  it  can  retain  more  information.  If  u  is 
made  large  then  the  emphasis  is  placed  on  the  gradient  and  the 
information  contained  in  previous  weight  vectors  is  soon  lost. 
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In  this  case  a  large  number  of  weights  would  not  be  beneficial. 

For  the  band  rate  analysis  system,  both  reasons  of  choice 
have  merit.  Since  the  message  sequence  is  p s eudo- r a ndom,  the 
mean  of  the  ANC  input  will  also  be  pseudo-random.  Making  u  large 
corresponds  to  attempting  to  obtain  good  noise  cancellation.  The 
reason  u  must  be  large  is  that  convergence  should  occur  before 
the  mean  value  changes.  The  output  would  then  be  the  clean 
message  sequence  which  would  easily  indicate  the  baud  rate. 
Making  u  small  will  not  provide  good  noise  cancellation,  but 
rather,  the  weight  vector  will  contain  the  information  about  the 
sequence.  Suppose  the  mean  and  variance  of  the  changing  weight 
vectors  is  determined.  An  intuitive  argument  is  that  the  mean 
should  resemble  the  autocorrelation  function  of  the  message  dis- 
placed by  the  delay  of  the  ANC.  Once  again  the  baud  rate  could 
be  easily  determined.  As  for  the  variance,  a  relative  maximum 
should  occur  at  the  first  sample  of  the  second  bit  interval 
(again  a  displacement  results  from  the  delay)  resulting  from  the 
indecision  of  whether  the  current  bit  has  been  repeated  or  not. 
Zero-Cross  ing  Algorithm 

The  zero-crossing  algorithm  is  implemented  as  follows. 
First  a  set  is  constructed  whose  members  are  the  number  of 
samples  between  each  zero  crossing.  Next  a  histogram  is  formed 
from  the  data  of  the  set.  Since  it  is  very  possible  that  the  bit 
time  will  not  relate  to  an  integral  number  of  samples,  bins  of 
the  histogram  around  the  nearest  integral  sample  which  represents 
the  bit  time  might  also  have  entries.   For  example,  if  the  baud 
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rate  is  2  bps  and  the  sampling  rate  is  19  Hz  the  resulting  histo- 
gram might  look  like  the  one  illustrated  in  Figure  12a.  Noise 
perturbations  may  also  cause  this  result,  especially  for  low 
signal-to-noise  ratios.  Noise  may  also  introdnce  false  crossings 
as     illustrated    in    bins    number    3     and    6. 

To  negate  this  spilling  action,  the  entries  of  each  group 
are  averaged.  A  group  is  defined  as  a  collection  of  bins  not 
containing  two  or  more  empty,  adjacent  bins.  These  averages  are 
then  plotted  on  the  real  line  as  impulses  with  weights  equal  to 
the  number  of  entries  that  contributed  to  the  respective  average. 
When  applied  to  the  histogram  of  Figure  12a,  this  procedure 
results     in    the    plot    of    Figure    12b.  Under    optimum    conditions, 

the  result  will  be  a  group  of  impulses  with  the  distance  between 
adjacent  impulses  being  the  bit  time  normalized  by  the  sampling 
interval  (or  conversely,  the  sampling  rate  divided  by  the  bit 
rate).  The  leftmost  impulse  will  be  located  at  the  normalized 
bit  time  and  the  others  at  integral  multiples  of  the  baud  rate. 
When  conditions  are  not  so  good,  as  in  Figure  12b,  some  sort  of 
decisions    must    be     made. 

There  are  several  criteria  which  the  decision  must  satisfy. 
Before  listing  the  criteria,  the  definition  of  consistent  and 
inconsistent  result  should  be  stated.  A  consistent  results 
occurs  when  the  location  of  an  impulse  to  the  right  of  the  con- 
sidered impulse  is  an  integral  multiple  of  the  location  of  the 
considered     impulse,     within     some     tolerance.  An     inconsistent 

result     occurs    when     the     tolerance     is     exceeded.        The     number    of 


33 


X  X  X  X  X  X 


X  X  X  X  X 

o 

x  x  x  x  x 

01 

o 

r*- 

X 

CO 

in 

T 

X 

m 

ru 

ru 

o 

ru 

01 

■n 
CD 

CO 

•rt 

in 

•H 
V 

ro 
ru 


co 


o 

01 


E 
<B 
C 
Dl 
O 


J3 
UJ 

a 
a 

L 
Q. 


in 

01 


c 
o 


3 
13 
■H 

c 
in 


a) 
in 

3 

a 

E 


CD 

e 


in 

0) 

n 


CO 


a 
en 


o 

ITJ 

E 


in 
ID 

O) 

C 

in 
in 
o 
c. 
u 
I 

o 
c 
a> 

N 

a 


a 
a 

•H 

in 

>>. 

i-* 

m 

c 

ID 

E 
ID 
C 
O) 

o 


UJ 

3 


34 


consistencies  is  the  sum  of  the  weights  of  consistent  impulses 
while  the  tin  of  all  weight  of  inconsistent  impulses  is  the 
number  of  inconsistencies. 

The  first  criterion  is  that  the  weight  of  the  considered 
impulse  plus  the  number  of  consistencies  most  be  as  large  as  some 
number  times  the  number  of  inconsistencies.  The  second  criterion 
is  that  the  weight  of  the  considered  impulse  most  be  at  least  a 
certain  proportion  of  the  sum  of  all  impulse  weights  (which  is 
the  total  number  of  zero  crossings).  The  final  criterion  is  that 
the  number  of  consistencies  must  be  at  least  a  certain  proportion 
of  the  sum  of  weights  (number  of  zero  crossings).  Obviously 
there  is  a  tradeoff  when  setting  the  tolerance  and  the  ratios. 
The  optimal  parameters  will  determined  during  initial  runs  of  the 
simulation.  If  no  impulse  meets  the  criteria,  the  impulse  used 
for  the  bit  time  estimate  is  the  one  of  maximum  weight. 

Consider  the  example  in  Figure  12  once  again.  The  first 
impulse  meeting  the  requirements  (upon  reasonable  choice  of  the 
parameters)  is  the  one  at  9.5.  The  bit  time  estimate  is  then 
found  by  multiplying  the  location  by  the  sampling  interval.  For 
the  example  the  bit  time  estimate  would  be  0.5  seconds  since  the 
sampling  interval  is  1/19  seconds.  Finally  the  baud  rate  esti- 
mate is  found  as  the  inverse  of  the  bit  time  which  would  be  2  bps 
for  the  example. 


IV.  Simnlat  ion  De  script  ion 

The  system  modeled  in  the  previous  section  will  be  snbjected 
to  a  Monte  Carlo  simulation  to  determine  its  performance.  The 
first  step  of  the  simulation  is  to  generate  an  FSK  signal.  Next 
noise  most  be  added  subject  to  the  signal-to-noise  ratio  para- 
meter. Then  the  response  of  the  system  to  the  signa 1-p 1 ns-no ise 
input  is  observed.  This  procedure  is  repeated  many  times  so  that 
the  error  between  the  system  response  and  the  known  signal  para- 
meters can  be  averaged  in  a  sound  statistical  sense.  The  results 
will  be  the  mean  of  the  baud  rate  estimate  and  the  variance  of 
the  estimate. 
FSK  Signa 1  Modnl a  tor 

In  order  to  simulate  the  baud  rate  analysis  system,  a  repre- 
sentation for  a  received  FSK  signal  must  be  formulated.  The 
scheme  used  for  the  simulation  is  based  on  work  done  by  Humme  1  s 
and  Ratcliffe  [41.  The  input  to  the  modulator  is  a  binary 
message  sequence  having  symbols  that  are  equally  likely  and 
statistically  independent.  To  meet  these  conditions  a  psuedo- 
random  bit  sequence  is  used  with  a  bit  time  of  T.  seconds.  Each 
T\  seconds  the  modulator  puts  out  one  of  two  different  waveforms 
in  accordance  to  the  FSK  scheme.   The  two  waveforms,  which  are  of 


s12(t)  =  A  cos[2n(fc  +  AfFSK)t]  (4.1) 

The  plus  sign  is  used  for  a  1  in  the  bit  sequence  and  the  minus 
is  used  for  a  0.  The  frequency  deviation  from  the  center  fre- 
quency is  denoted  by  AfpsK.     The  FSK  signal  is  then  a  composite 
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of  the  separate  waveforms.   It  can  be  expressed  as 


>(t)  =  Y.      sn(t-nTb) 


(4.2) 


,th 


where  the  n  n    transmission  !B(t-Ij))  is  one  of  the  two  waveforms 
of  Eqn.  (4.1). 

For  this  FSK  representation,  the  signal  will  have  a  phase 
discontinuity  at  the  end  of  each  bit  interval  unless  the  two 
signaling  frequences,  fc+fpSK  and  tg-fpsZ'  contain  an  integral 
number  of  whole  cycles  per  bit  interval.  This  discontinuity  will 
occur  even  when  one  bit  is  repeated.  Generally  FSK  is  generated 
using  two  oscillators,  one  operating  at  each  of  the  different 
signalling  frequencies.  The  output  is  produced  with  a  switch 
that  connects  to  either  one  oscillator  or  the  other.  The  switch 
is  contolled  by  the  bit  sequence.  In  this  manner,  the  phase  of 
each  oscillator  output  will  be  continuous  when  a  bit  is  repeated. 
This  form  was  used  for  the  simulation  since  it  is  a  commonly  used 
scheme.   The  equation  used  to  implement  this  FSK  scheme  is 


l(t)  =   Z  *B(t){u(t-nTb)[l-u(t-nTb-Tb)]J 


(4.3) 


This  time  sn(t)  is  not  limited  to  a  duration  of  Tb  seconds, 
instead  its  duration  will  be  for  the  entire  signalling  interval. 
The  unit  step  functions,  represented  by  u(t),  create  a  window 
from    time    nTb    to     (n+l)Tb     »nd     thus     implement     the     switching. 

To  be  useful  the  resulting  bandpass  FSK  signal,  s(t),  should 
be  expressed  in  its  low-pass  representation,  s(t).  To  begin  this 
development    each  waveform    s     (t)    has    a    complex    envelope    that    is 
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easily  shown  to  be 

+j2"AfFSKt 
sl,2(t)  =  Ae  (*.*) 

Then  by  definition  and  substitution  each  waveform  can  be  expressed 

in  terns  of  its  complex  envelope  as 

j2nfct 
sn(t)  =  Rets(t)e       }  (4.5a) 

+j2nAfFSKt  j2nfct 
-  Re{Ae  e       )        (4.5b) 

Substituting  the  definition  result  into  the  FSK  repre sent  a  ion  gives 

s(t)  -  2.  Ref»„(t)e       )u(t-nTb) [l-u(t-nTb-Tb)]     (4.7) 
n  =  -<» 

Moving  the  summation  inside  the  equaton  yields 

J2nfct   - 
s(t)  =Re{e         2.  sn(t)u(t-nTb) [l-u(t-nTb-Tb)])     (4.8) 
n=-w 

Not  the  complex  envelope  of  s(t)  is  seen  to  be 

■ 

*(t>  =     I  sn(t)n(t-nTb) [l-u(t-nTb-Tb)]  (4.9) 

n=-<» 

Then  substituting   Eqn.  (4.4)  into  this  result  yields  a  low-pass 

form  of  an  FSK  signal  that  can  easily  be  implemented  in  the 

simulation.   The  equation  for  implementation  is 

-    +2nAfFSKt 
s(t)  =     2,  Ae  u(t-nTb)  [l-n(t-nTb-Tb)]     (4.10) 

n=-"> 

Obviously  the  limits  of  the  summation  become  finite  for  a  finite 

length  message.   Once  again  the  +  distinguishes  between  either  a 

0  or  a  1  being  sent. 

No  ise  Mode  1 

The  FSK  signal  generated  for  simulation  purposes  by  the 
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expression  above  will  have  noise  added  to  it  before  entering  the 
system.  Since  the  object  of  the  investigation  is  to  determine 
how  well  one  method  of  band  rate  analysis  performs  at  low  signal- 
to-noise  ratios  (SNR),  the  manner  of  setting  the  SNR  will  be 
discussed.  The  formal  definition  of  the  SNR  is  the  ratio  of 
signal  power  to  noise  power  as  measnred  after  the  pre-filter. 
Note  that  SNR  is  not  meaningful  until  the  pre-filter  bandwidth  is 
specified.  The  noise  is  assumed  to  be  from  a  white  Gaussian 
random  process. 

Suppose  a  random  sample,  i(m),  of  length  N  is  drawn  from  a 
statistically  independent,  Gaussian  random  process  with  zero 
mean.  The  average  power  in  the  sample,  Px>  is  a2,  the  variance 
of  each  sample.   This  can  be  expressed  as 

N-l 
=  0 


5  I    =* 

N  „4rn 


(m) 


(4.11) 


where  the  overbar  denotes  mean  value.   Parseval's  theorem  for  the 
discrete  Fourier  transform  [5]  which  is  defined  by 


N-l  -j2irmn/N 

X(n)     =   V   x(m)e 
m=0 


(4.12) 


•2-^5>Ha 


(4.13) 


states    that    (4.11)     can    also    be    expressed    a: 

N-l 

X 

It  can  be  shown  that  the  mean  value  of  |x(n)|2  is  independent  of 
n  and  thus  is  No2.  Suppose  the  noise  sample  is  passed  through  a 
low-pass  filter  as  illustrated  below. 


x  (m) 

h  fa) 

H(n) 

y  (m) 

X(n) 

Y(n) 
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The  power  response  of  the  filter  is 


|H(nf0) 


k=0,l N/2 


JH((n-N)f0)|  2    k=N/2+l,....N-l 


(4.14) 


where  fg  is  the  fundamental  frequency  of  the  time  sample.   Then 


the  output  power,  P  ,  is 


„2  N-l 

—  r  lHn 2 


(4.15) 


To     simplify     this     expression,      the     equivalent     noise     bandwidth     for 
the    filter   can  be    defined    as 

N-l 


2B_    = 


51      LIHn|2f0 
"0    n=0 


N-l 


^  n-j. 


(4.16a) 


(4.16b) 


This  is  actually  an  approximation  which  grows  better  as  N  grows 
larger.  Note  that  f()-l/(T1N)  where  T,  is  the  sampling  interval. 
Substituting  this  result  into  (4.15)  and  solving  for  a2    gives 


2H6TsBn 


(4.17) 


These     results     are     for     the     low-pass     case.        Since    noise     enve- 
lopes   used    for    the    bandpass    case    are    complex    and    have,    in    gene- 
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ral,  both  real  and  imaginary  parts,  it  will  be  necessary  to 
extend  these  concepts  to  the  low-pass  equivalent  model  as  shown 
below. 


For    the    low-pass  model,    the      powers  P~   and  P      are   related  by 


p      =    T    P- 
y         2      y 


(4.18) 

This  can  be  shown  by  taking  the  mean  squared  value  of  the  expres- 
sion relating  y(t)  and  y(t). 

Now  each  sample  of  the  complex  noise  envelope  x(m)  requires 
two  numbers  from  the  random  process,  one  for  the  real  part  and 

one  for  the  imaginary  part.   The  object  is  to  find  a  relation 

2         2 
between  aT    and  o^  to  the  output  power  P   of  the  bandpass  process 

y(t).   The  case  of  interest  is  where  *r(t)  and  x.(t)  are  indepen- 
dent processes  with  identical  statistical  properties  so  that  o^  = 

2     2 
a j  =  a    .   Since 

y(t)  =  x(t)*E(t)  (4.19a) 

=  [5r(t)+jxi(t)]*[nr(t)+jE1(t)]         (4.19b) 

(*  denotes  convolution)  it  is  obvious  that 


7r(t)  =  xr(t)»Kr(t)  -  SjUJ'E.U) 


and 


yi(t)  -  51(t)*Er(t)  +  5r(t)»Ei(t) 

In  many  cases  E(t)  will  be  a  real  function  so  that 


(4.20a) 


(4.20b) 
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yr(t)  =  xr(t)*Er(t)  (4.21.) 

y^t)  =  ii(t)*Er(t)  (4.21b) 

Recall  that  the  low-pass  procedure  derived  earlier  can  be  applied 

to  these  expressions.  Hence 


»l      -  P~ 


*r     yr  2TsBnH8 


(4.22a) 


,2      m   p_   _* 
*i     yi  2TsBnH0 


(4.22b) 


Note  that  Bn  is  the  equivalent  low-pass  bandwidth  which  is  one- 
half  the  bandpass  bandwidth,  i.e.  Bn  =  BRp/2. 

Then  since  the  variances  have  been  assumed  equal,  it  follows 
from 

P_  =  P_  +  P„  =  2P  (4.23) 

y    yr   yi    y 

that 

(4.24) 


P~  =  P™  =  P 

yr   yi   y 

The  conclusion  is  to  set  or2  =  cr2  =  »2  where 

xr     ii 


2T,B^ 


(4.25) 


The  final  step  is  to  relate  the  noise  power  P  ,  the  noise 
variance  a2,  and  the  SNE.  Since  the  FSK  signal  is  a  sinusoid  of 
amplitude  A  at  one  of  two  frequencies,  the  signal  power  is  A2/2. 
Therefore 


SNR  = 


A2/ 2 


(4.26) 
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Or 


r    -  til 

t  SNE 


(4.27) 


The  final  result  is  obtained  by  substituting  this  expression  for 


P   into  (4.24).   Thus,  the  noise  should  have  zero  mean  and  var- 
iance defined  by 


Az/2 


TsBn(SNR) 


(4.28) 


Or  since  Bn  =  BRp/2, 


TSBRF(SNR) 


(4.29) 


V.   Prel iminarv  Observa t ions  and  Es  t ima tor  Definitions 

The  object  of  the  investigation  is  to  determine  the  baud 
rate  of  an  FSK  signal  using  the  adaptive  noise  canceller  (ANC). 
In  order  to  do  so,  first  a  decision  has  to  be  made  as  to  how  the 
baud  rate  can  be  obtained  from  the  information  provided  by  the 
ANC.  Since  the  ANC  is  a  non-linear  device  (linearity  only  occurs 
once  convergence  of  the  weight  vector  is  reached  and  maintained), 
this  decision  should  be  made  after  the  outputs  of  the  ANC  for 
different  parameters  and  situations  are  observed.  The  observa- 
tions will  be  presented  in  this  section  along  with  the  decision 
of  how  to  obtain  the  baud  rate  from  the  ANC  information. 

To  generate  the  signal,  the  baud  rate  for  the  FSK  generation 
is  arbitrarily  set  to  300  Hz.  The  frequency  deviation  will  be 
initially  set  to  300  Hz.  Generally  for  FSK  systems  the  frequency 
deviation  must  be  as  large  as  the  band  rate  to  insure  quality 
reception.  As  deviation  increases,  performance  does  also. 
Therefore  this  represents  the  worst  case  situation  for  a  given 
SNR.  The  signal  is  sampled  at  a  rate  of  9600  Hz  and  1024  samples 
are  recorded.  This  means  that  32  bits  will  be  seen  in  the 
sampled  sequence  and  each  bit  will  consist  of  32  samples.  As 
discussed  earlier,  the  system  is  designed  to  operate  for  FSK 
signals  with  signal  bandwidth  of  6000  Hz  or  less  and  baud  rates 
of  1000  Hz  or  less.  Consequently  the  pre-filter  bandwidth  is  set 
to  6000  Hz  and  the  post-filter  bandwidth  is  set  to  1000  Hz.  In 
order  that  the  DLD  operates  in  the  linear  region  for  this  band- 
width, the  parameters  of  the  DLD  are  T  =  13.854  msec    and  n  = 
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266  (see  Egns.  (2.3)  and  (2.6)).  As  for  the  limiter  and  post- 
filter,  they  will  he  included  in  the  demodulator  system.  The 
reason  for  doing  so  will  be  substantiated  shortly.  These 
parameters  are  held  constant  throughout  the  observations  unless 
otherwise  noted.  Furthermore,  unless  otherwise  specified,  the 
signal-to-noise  ratio  is  set  to  8  dB. 
DLD  Output 

It  was  previously  argued  that  the  post-filter  should  be 
omitted  when  the  ANC  approach  for  band  rate  determination  is  used 
(see  the  ANC  discussion  in  section  II).  Figure  13  suggests  that 
the  post-filter  should  be  included.  The  four  plots  in  Figure  13 
are  the  DLD  output  for  the  aforementioned  parameter  set  except  as 
noted  on  each  plot.  The  noise  sequence  added  to  the  FSK  signal 
is  the  same  for  each  case.  These  plots  investigate  the  effect  of 
the  post-filter  and  the  limiter  upon  the  output.  Even  for 
modest  bandwidth  of  1000  Hz,  the  post-filter  clearly  improves  the 
output.  As  for  the  limiter,  when  the  post-filter  is  included,  it 
too  is  desirable.  When  the  post-filter  is  omitted  the  use  of  the 
limiter  is  questionable.  Based  on  these  observations,  the  post- 
filter  and  limiter  should  be  included  in  the  system. 

The  output  of  the  DLD  (post-filter  and  limiter  included)  for 
four  cases  is  plotted  in  Figure  14.  The  four  cases  represent  a 
noise-free  FSK  signal,  and  signals  with  SNR  of  12  dB,  8  dB  and  4 
dB.  As  should  be  expected,  the  signal  is  degraded  as  the  signal- 
to-noise  ratio  decreases.  For  SNR  -  12  dB  the  message  sequence 
can  be  easily  identified.   When  the  SNR  drops  to  8  dB  the  dis- 
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FIGURE  13:   Effect  of  post-filter  and  limiter  on  DLD 
output . 
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FIGURE  14:   DLD  output  for  signals  of  small  deviation 
and  varying  SNH. 
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tinction  between  the  ones  and  zeroes  is  less  obvious.  At  SNR  =  4 
dB,  the  bits  are  indistinguishable. 

Figure  IS  shows  the  OLD  output  for  the  same  signal-to-noise 
ratios  but  with  the  frequency  deviation  set  to  2000  Hz.  This  is 
more  likely  the  case  for  practical  systems  than  the  minimum 
deviation  as  shown  in  Figure  14.  The  ones  and  zeroes  are  now 
clearly  discernible  at  SNR  =  8  dB  but  only  vaguely  discernible 
when  the  SNR  is  4  dB. 
ANC  Invest  i«a t ion 

The  input  of  the  ANC,  implemented  as  in  Figure  7,  consists 
of  the  output  sequence  from  the  DLD.  The  output  of  the  ANC  can 
be  thought  of  as  consisting  of  three  sequences,  the  error 
sequence,  the  adaptive  filter  output  sequence,  and  the  weight 
vector.  The  object  is  to  fix  the  parameters  of  the  ANC  such  that 
the  baud  rate  can  be  estimated  from  one  of  these  three  sequences 
more  accurately  than  it  can  be  from  the  input  sequence. 

The  method  postulated  in  the  previous  section  is  that  the 
message  and  the  noise  can  be  separated  by  the  ANC.  This  relies 
on  the  earlier  made  observation  that  the  noise  portion  of  the 
input  sequence  when  delayed  8  or  more  samples  (with  the  post- 
filter  and  limiter)  is  uncorrelated  with  the  noise  of  the  non- 
delayed  input  but  that  the  message  portion  is  still  correlated. 

Recall  the  LMS  equation  used  to  implement  the  method.   It  is 


f  j  +  1  '-  Wj 


■*. 


(5.1) 


where  the  estimate  for  the  gradient  is 


Vj  "  -2eJXj 


(5.2) 
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In  order  for  the  cancelling  (actually  predicting  in  this  situa- 
tion) to  perforin  well  the  gradient  constant  u  must  be  large  so 
that  the  weight  vector  converges  rapidly.  Also,  since  the  mes- 
sage is  ps  ue  do- random,  the  ANC  will  not  be  able  to  predict  past 
one  bit.  Therefore,  convergence  needs  to  occnr  within  one  bit 
time.  Thus,  it  is  best  to  pick  u  as  large  as  possible  so  that 
stability  is  maintained.  Furthermore,  the  bias  weight  shonld 
also  be  included  since  it  allows  the  weight  vector  to  respond 
quickly  to  changes  in  the  input  level.  There  are  methods  to 
determine  the  proper  u  [1]  but  the  choice  is  easier  made  with 
tr ia 1-and-error  procedure.  For  a  weight  vector  composed  of  128 
weights  u  can  be  as  large  as  0.001  without  chancing  instability 
(maximum  u  is  dependent  of  the  weight  vector  length  in  an  inverse 
relat  ion) . 

The  error  sequence  can  be  inspected  to  see  how  well  the  ANC 
performed  the  cancelling.  Since  the  message  is  the  portion  being 
cancelled,  under  optimal  cancelling  the  error  would  be  the  noise 
portion  of  the  input.  Therefore,  if  the  error  seems  to  have  a 
stationary  mean,  the  ANC  provided  good  cancellation.  This  is 
synonymous  to  good  prediction,  i.e.  the  adaptive  filter  output 
should  be  the  message  portion  of  the  input. 

Recall  that  for  the  sampling  frequency  of  9600  Hz  and  a  baud 
rate  of  300  Hz  there  will  be  32  samples  per  bit.  However  the 
system  is  designed  for  band  rates  up  to  1000  Hz  which  would  have 
9.6  samples  per  bit.  Besides  the  problem  of  attaining  conver- 
gence in  the  time  of  9.6  samples  another  problem  exists.   Since  a 
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signal  with  a  1000  Hz  baud  rate  will  be  nncorrelated  after  the 
tenth  sample  for  the  given  sampling  freqneney,  the  delay  of  8 
samples  might  decorrelate  the  message  so  much  that  cancelling  can 
not  be  performed.  A  closer  examination  of  the  results  should  be 
made  in  setting  this  delay.  Figures  16a  and  16b  provide  addi- 
tional insight.  Figure  16a  consists  of  the  three  ouput  sequences 
of  the  ANC  for  the  case  of  A  =  8.  The  other  parameters  are  the 
standard  set  with  u  =  0.001  and  N  ■  128.  In  Figure  16b  the  delay 
has  been  changed  to  A  =  1  while  all  other  parameters  are  the  same 
as  in  Figure  16a. 

The  main  conclusion  to  draw  from  Figures  16a  and  16b  is  that 
the  output  sequences  for  the  different  choices  of  delay  are  very 
similar.  This  means  that  the  delay  can  be  set  to  one  sample  and 
the  results  should  not  be  degraded.  By  being  able  to  do  this, 
the  problems  of  convergence  associated  with  the  larger  baud  rates 
should  not  be  as  severe.  One  should  recall  that  decreasing  the 
delay  increases  the  correlation  of  the  noise  and  that  this  should 
make  matters  worse.  However  the  correlation  of  the  message  is 
also  increased.  Referring  back  to  Figure  6  for  the  post-filter 
and  limiter  case,  the  slope  of  the  message  autocorrelation  is 
larger  than  the  slope  of  the  noise  autocorrelation.  Thus,  the 
correlation  of  the  message  increases  faster  than  that  of  the 
noise  as  the  delay  is  reduced. 

Referring  again  to  Figures  16a  and  16b,  note  how  the  error 
sequence  dictates  the  degree  of  convergence.  In  both  cases  the 
error  initially  has  distinct  mean  values.   After  the  weights  are 
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given  some  time  to  adjust  from  their  all-zero  initial  value,  the 
mean  of  the  error  does  not  vary  as  much.  Although  the  two  error 
sequences  are  nearly  the  same,  the  end  of  the  error  sequence  for 
the  delay  of  8  seems  to  fluctuate  slightly  more  than  when  the 
delay  is  1.  This  suggests  convergence  is  more  complete  for  the 
delay  of  1.  However,  the  ADF  output  is  a  bit  more  ragged  for  the 
smaller  delay  than  for  the  larger  delay.  Still  this  is  a  consi- 
derable improvement  over  the  original  input  which  is  depicted  in 
the  top  plot  of  Figure  13.  Another  point  is  that  the  weight 
vectors  appear  to  be  distinctly  different.  (Weight  vector  in 
this  context  implies  the  mean  and  variance  of  the  individual 
weight  sequences  of  the  iterative  update.)  A  second  look  might 
suggest  that  the  weight  sequence  for  the  larger  delay  is  just  the 
weight  sequence  of  the  smaller  delay  advanced  by  the  difference 
in  the  delay. 

The  question  impending  is  how  to  extract  the  baud  rate  from 
this  data.  As  previously  noted,  the  ADF  output  should  be  the 
message  portion  of  the  input  upon  ideal  cancellation.  Although 
the  cancellation  is  not  ideal,  the  result  is  still  good.  By 
inpection  of  the  ADF  output  and  the  corresponding  input  (Fignres 
13  and  16b)  it  can  be  concluded  that  the  zero- crossing  analysis 
would  yield  better  results  when  performed  on  the  ADF  output  than 
when  performed  on  the  ANC  input. 

The  weight  vector  can  also  be  used  to  determine  the  baud 
rate.  Since  the  message  sequence  is  uncorrelated  from  the  first 
sample  of  the  second  bit  onward,  one  should  suspect  some  special 
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behavior  at  this  point  in  the  mean  weight  vector.  By  knowing  the 
number  of  samples  per  bit  to  be  32,  the  weight  number  to  look  at 
would  be  31  (since  the  delay  is  one).  Figure  17  gives  the  mean 
and  variance  of  the  first  40  weights  of  two  other  trials  produced 
by  different  runs  with  different  noise  sequences.  The  trend  for 
the  mean  of  this  specific  weight  is  for  it  to  be  a  relative 
minimum.  To  distinguish  it  from  the  other  relative  minima  of  the 
mean  weight  vector,  two  other  observations  can  be  made.  One, 
this  minimum  is  the  first  relative  minimum  less  than  zero.  Two, 
this  miminum  coincides  with  a  maximum  in  the  variance  of  the 
weight  vector. 

These  criteria  were  not  conceived  through  analysis  but  were 
noticed  by  inspecting  many  weight  vectors  for  many  cases.  Some 
of  the  observations  of  this  investigation  should  be  mentioned. 
Recall  the  LMS  equation  (Eqn.  (5.1)).  If  the  information  is  to 
come  from  the  weight  vector  it  seems  that  the  update  procedure 
should  emphasize  the  weight  vector,  i.e.  make  n  small.  Figure  18 
presents  two  weight  vector  sequences  for  small  u.  They  were 
calculated  for  the  same  inputs  as  those  of  Figure  17.  In  both 
comparisons  between  large  and  small  u  the  weight  vector  has  the 
nearly  the  same  shape.  The  only  basic  difference  is  the  magni- 
tude for  the  two  cases.  It  seems  that  selection  of  ji  for  the 
weight  vector  baud  rate  estimate,  so  long  as  stability  is 
sustained,  is  not  critical.  This  was  confirmed  by  preliminary 
runs  of  the  simulation. 

Another  observation  which  was  previously  mentioned  is  that 
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the  weight  vector  characteristics  do  not  change  for  different 
delays.  The  delay  only  influences  the  amount  the  weight  vector 
is  shifted  with  respect  to  the  weight  vector  of  the  non-delayed 
condition.  Supporting  this  observation  is  Figure  19  which  pre- 
sents three  results  for  the  same  input.  The  only  parameter  which 
differs  is  the  delay  chosen. 

One  interesting  observation  is  that  when  the  post-filter  of 
the  DLD  is  omitted  the  weight  vector  mean  and  variance  are  more 
spurious.  This  makes  identification  of  the  minimum  and  maximum 
easier  which  should  lead  to  better  estimates.  However  prelimi- 
nary simulation  runs  indicate  better  estimates  are  obtained  with 
the  post-filter  included.  Evidently  the  desirable  noise  reduc- 
tion provided  by  the  post-filter  outweighs  the  undesirable 
smoothing  effect  of  the  filter  on  the  weight  vector. 

The  weight  vector  estimate  described  herein  is  to  make 
estimates  from  the  mean  and  variance  of  the  individual  weight 
vectors  of  each  update  iteration.  Keep  in  mind  that  the  estimate 
result  will  be  expressed  in  number  of  samples  therefore  intro- 
ducing a  bias  if  the  bit  time  does  not  equate  to  an  integral 
number  of  sampling  intervals.  Another  approach  would  be  to  apply 
the  estimate  procedure  to  each  separate  weight  vector  then 
average  the  results.  This  way  the  estimate  is  not  limited  to  an 
integral  number  of  sampling  intervals.  However  preliminary  runs 
once  again  suggest  that  better  results  are  produced  when  the 
estimates  are  made  from  the  mean  and  variance  of  the  weight 
vectors. 
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FIGURE  19:   Effect  on  weight  vector  of  varying  delay. 
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Another  attempt  at  estimating  the  band  rate,  which  was  not 
too  successful,  employed  the  similarities  of  the  weight  vector 
mean  and  variance  with  the  autocorrelation  of  the  message. 
Recall  that  a  DFT  of  the  autocorrelation  of  a  sequence  provides  a 
power  spectral  density  (PSD)  estimate  of  the  sequence.  It  is 
known  that  the  PSD  of  the  message  will  contain  nulls  at  integral 
multiples  of  the  baud  rate.  The  strategy  is  to  construct  a 
pseudo-autocorrelation  sequence  from  either  the  mean  or  variance 
of  the  weight  vector,  perform  a  DFT,  and  look  for  the  identifying 
nulls.  Note  that  folding  the  weight  vector  mean  about  the  zero 
axis  gives  a  rough  looking  autocorrelation  function.  For  the 
variance,  the  triangular  shape  is  also  present  but  with  a  base  of 
one  bit  time  instead  of  the  expected  two  bit  times.  Once  again 
preliminary  runs  discredited  the  credibility  of  this  method. 

Recall  that  one  of  the  problems  with  this  receiver  is  in 
setting  the  bandwidth  of  the  post-filter  of  the  DLD.  Obviously 
if  the  baud  rate  is  300  Hz  and  the  bandwidth  is  1000  Hz  the  baud 
rite  estimates  will  not  be  as  accurate  as  if  the  bandwidth  were 
set  to  300  Hz.  Thus  it  makes  sense  to  try  to  determine  the 
bandwidth  in  order  to  improve  results.  This  is  synonymous  to 
determining  the  baud  rate.  Therefore  any  of  the  previously 
mentioned  methods  of  estimating  the  baud  rate  could  be  used  to 
set  the  bandwidth  of  a  second  low-pass  filter.  This  narrower 
filtering  would  permit  better  results  to  be  achieved.  This 
modified  system  is  shown  in  Figure  20.  The  intermediate  filter  is 
a  Butterworth  filter  with  4  poles. 
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The  bandwidth  determination  of  the  filter  is  performed  nsing 
a  DFX  analysis  of  the  weight  vector  instead  of  one  of  the  baud 
rate  estimates  previously  discussed.  One  reason  for  this  choice 
is  that  the  baud  rate  estimates  are  either  fairly  good  or  very 
bad.  For  example  consider  the  estimate  that  is  based  on  the 
minimum  of  the  mean  weight  vector  coinciding  with  the  maximum  of 
the  variance  of  the  weight  vector.  If  the  points  coincide  the 
estimate  will  be  very  accurate,  if  they  do  not  the  estimate  will 
depend  on  if  the  condition  is  met  somewhere  else  in  the  sequence. 

Another  reason  for  using  the  DFT  technique  is  that  it  is 
very  stable.  Consider  the  ADF  filter.  Its  input  is  the  message 
plus  noise  sequence  and  its  output  is  the  forward  shifted  message 
sequence.  The  weight  vector  is  the  impulse  function  which  accom- 
plishes this  to  the  best  degree.  As  such  the  DFT  of  the  weight 
vector  provides  the  transfer  function  of  this  impulse  function. 
If  the  task  is  performed  well  the  transform  should  indicate  the 
bandwidth.  This  bandwidth  estimator  might  not  provide  the  exact 
baud  rate  every  time  but  should  be  very  consistent  and  should  not 
fluctuate  wildly. 

The  bandwidth  estimate  used  is  twice  the  3  dB  bandwidth  of 
the  transform  of  the  weight  vector.  This  choice  was  arbitrarily 
made  and  found  to  give  good  results.  The  justification  for  the 
choice  is  based  on  the  transform  of  a  message  sequence.  For  this 
transform  the  magnitude  drops  about  3  dB  at  one  half  the  baud 
rate  which  coincides  with  the  arbitrary  choice.  It  was  found 
that  the  bandwidth  estimate  is  more  accurate  when  u  is  small. 
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This  is  hard  to  explain  since  it  was  previously  shown  that  the 
weight  vector  shape  was  more  or  less  independent  of  the  choice  of 
u.  It  could  be  possible  that  as  the  SNR  gets  lower  the  weight 
vector  values  become  dependent  upon  the  choice  of  u.  Further 
simulation  will  be  needed  to  resolve  this  issue. 

Figure  21  shows  the  input  to  the  ANC  when  it  is  taken 
directly  from  the  output  of  the  DLD.  This  is  the  top  plot.  The 
bottom  plot  shows  the  input  to  the  ANC  after  the  intermediate 
filtering  has  been  performed.  Figure  22  shows  the  ANC  output 
sequences  when  the  input  was  the  sequence  of  the  bottom  plot  of 
Figure  21.  It  was  previously  indicated  that  the  ADF  output 
should  be  a  clean  version  of  the  message,  which  it  is,  and  that 
the  baud  rate  estimate  could  be  made  with  a  zero-crossing  analy- 
sis. Unfortunately,  the  ADF  output  appears  to  be  too  heavily 
filtered;  therefore,  the  zero  crossings  might  be  displaced  due  to 
the  smoothing.  The  best  results  should  come  from  using  the  zero- 
crossing  analysis  on  the  ANC  input  after  it  has  been  intermedi- 
ately filtered.  With  intermediate  filtering,  the  weight  vector 
mean  and  variance  of  Figure  22  are  beginning  to  look  more  like 
the  autocorrelation  function.  Again  the  smoothing  is  undesirable 
in  the  sense  that  it  will  probably  be  harder  to  determine  the 
baud  rate  from  the  minimum  of  the  mean  and  maximum  of  the 
variance . 
Est  imator  Def init  ions 

Based  on  the  preceding  observations,   there  are  several 
classes  of  estimators  worth  considering.   They  are  listed  and 
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described  in  Table  1.  Notice  should  be  taken  that  the  parameters 
of  the  two  separate  ANC's  during  the  simulation  are  the  same  with 
the  exception  of  u.  For  the  intermediate  ANC  (the  one  making  the 
bandwidth  estimate)  u  is  small.  A  large  valne  of  u  is  used  for 
the  final  ANC.  Therefore  any  estimates  obtained  from  the  inter- 
mediate ANC  (i.e.,  any  estimates  made  independent  of  the  second 
ANC)  should  be  associated  with  small  u.  The  estimates  made  with 
the  final  ANC  are  accordingly  associated  with  large  u. 
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ACRONYM 


MMI* 
MMF* 


M2MI 
H2MF 


AMI 

AHF 


ML  TO  I 
MLTOF 


MZCI 
MZCF 


XZCF 


EZCI 
EZCF 


TZCF 


MAI 

MAF 


VAI 
VAF 


MEANING 


Min-Mai 


Min-2- 
Max 


Absolute 

minimum 


Min<0 


Mean<0 


ZC  of  ANC 

input 

ZC  of  ANC 
error 

ZC  of  ADF 
output 

mean/auto- 
corr . 

var. /auto- 
corr . 


DEFINITION 


first  minimum  of   the  weight 
vector  mean  that  coincides  with 
a  maximum  of  the  vector  variance 

first  minimum  such  that  a  maximum 
is  within  two  weights  of  the  minimum 

absolute  minimum  of  the  weight 
vector  mean' 

first  minimum  of  the  weight 
vector  mean  that  is  less  than  0 

first  weight  of  the  mean  that  is 
less  than  0 

zero-crossing  analysis  of  ANC  input 
zero-crossing  analysis  of  ANC  error 
zero-crossing  analysis  of  ANC  output 


similarity  of  mean  weight  vector  to 
the  message  autocorrelation 

similarity  of  variance  of  vector  to 
the  message  autocorrelation 


The  I  indicates  the  estimate  was  made  from  the  intermediate 
ANC  and  the  F  indicates  it  was  made  from  the  final  ANC. 


TABLE  1:  Estimator  Definition 


VI.   Re  sal  t  s  jind  Conclns  ions 

All  results  presented  in  this  section  were  generated  with 
simulation  runs  of  100  trials  each.  The  estimate  of  each  trial 
was  normalized  as  such  to  give  a  fractional  error 

baud  rate  estimate  -  baud  rate 


baud  rate 


(6.1) 


In  the  simulation  the  baud  rate  was  actually  estimated  in  number 
of  sampling  intervals.  The  normalization  was  done  taking  this 
into  account  so  the  result  is  the  same  as  if  the  estimate  was 
made  in  bps.  The  bandwidth  estimates,  if  normalized,  were  also 
done  in  this  manner. 

The  results  are  presented  as  the  mean  error,  e,  and  the 
root-mean-square  (rms)  error,  erms.  They  are  the  sample  mean  and 
sample  standard  deviation  of  the  error.   That  is 


N 


and 


s  £«•»  - 5) 


(6.2) 


(6.3) 


where  e.  is  the  error  of  the  i  iteration  and  N  is  the  number  of 
trials.  For  N  =  100  the  sample  mean  and  standard  deviation 
should  be  a  good  estimate  of  the  distribution  mean  and  standard 
deviation.  The  results  should  be  judged  according  to  both  the 
parameters  since  it  is  possible  for  an  estimator  to  give  a  con- 
sistent but  extremely  biased  estimate  or  to  give  a  highly  fluc- 
tuating estimate  whose  mean  is  within  tolerable  limits. 
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Eva lnat ion  of  Est  ima tors 

The  results  using  the  estimators  defined  in  the  previous 
section  are  presented  by  estimator  class  in  Figures  23a-d.  For 
the  weight  vector  estimators  without  intermediate  filtering  shown 
in  Figure  23a,  only  two  estimators  give  good  resnlts  down  to  8 
dB.  They  are  the  HLTOI  and  the  HZCI  estimators.  The  others  are 
too  inconsistent  to  be  of  any  good.  When  the  intermediate  fil- 
tering is  incorporated,  the  overall  results  for  this  group  of 
estimators  do  not  get  much  better.  This  time  the  good  estima- 
tors, as  seen  in  Figure  23b,  are  the  AMF  and  the  MLTOF  estima- 
tors. The  filtering  has  introduced  a  distinct  bias  in  the  MZCF 
estimator  and  the  rms  error  for  the  others  are  intolerable. 

In  Figure  23c  the  results  for  the  zero  crossing  analysis 
using  the  ANC  data  are  presented.  Obviously,  any  estimator 
should  work  fairly  well  for  larger  SNR,  e.g.,  12  dB.  Thus,  all 
the  estimators  except  the  XZCF  estimator  should  be  ruled  out. 
The  XZCF  estimator  gives  the  best  results  for  SNR  of  12  and  16  dB 
of  any  estimator  considered  thns  far.  Its  performance  for  SNR  ■ 
8  dB  is  slightly  worse  than  the  MLTOI  estimator. 

Figure  23d  presents  the  pseudo-correlation  sub-class  of  the 
weight  vector  class.  The  rms  error  for  this  group  is  the  best  of 
any  group  considered.  However,  the  mean  error  discredits  every 
one  of  the  estimators.  Note  the  VAI  estimator.  It  is  very 
consistent  and  has  a  steady  mean  error.  It  would  be  very  useful 
if  the  bias  could  be  directly  linked  to  the  baud  rate.  Until 
this  connection  can  be  made  it  should  be  considered  as  a  poor 
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FIGURE  23a:   Evaluation  of  estimators. 
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FIGURE  23b:   Evaluation  of  Estimators. 
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FIGURE  23c:   Evaluation  of  estimators. 
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FIGURE  23d:   Evaluation  of  estimators. 
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e  *t ima tor. 
Demodulator    Evaluation 

The  limiter  and  the  post-filter  were  included  in  the  FSK 
demodulator  when  the  previous  estimator  results  were  obtained. 
It  was  previously  argned  that  this  combination  would  provide  the 
best  results.  In  an  effort  to  substantiate  this,  the  different 
combinations  of  post— filter  and  limiter  were  used  in  obtaining 
the  estimates.  Figure  24a  shows  the  results  for  the  XZCF  estima- 
tor. Clearly  the  post-filter  should  be  included.  The  presence 
of    the    limiter    also    seems    to    improve     the    results. 

Another  way  to  evaluate  the  post-filter  and  limiter  is  to 
observe  the  accuracy  of  the  bandwidth  estimate  when  the  post- 
filter  and  limiter  combination  is  varied.  Figure  24b  shows  these 
accuracy  comparisons.  Obviously  the  post-filter  and  limiter 
should    be     included. 

Since  the  XZCF  estimator  is  critically  dependent  on  the 
noise  reduction  performed  by  the  intermediate  filter,  it  is  also 
dependent  on  the  bandwidth  estimate.  The  previous  comparisons 
will  be  made  for  the  HLTOI  estimator.  This  estimator  does  not 
depend  on  the  intermediate  filtering.  Figure  24c  gives  the  error 
results  for  this  estimator  under  the  same  situations.  Once  again 
the  resnlts  support  inclusion  of  the  post-filter  and  limiter. 
ANC    Parameter    Evalnat  ion 

Another  parameter  set  used  on  the  basis  of  assumption  and 
intuitive  feeling  is  the  ANC  parameter  set.  It  was  previously 
argued    that    the    delay    should    be    1    and    the    gradient    constant,     u 


74 


SNR  (dB) 


POST-FILTER.  LIMITER 
POST-FILTER  NO  LIMITER 
NO  POST-FILTER.  LIMITEP, 
NO  POST-FILTER  NO  LIMITEP. 


X2CF  ESTIMATOR 
INTERMEDIATE  FILTERING 


R  - 

if 

f. 


300  Hz 
■  300  Hz 
•  9600  Hz 


DELAY  -  1 

MU  -  0.1E-8 

N  -  128 

BIAS 


FIGURE  24a:   Evaluation  of  FSK  demodulator. 
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FIGURE  24b:   Evaluation  of  FSK  demodulator, 
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FIGURE  24c:   Evaluation  of  FSK  demodulator. 
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should  be  large.  Table  2  presents  the  errors  for  a  variety  of 
combinations  of  these  two  painters  for  a  SNR  of  8  dB.  The 
frequency  deviation  is  300  Hz  and  the  band  rate  is  300  bps.  The 
weight  vector  length  in  all  cases  is  128.  This  night  be  unneces- 
sarily long  but  will  not  be  detrimental.  The  weight  vector  bias 
has  been  included  since  it,  neither,  will  degrade  the  resnlts. 
The  estimates  were  made  with  the  ANC  following  the  intermediate 
filter.  The  ANC  making  the  bandwidth  estimate  had  u  set  to  0.1E- 
8  and  the  delay  as  specified  in  the  table.  The  resulting  esti- 
mate of  filter  bandwidth  is  specified  in  the  table. 


ANC 
Paramt er s 

Estimator 

ML  TOP 

MZCF 

Bandwidth 

e 

e  rms 

e 

erms 

B 

var 

A  =  1 
M  =  0.1E-2 

0.17 

0.20 

-0.2  8 

0.05 

352 

41 

A  -  1 
u  =  0.1E-8 

0.41 

0.21 

-0.15 

0.07 

A  =  8 
u  =  0.1E-2 

0.17 

0.16 

-0.25 

0.05 

385 

76 

A  =  8 
u  =  O.lE-8 

0.44 

0.37 

-0.06 

0.47 

Table  2:   Ev 


of  ANC  parameters 


The  reason  for  nsing  these  two  estimators  is  that  they  are 
the  better  of  the  class  as  noted  in  the  estimator  evaluation. 
For  both  estimators  the  different  delays  give  nearly  the  same 
resnlts  when  u  is  large.   When  it  is  small,  the  only  distinct 
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difference  is  that  the  rms  error  is  worse  for  the  delay  of  8  for 
the  HZCF  estimator.  As  for  the  bandwidth  estimate,  the  delay 
should  be  1  in  order  to  obtain  the  better  estimate. 

The  bandwidth  estimate  gives  better  results  for  the  shorter 
delay.  Hence  the  delay  is  best  set  when  it  is  one.  This  will 
also  eliminate  the  possible  loss  of  correlation  when  the  band 
rate  is  large.  The  problem  results  from  the  delay  extending  past 
the  last  bit  difference  for  which  correlation  exists. 

One  explanation  for  the  similarities  of  results  is  that  the 
gradient  constant  might  not  have  been  large  enough  to  make  a 
distinct  difference.  Perhaps  it  can  be  made  larger  by  sacri- 
ficing weight  vector  length  in  an  attempt  for  better  results.  As 
for  a  decision,  u  will  be  left  at  the  large  value  of  0.1E-2  in 
keeping  with  the  philosophy  that  this  will  provide  better 
cancel  la  t  ion. 
Varying  Frequency  Devia  t  ion 

All  of  the  previous  results  have  been  generated  for  a  fre- 
quency deviation  of  300  Hz.  This  is  the  worst  case  possible 
since  the  baud  rate  is  also  300  bps.  Therefore,  results  should 
improve  when  the  deviation  is  increased  and  the  baud  rate  is  kept 
constant.   Figures  25a  and  25b  substantiate  this  conclusion. 

In  Figure  25a,  the  results  for  the  XZCF  estimator  and  the 
HLTOI  estimator  are  presented  for  the  cases  of  Af  equal  to  300  Hz 
and  2000  Hz.  According  to  the  mean  error,  both  estimators  do  a 
better  job  for  the  larger  deviation  at  lower  SNR.  However,  the 
rms  error  has  become  slightly  worse  for  the  XZCF  estimator  and 
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FIGURE  25a:   Effect  of  varying  frequency  deviation 
of  the  FSK  signal  on  the  estimators. 
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slightly  better  for  the  MLTOI  estimator.   This  unsuspected  degra- 
dation for  the  XZCF  estimator  is  nnexp 1  a inab 1 e. 

Figure  25b  shows  how  the  bandwidth  estimate  is  affected  by 
the  frequency  deviation  of  the  FSK  signal.    The  estimate  is 
substantially  better  for  the  larger  deviation. 
Different  Band  Rates 

As  stated  in  the  introduction  the  design  should  work  for 
baud  rates  up  to  1000  bps.  Therefore,  the  results  of  the  estima- 
tors at  different  baud  rates  should  be  evaluated.  Figures  26a-c 
give  the  results  of  this  evaluation. 

The  estimates  of  the  XZCF  estimator  made  for  three  different 
baud  rates  are  plotted  in  Figure  26a.  The  three  baud  rates  are 
150,  300,  and  600  bps.  Since  the  deviation  must  be  at  least  as 
large  as  the  baud  rate,  it  was  set  to  2000  Hz  for  each  case.  The 
conclusion  to  make  is  that  the  XZCF  estimator  performs  well  for  a 
variety  of  baud  rates  down  to  a  SNR  of  8  dB.  The  150  bps  baud 
rate  has  a  significantly  smaller  rms  error  at  8  dB  which  perhaps 
indicates  the  estimator  will  work  better  for  lower  baud  rates. 
This  is  agreeable  with  the  manner  that  the  zero  crossing  analysis 
is  performed  and  with  the  way  that  errors  may  be  introduced. 

The  MLTOI  estimator  also  performs  slightly  better  at  lower 
SNR  for  the  lower  baud  rates  as  evidenced  by  Figure  26b.  There 
is  a  trend  in  the  mean  error  to  have  an  increasing  bias  as  the 
baud  rate  increases.  The  rms  error  follows  a  similar  increasing 
pattern  although  not  as  severely  as  the  mean  error. 

One  explanation  for  the  estimators  performing  better  at 
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lower  band  rates  is  that  the  sampling  frequency  is  too  low.  If 
the  sampling  frequency  is  larger,  more  samples  per  bit  would  be 
observed.  This  would  have  a  direct  effect  on  any  estimator  from 
the  weight  vector  class  of  which  the  MLTOI  estimator  is  a  member. 
It  would  also  affect  the  zero  crossing  class. 

The  resnlts  of  the  bandwidth  estimate  are  shown  in  Figure 
26c.  The  tendency  this  time  is  for  the  bias  to  increase  as  the 
baud  rate  decreases.  This  shows  that  the  bandwidth  estimate  made 
in  this  way  would  not  be  a  good  baud  rate  estimate.  However,  the 
purpose  of  the  intermediate  filter  is  to  reduce  the  noise  power 
which  the  post-filter  of  the  DLD  let  pass.  The  fluctuating  bias 
should  not  matter  since  in  all  cases  the  resulting  bandwidth  is  a 
significant  improvement  over  the  1000  Hz  bandwidth  of  the  post- 
filter,  e.g.,  the  bandwidth  estimate  for  the  150  bps  baud  rate 
case  is  225  Hz.  The  best  feature  of  the  bandwidth  estimator  is 
that  it  is  very  consistent.  As  long  as  the  SNR  is  4  dB  or  more, 
the  intermediate  filter  will  reduce  noise  power  in  a  consistent 
manner. 
Evaluation  Against  the  Traditional  Method 

The  traditional  method  for  determining  the  baud  rate  of  an 
FSK  signal  is  the  use  of  zero  crossing  methods.  Therefore  the 
proposed  method  should  be  compared  to  this  traditional  method  to 
determine  if  the  proposed  method  has  merit.  Note  that  the  only 
difference  between  the  XZCF  estimator  and  the  traditional  method 
is  the  inclusion  of  an  intermediate  filter  to  further  reduce  the 
noise  power.   The  bandwidth  of  the  filter  is  set  using  the  ANC. 
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In  order  to  make  the  judgement,  the  results  of  these  two  methods 
along  with  the  results  nsing  the  MLTOI  estimator  are  plotted  in 
Figures  27a-d  for  various  baud  rates  and  signal  deviations. 

In  Figure  27a  the  baud  rate  is  300  bps  and  the  deviation  is 
300  Hz.  The  accuracy  of  the  traditional  method  drops  somewhere 
above  12  dB  while  both  of  the  proposed  methods  are  good  down  to 
nearly  8  dB. 

The  baud  rate  in  Figure  27b  is  600  bps  and  the  deviation  is 
2000  Hz.  This  time  the  MLTOI  estimator  is  the  poorest  estimator. 
The  ZC  and  the  XZCF  estimator  give  nearly  identical  results. 
This  could  be  for  two  reasons.  One,  the  addition  of  the  interme- 
diate filter  at  about  600  Hz  bandwidth  did  not  make  a  significant 
reduction  in  noise  power.  Two,  the  larger  baud  rate  means  more 
samples  per  bit  time,  and  thus,  a  steadier  zero  crossing  analysis 
can  be  made. 

When  the  baud  rate  is  decreased  to  150  bps  with  the  same 
deviation,  as  in  Figure  27c,  the  XZCF  estimator  does  indeed 
perform  better  than  the  traditional  ZC  estimator.  However,  the 
MLTOI  estimator  performs  the  best  of  the  three.  Its  results  are 
good  to  4  dB. 

Figure  27d  prompts  the  same  conclusions  for  a  baud  rate  of 
300  bps  and  a  signal  deviation  of  2000  Hz.  This  time,  however, 
the  MLTOI  estimator  contains  a  bias  which  is  une  xp  1  a  inab  1  e. 
Still,  its  estimates  are  relative  accurate  to  4  dB.  Again,  the 
traditional  ZC  method  is  useless  below  12  dB. 
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FIGURE   27a:.     Evaluation   of    proposed'  estimators    against 
the   traditional   zero-crossing   estimator 
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FIGURE  27b:   Evaluation  of  proposed  estimators  against 
the  traditional  zero-crossing  estimator. 
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FIGURE  27c:       Evaluation   of   proposed   estimators    against 
the   traditional   zero-crossing   estimator. 
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FIGURE  27d:   Evaluation  of  proposed  estimators  against 
the  traditional  zero-crossing  estimator. 
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Cone  Ins  ions 

Two  classes  of  estimators  proposed  provide  better  results 
than  the  traditional  zero  crossing  method.  Both  methods  rely  on 
the  ANC,  one  directly  and  the  other  indirectly. 

The  weight  vector  class  of  estimators  uses  the  mean  and 
variance  of  the  weight  vectors  of  the  iterative  convergence 
procedure  to  obtain  the  band  rate.  Of  the  estimators  of  this 
class,  the  MLTOI  estimator  was  found  to  be  one  of  the  better 
estimators.  Figure  28  shows  how  the  MLTOI  estimate  is  acquired. 
The  estimate  defines  the  bit  time  in  number  of  sampling  intervals 
as  being  the  delay  of  the  ANC  plus  the  number  of  the  weight  which 
represents  the  first  relative  minimum  of  the  weight  vector  which 
is  less  than  zero.  The  bit  interval  is  then  obtained  from  the 
weight  number  plus  delay  by  multiplying  by  the  sampling  interval. 
Then  the  baud  rate  is  found  by  inverting  the  bit  interval.  The 
estimator  gives  good  results  over  a  large  part  of  the  desired 
range  of  baud  rates  and  frequency  deviations.  Its  estimates  are 
reliable  down  to  8  dB  in  most  instances.  This  estimate  performs 
better  than  the  traditional  method  for  nearly  all  baud  rates  and 
appreciably  better  for  low  baud  rates. 

The  other  class  of  estimators  relies  indirectly  on  the  ANC. 
This  time  the  ANC  is  used  to  obtain  a  bandwidth  estimate  of  the 
message.  Then  intermediate  filtering  is  performed  using  the 
bandwidth  estimate  to  further  reduce  the  noise  power.  The  baud 
rate  estimate  is  then  obtained  using  zero  crossing  analysis  just 
as  the  traditional  method  does.   The  analysis  is  performed  on 
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either  the  filtered  ANC  input  or  on  the  various  outputs  from 
either  one  of  the  two  ANC's.  The  best  results  are  obtained  for 
the  XZCF  estimator  which  is  when  the  zero  crossing  analysis  is 
performed  on  the  intermediate  filter  output.  This  is  illustrated 
in  Figure  28.  Obviously,  this  estimate  will  perform  at  least  as 
well  as  the  traditional  method.  When  the  band  rate  drops  to  the 
lower  end  of  the  range,  this  estimator  performs  much  better  than 
the  traditional  method.  This  is  a  direct  result  of  the  addi- 
tional noise  power  cancellation  attributed  to  good  bandwidth 
estimation  by  the  ANC.  However,  for  the  upper  range  of  baud 
rates  and  the  lower  range  of  frequency  deviations,  this  estimator 
will  degrade  just  as  fast  as  the  traditional  method  when  the  SNR 
is  decreased. 

The  bandwidth  estimation  provided  by  the  ANC  should  be 
considered  as  a  significant  achievement.  The  estimate  rarely 
underestimates  the  bandwidth  and  the  accuracy  remains  stable  down 
to  4  dB  SNR.  Obviously  this  addition  to  the  system  is  only 
valuable  for  the  surveillance  situation  since  the  bandwidth  would 
be  known  in  other  aplications. 

One  factor  which  definitely  influences  the  estimators  is  the 
sampling  frequency.  When  set  relatively  low  in  regards  to  the 
larger  baud  rates,  the  number  of  samples  per  bit  is  small.  This 
can  create  errors  and  biases  in  either  class  of  estimators.  The 
sampling  frequency  should  be  set  perhaps  at  least  20  times  the 
largest  expected  baud  rate.  Unfortunately,  this  means  the  weight 
vector  must  contain  a  large  number  of  weights. 
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Another  item  which  deserves  attention  is  the  weight  vector. 
The  weight  vector  is  supposed  to  converge  to  s  stable  weight 
vector  which  will  produce  an  error  output  of  the  ANC  of  the  least 
power.  In  this  sense  the  ADP  filter  weights  will  have  converged 
to  the  optimal  Wiener  filter  [1].  An  analysis  of  the  weight 
vector  would  be  beneficial,  especially  in  the  case  of  a  random 
bit  input  sequence  as  modeled  in  the  simulation.  Once  the  analy- 
sis is  performed  it  will  be  possible  to  decide  if  the  pseudo- 
correlation  sub-class  of  estimators  contain  merit,  and  if  they 
do,  how  well  they  will  perform  in  comparison  to  the  other  pro- 
posed estimators. 

Another  aspect  the  analysis  be  concerned  with  is  in  the 
method  of  noise  (or  signal)  cancelling.  The  method  used  here 
relies  on  the  ANC  and  can  be  shown  to  be  in  the  class  of  auto- 
regressive  predictors  [5].  This  class  of  estimators  is  known  for 
a  spectrum  which  is  characterized  by  peaks.  The  peaks  result 
from  a  rational  model  for  this  estimator  involving  an  all-pole 
transfer  function.  However,  since  the  spectrum  of  the  bit 
sequence  is  characterized  more  so  by  nulls  in  the  frequency 
domain  than  by  peaks,  the  bit  sequence  should  be  considered  as  a 
process  other  than  an  auto-regressive  process.  Therefore,  it 
makes  sense  to  use  an  estimator  which  is  compatible  with  nulls, 
or  in  general,  nulls  and  peaks.  The  a u t o- r e g r e s s i v e ,  moving 
average  estimator  meets  this  requirement  since  its  rational  model 
is  characterized  with  a  transfer  function  containing  both  zeros 
and  poles. 


Appendix  A:   Computer  Program  Listin 


gs 


Program  page 

SBAUD  9$ 

SADF  13  2 

SSADF  145 

BITLENGTH  151 

ZCROSS  155 

DLD  160 
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VAX-11  FORTRAN  SOURCE  FILENAME: 


DEPARTMENT  OF  ELECTRICAL  ENGINEERING 


SBAUD.FOR 


KANSAS  STATE  UNIVERSITY 


REVISION 


00.0 


DATE 


MAT  6,  1983 


PROGRAMMER (S) 


MARC  D.  BRACE 


£»•••»•••*»•••••»*•»»•••*»»«»••»»•*•»•«»•»•»••»»••»•••*»*«»•«••»••»»««»» 

C 

C       SBAUD 

C 

C 

c 
c 
c 

C 

c 
c 
c 

CO********************************************************************* 

c 

C       PURPOSE 

C 

C 

c 

C 
C 
C 

c 
c 
c 
c 

C  DLD 

C  SADF 

C  SSADF 

C  ZCROSS 

C 

C  SGOPEN 

C  SGTRAN 

C  SPLOT 

C  BUTER 

C  DFT 

C  GAUSS 

C 

C 

C 

C 

C 

C 

IMPORTANT:   Weight  vector  length  must  be  a  power  of  two  if 
the  intermediate  filter,  the  psendo-correlation 
estimates,  or  the  frequency  response  of  the  weight 
vector  data  are  to  be  implemented.   The  DFT  routine 
used  in  these  applications  requires  a  vector  of 
length  of  a  power  of  two. 


This  program  is  a  simulation  for  a  system  that  estimates 
the  band  rate  of  an  FSK  signal.   The  system  consists  of 
a  delay-line  discriminator  (DLD)  to  demodulate  the  signal 
and  of  an  adaptive  noise  canceller  (ANC)  to  make  the 
estimate.  The  results  are  compared  to  the  results  of  the 
traditional  method  of  zero  crossing  estimation. 

ROUTINE(S)  ACCESSED  OR  CALLED  BY  THIS  ROUTINE 


These  routines  are  specialized  routines  written 
only  for  the  purpose  of  the  simulation.   They 
are  compatible  with  SBAUD  and  each  other  and  are 
meant  to  be  run  in  a  single  package. 

These  routines  are  general  purpose  routines  that 
perform  such  functions  as  writing  to  disk, 
plotting,  filtering,  and  making  FFT's.   Any  other 
routines  performing  identical  functions  may  be 
substituted  for  them. 


ROUTINE  (S)  ACCESSED  OR  CALLED  BY  EMBEDDED  ROUTINE  (S) 

BITLENGTH  This  is  another  specialized  routine. 
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C       Variables  used  in  main  routine: 

C 

C       ADF     TRUE  if  ADF  is  to  be  included  in  simulation 

C 

C       AI      Amplitude  of  FSE  signal  envelope 

C 

C       ASG     Argument  of  FSE  signal  envelope 

C 

C      AUTO    TRUE  if  autocorrelation  estimate  is  enabled 

C 

C  BITS  Number  of  bits    in  message 

C 

C       BN      Blank  character  for  plotting  purposes 

C 

C       BOTH    TRUE  if  two  files  are  to  be  written  to  disk 

C 

C       BUTER   Complex  butterworth  filter  rontine 

C 

C       BW      Max  bandwidth  of  signal  according  to  sampling  freqency 

C 

C       C       User  response  as  Y  or  N 

C 

C       CHOICE  Pointer  of  directed  GOTO  statements 

C 

C       CI      First  estimate  to  be  written 

C 

C       C2      Number  of  estimates  to  be  written 

C 

C       DELAY   Number  of  samples  to  delay  in  ANC 

C 

C       DW      FSE  frequency  shift  in  radians 

C 

C       FC      Center  frequency  of  FSE  signal 

C 

C       FFT     TRUE  if  DFT  is  to  be  made  of  weight  vector 

C 

C       FLAG    Indicates  sin-  or  cosine-type  discriminator 

C 

C       FO      Fundamental  frequency  of  entire  message 

C 

C       FREQP   Max  frequency  of  frequency  response  to  be  plotted 

C 

C       FS      Sampl ing  frequency 

C 

C       FSE     TRUE  if  FSE  signal  information  is  to  be  obtained 

C 

C       FW      FSE  frequency  shift  in  Hertz 

C 

C       HT      Horizontal  title 

C 

C       HU      Horiaontal  units 

C 
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C       I       Loop  counter  variable 

C 

C       II      Loop  counter  variable 

C 

C       ITER    Simulation  loop  counter  variable 

C 

C       J       Loop  counter  variable 

C 

C       LIM     TRUE  if  limiter  is  to  be  included  in  sumnlation 

C 

C       LTI     Left  title 

C 

C       LUN     Left  units 

C 

C       MAXF    Max  of  DFT  of  weight  mean  or  variance 
C 

C       MU      Gradient  convergence  rate  constant 

C 

C       N       Number  of  samples  in  signal 

C 

C       NN      Increment  factor  for  plotting  purposes 

C 

C       ND      Delay  for  DLD 

C 

C       NOISE   TRUE  if  noise  is  to  be  added  in  respective  iteration 
C 

C       NSBWM   Normalized  bandwidth  estimate  mean 
C 

C       NSBWS   Normalized  bandwidth  standard  deviation 
C 

C       OFILT   TRUE  if  intermediate  filtering  is  to  be  done 
C 

C       P       Factor  for  either  shifteing  up  or  down  in  FSK  signal 

C       PI      The  constant  pi 

C 

C       PNUM    Number  of  points  to  plot 

C 

C       POSTA   TRUE  if  post-filtering  is  to  be  done  for  ANC 
C 

C       POSTBW  Post-filter  bandwidth  for  ANC 

C 

C       POSTBZ  Post-filter  bandwith  for  ZC  analyisis  (if  POSTA  false) 

C       PRE     TRUE  if  pre-filtering  is  to  be  done 
C 

C       PREBW   Pre-filter  bandwidth 

C 

C       R       Baud  rate 

C 

C       RTI     Right  title 

C 
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C       HUN     Right  units 

C 

C       SANCZC  TRUE  if  ZC  analysis  is  to  be  done  in  ANC 

C 

C       SBW     Bandwidth  estimate 

C 

C       SBWM    Bandwidth  estimate  mean 

C 

C       SBWS    Bandwidth  estimate  mean  square 

C 

C       SBW    Bandwidth  estimate  variance 

C 

C       SCORR   TRUE  if  pseudo-correlation  analysis  is  to  be  done 

C 

C       SMM     TRUE  if  intermediate  estimates  are  to  be  made 

C 

C       SNR     TRUE  if  noise  is  to  be  added  to  signal 

C 

C       SNRDB   Signal-to-noise  ration  in  dB 

C 

C       SFB     Samples  per  bit 

C 

C       STITLE   Sub-title  for  plot 

C 

C       TB      Bit  interval 

C 

C       ID      Time  delay  for  discriminator 

C 

C       TF      Time  interval  for  FSK  signal 

C 

C       THETA   Phase  offset  of  FSK  signal 

C 

C       TITLE   Major  title  for  plot 

C 

C       TRIALS  Number  of  iterations  for  simulation  loop 

C 

C       TS      Sampling  interval 

C 

C       VAR     Variance  for  noise 

C 

C       WBIAS   TRUE  if  bias  is  to  be  ioncluded  in  weight  vector 

C 

C       WC      Center  frequency  in  radians  per  second 

C 

C       WINDOW  TRUE  if  window  is  to  be  used  in  DFT  situations 

C 

C       WNUM    Number  of  weights  in  weight  vector 

C 

C       WRIT    TRUE  if  data  is  to  be  written  to  disk 

C 

C       WTNUM   Length  of  zero  padded  weight  vector  in  DFT  situations 
C 
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C       XSCALE   Scaling  of  abscissa  values  for  plotting  purposes 

C 

C       ZC      TRUE  if  traditional  ZC  anlysis  is  to  be  done 

C 

(;•«•«••••»*«••«»••«•»•»»»»••••«»**»•••«»»»••»••««»••»»»»•»««»»««»,,,»»». 

C 

C       Arrays  used  in  the  main  routine: 

C 

C       AUTOC   Used  in  obtaining  autocorrelation 

C 

C       CNOISE  Noise  to  be  added  to  signal 

C 

C       DLDNN   Output  of  DLD  for  noiseless  input 

C 

C       EROR    Error  output  of  secondary  ANC 

C 

C       ERORI   Error  output  of  intermediate  ANC 

C 

C       FILT    DLD  output  intended  of  tradistional  zc  analysis 

C 

C       FSIG    FSK  signal  magnitude 

C 

C       LEFT    Left  array  for  plotting 

C 

C       NDLD    DLD  output,  general  case 

C 

C       NAUTO   Noise  autocorrelation 

C 

C       PFSK    FSK  signal  phase 

C 

C       RIGHT   Right  array  for  plotting 

C 

C       S       FSK  envelope  and  general  purpose  complex  array 

C 

C       SAUTO   Message  autocorrlation 

C 

C       SEQ     Bit  sequence  for  of  message 

C 

C       SFEST   Secondary  baud  rate  estimates 

C 

C       SFESTH  Mean  of  secondary  baud  rate  estimates 

C 

C       SFESTS  Mean  square  of  secondary  baud  rate  estimates 

C 

C       SFESTV  Variance  of  secondary  baud  rate  estimates 

C 

C       SIEST   Intermediate  baud  rate  estimates 

C 

C  SIESTM     Mean  of    intermediate  baud  rate   estimates 

C 

C       SIESTS  Mean  square  of  intermediate  baud  rate  estimates 
C 
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C       SIESTV  Variance  of  intermediate  band  rate  estimates 

C 

C       WMEAN   Weight  vector  mean 

C 

C       WMEANA  Average  of  all  weight  vector  means 

C 

C       WVAR    Weight  vector  variance 

C 

C      WVARA   Average  of  all  weight  vecto  variances 

C 

C       Y       ADF  output  of  secondary  ANC 

C 

C       YI      ADF  output  of  intermediate  ANC 

C 

C       ZEST    Traditional  ZC  band  rate  estimate 

C 

C       ZESTH   Mean  of  traditional  ZC  band  rate  estimate 

C 

C       ZESTS   Mean  square  of  traditional  ZC  baud  rate  estimate 

C 

C       ZESTV   Variance  of  traditional  ZC  baud  rate  estimate 

C 

get********************************************************************* 

c 

IMPLICIT  NONE 

REAL  THETA,BW,AI,FSIG(2048) ,FILT(2048> ,MDLD(2048) ,MD,ARG 

REAL  TD,FC,TF,PREBW,POSTBW,PFSK(2048),FW,DW,R,FO,DLDNN(2048) 

REAL  VAR,SNRDB,SAUTO(2048),NAUTO(2048),SBWV 

REAL  WMEAN(0:256,200),WVAR(0:256,200),WMEANA(0:256),WVARA(0:256) 

REAL  XSCALE,LEFT(2048),RIGHT(2048),SBW,SBWS,SBWM,MAXF 

REAL  ZEST(10,200),SIEST(10,200),  SFEST(10,200)  .POSTBZ.FREQP 

REAL  ZESTJC(IO)  ,SIESTM(10)  ,SFESTM(10)  ,NSBWM,NSBWS 

REAL  ZESTS(10),SIESTS(10),SFESTS(10) 

REAL  ZESTV(10),SIESTV(10),SFESTV(10) 

REAL  EROR(2048) ,Y(2048) ,ERORI(2048) ,YI(2048) 

REAL*8  PI,WC,FS,TS,TB 

INTEGER  SFB,P,I,J,SEQ(64),BrTS,N,NN,FLAG,WNUM,DELAY 

INTEGER  ITER,ND,CHOICE,TRIALS,PNnM,Cl,C2,II,WTNUM 

COMPLEX  S(2048),BUTER,CNOISE(2048),AUT0C(2048) 

CHARACTER»80  BN, TITLE, LTI,RTI,STITLE,HT,HU,LDN, RUN 

CHARACTER  C»l 

LOGICAL  NOISE,  WBIAS,  AUTO, ZC,  ADF,  FSK.POSTA,  WINDOW,  SANCZC 

LOGICAL  PRE,SNR.OFILT,LIM,SMM,BOTH,WRIT,FFT,SCORR 

NN  =  1 

BN  =    *    ' 

PI  -  3.141592653589793D0 

N  «  1024 

FS  =  9600 

FO  =  FS/N 

TS  =  l./FS 

BW  =  FS/2. 
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FC  =  500*FS 
WC  -  2»PI»FC 
BITS  =32 
SFB  =  N/BITS 
TB  =  TS'SPB 
R  =  1.0/TB 
AI  =  1.0 
THETA  =  0.0 
NOISE  =  .FALSE. 
C 

CM********************************************************************* 

c 

C  PARAMETER  INPUT 

C 

C  AUTOCORRELATION  ESTIMATION  ENABLE 

C 

CALL  PAUTO(AUTO) 
C 

C  FSE  SIGNAL  INFORMATION  ENABLE 
C 

CALL  PFSKINFO(FSK) 
C 

C   SIMULATION  ADAPTIVE  FILTERING  SUBROUTINE  ENABLE 
C 

CALL  PADF(ADF) 
C 

C  ZERO  CROSSING  ALGORITHM  ENABLE 
C 

CALL  PZCROSS(ZC) 
C 

C  INPUT  SIGNAL  PARAMETERS 
C 
C 

TYPE  •,'  ' 

TYPE  *,  'INPUT  SIGNAL  PARAMETERS:' 
C 

C   SEQUENCE  FORMATION 
C 

CALL  PSEQUENCE(SEQ.BITS) 
C 

C  FSE  PARAMETERS 
C 

CALL  PFSKSHIFT(FW) 

DW  =  2.0»PI»FW 
C 

C  NOISE  SELECTION 
C 

CALL  PNOISE(SNRDB,SNR) 
C 

C  DELAY  LINE  DESCRIMINATOR  PARAMETERS 
C 

TYPE  •,  '  ' 
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TYPE  *,    'DLD  RECEIVER  PARAMETERS:' 
C 

C     FILTER  SELECTION 
C 

CALL  PPREFILTER(PREBW,PRE) 

CALL  PPOSTFILTER(POSTBW,POSTA,ZC,POSTBZ) 
C 

C  LIMITER  ENABLE 
C 

CALL  PLIM(LIM) 


C 

c 

C 
C 

DLD  SHIFT  PARAMETER 

CALL  PDELAT(ND) 

c 

C 

c 

ND  =  INT(FC/(6»PREBW)) 

ADAPTIVE  FILTER  PARAMETERS 

IF(ADF)  THEN 

TYPE  •,  '  ' 

c 

C 

C 

C 
C 

c 

C 
C 
C 

c 

C 

c 

TTPE  •,  'ADAPTIVE  FILTERING  PARAMETERS:* 

ADF  DELAY  SELECTION 

CALL  PADFDELAY (DELAY) 

NUMBER  OF  ADF  WEIGHTS  SELECTION 

CALL  PWEIGHTS'WNUM) 

GRADIENT  CONSTANT  SELECTION 

CALL  PCONSTANT(MD) 

WEIGHTING  BIAS  SELECTION 

CALL  PBIAS(WBIAS) 

TYPE  *,  '  ' 

TYPE   •,    'ADF  ESTIMATOR  SELECTION:  ' 
C 

C     INTERMEDIATE  FILTERING  ENABLE 
C 

CALL  PFTLTl(OFILT) 
C 

C  MAX/MIN  ESTIMATE  ENABLE 
C 

CALL  PHM(SMM) 
C 

C  PSEUDO-CORRELATION  ESTIMATE  ENABLE 
C 

CALL  PCORR(  SCORE) 
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C 

C  ANC-ZC  ESTIAMTE  ENABLE 

C 

CALL  PANCZC(SANCZC) 
C 

ENDIF 
C 

C  NUMBER  OF  TRIALS  SELECTION 
C 

CALL  PTRIALS(TRIALS) 
C 
(■**«•«»••***•**«***»**«*••*»»««»*»*«****•»•»»*»•**•*»»*«****»»»»««»««*»• 

C 

C  RUN  INITIALIZATION 

C 

1000    CONTINUE 

SBWM  =  0 
SBWS  =  0 
NSBWM  =  0 
NSBWS  -  0 
DO  140  1=1,10 

SIESTM(I)  -  0 

SIESTS(I)  =  0 

ZESTM(I)  =  0 

ZESTS(I)  =  0 

SFESTM(I)  =  0 

SFESTS(I)  =  0 
140     CONTINUE 

NOISE  =  .FALSE. 

DO  150  J=0,WNUM 
TOEANA(J)  =  0. 
WVARA(J)  =  0. 
150     CONTINUE 

IF  (AUTO)  THEN 
DO  160  1=1, N 
160         NAUTO(I)  =  0. 

ENDIF 
C 

£•***•••*•*•*****•***•••******•••**«**•*****»********••***•••*******•**» 
c 

C  SIMULATION  DO  LOOP 
C 

DO  1500  ITER=1, TRIALS 

TYPE  •,  '  ' 

TYPE  •,  'IN  ITERATION  NUMBER  '.ITER 

1100    CONTINUE 
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C 

c 

C  FSK  GENERATION 
C 

DO  180  1=1, BITS 

IF  (SEQ(I).EQ.l)  THEN 

P  =  1 
ELSE 

P  =  -1 
ENDIF 
C 

C  MODULATING  FUNCTION 
C 

DO  170  J=1.SPB 

S((I-1)*SPB+J)  =  AI»CDEXP(DCMPLX(0.0D0,P»((I-1)*SPB+J-1)«TS*DW)) 
170       CONTINUE 

180     CONTINUE 

C 

C 

C  NOISE  ADDITION 

C 

IF  (NOISE)  THEN 

CALL  SNOISE(CNOISE,N,AI,TS,PREBW,SNEDB,SNR) 

DO  190  1=1, N 
190         S(I)  =  S(I)  +  CNOISE(I) 

ENDIF 
C 

C  OBTAIN  FREQUENCY  COMPONENTS  OF  FSK 
C 

IF(FSK)  THEN 

CALL  DFT(S,N,0) 

DO  200  1=1, N 

FSIG(I)  =  CABS(S(D) 
200       CONTINUE 

CALL  DFT(S,N,1) 
C 

C  OBTAIN  PHASE  OF  FSK  ENVELOPE 
C 

DO  210  1=1, N 
210         PFSK(I)  =  ATAN2(AIMAG(S(I)),REAL(S(I)))/PI 

ENDIF 
C 
£**•**•*•*•••**«»•*********••* 

C 

C  DELAY-LINE  DISCRIMINATOR 

C 
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C 

C   PRE-FILTERING 

C 

IF  (PRE)  THEN 
CALL  DFT(S,N,0) 
DO  300  1=1, N/2 

S(I)  =  S(I)*BUTER(4,PREBW,FLOAT(I-l)*FO) 
SCN+l-I)  =  S(N+l-I)*BUTER(4,PREBW,FLOAT(-I)*FO) 
300       CONTINUE 

S(N/2+l)  -  CMPLX(0.,0.) 
CALL  DFT(S,N,1) 
ENDIF 
C 

C  LIMITER 
C 

IF  (LIM)  THEN 
DO  305  1=1. N 

IFC.NOT.  (REAL(S(D).EQ.O  .AND.  AIHAG(Sd) )  .EQ.O)  )THEN 

ARG  =  ATAN2(AIHAG(S(I)),REAL(S(I))) 
ELSE 

ARG  =  0 
ENDIF 

S(I)  =  AI*CMPLX(COS(ARG>,SIN(ARG)) 
305       CONTINUE 

ENDIF 
C 

C  DLD 
C 

TD  =  ND/(2.0*FC) 
FLAG  =  1 
TF  =  N*TS 

CALL  DLD(TD,FC,FLAG,S,N,TF) 
C 

C  POST-FILTERING  FOR  ADAPTIVE  FILTERING 
C 

IF  (POSTA)  THEN 
CALL  DFT(S,N,0) 
DO  310  1=1, N/2 

S(I)  =  S(I)*BUTER(4,POSTBW,FLOAT(I-l)*FO) 
S(N+1-I)  =  S(N+1-I)*BUTER(4,POSTBW,FLOAT(-I)*FO) 
310       CONTINUE 

S(N/2+l)  =  CMPLX(0.,0.) 
CALL  DFT(S,N.l) 
ENDIF 
C 
£••***•»***•••*•»••»*•**•***•« 

C 

C  INITIAL  NOISELESS  CASE 

C 

IF  (.NOT. NOISE. AND. SNR)  THEN 
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NOISE  =  .TRUE. 
DO  400  1=1, N 
400         DLDNN(I)  =  REAL(SU)) 
C 

C  MESSAGE  AUTOCORREALATION 
C 

IF  (AUTO)  THEN 
DO  460  1=1,  N 
460  AUTOC(I)  =  CMPLX(DLDNN(I),0.) 

CALL  DFT(AUTOC,N,0) 

DO  470  1=1, N 
470  AUTOC(I)  =  CMPLX(CABS(AUTOC(I))**2,0.) 

CALL  DFT(AUT0C,N,1) 

DO  480  1=1. N 
480  SAUTO(I)  =  CABS(AUTOC(I)) 

ENDIF 

GOTO  1100 
ENDIF 
C 

C 

C  NOISE  AUTOCORRELATION  ESTIMATE 

C 

IF  (AUTO)  THEN 
DO  410  1=1,  N 
410         AUTOC(I)  =  CMPLX(REAL(S(I))-DLDNN(I),0.) 
CALL  DFT(AUTOC,N,0) 

DO  420  1=1, N 
420         AUTOC(I)  =  CMFLX(CABS(AUTOC(I))*'2,0.) 
CALL  DFT(AUTOC,N,l) 

DO  430  1=1, N 
430         NAUTO(I)  =  NAUTO(I)  +  CABS(AUTOC(I) ) /TRIALS 
ENDIF 

DO  440  1=1, N 

440  MDLD(I)  =  REAL(SU)) 
C 

C  PERFORM  ZERO  CROSSING  ESTIMATE,  IF  NO  ADF  POST  FILTER,  INCLUDE  FILTERING 
C 

IF(ZC)  THEN 

IF(.NOT.POSTA)  THEN 
DO  441  1=1, N 

441  CNOISE(I)  =  S(I) 
CALL  DFT(CNOISE,N,0) 
DO  442  1=1, N/2 

CNOISE(I)  =  CNOISE(I)»BUTER(4,POSTBZ,FLOAT(I-l)»FO) 
CNOISE(N+l-I)  «  CNOISE(N+l-I)*BUTER(4,POSTBZ,FLOAT(-I)*FO) 
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442  CONTINUE 

CNOISE(N/2+l)  =  CMPLX(0.,0.) 
CALL  DFT(CNOISE,N.l) 

DO  443  1=1, N 

443  F1LT(I)  =  REAL(CNOISEd)) 

CALL  ZCROSS(FILT,N,2TCST(l,rTER),l) 
ELSE 

CALL  ZCROSS(MDLD,N,2EST(l,ITER),l) 
ENDIF 
ENDIF 
C 

C 

C  ADAPTIVE  NOISE  FILTERING  SUBROUTINE 

C 

IF  (ADF)  THEN 

CALL  SSADF(MDLD,N,MU,WNUM,DELAT,WBIAS,WMEAN(0, ITER)  ,WVAR(0, ITER), 
X  OFTLT,SFEST(l,  ITER)  ,TS,SIEST(1,  ITER),  SBW, SMM.ERORI.YI, 

X  EROR,T,SCORR,SANCZC) 

C 

C     BANDWIDTH  ESTIMATE 
C 

IF(OFILT)    THEN 

SBWM  =  SBWM  +  SBW/TRIALS 
SBWS  =  SBWS  +  SBW*»2/TRIALS 
NSBWM  =  NSBWM  +  ((SBW  -  R)/R) /TRIALS 
NSBWS  =  NSBWS  +  ((SBW  -  R)  /R)  **2/TRIALS 
ENDIF 
C 

C  BAUD  RATE  ESTIMATE  STATISTICAL  UPDATE 
C 

IF(SMM)  THEN 
DO  444  1=1,10 

SIESTM(I)  =  SIESTM(I)  +  ( (SIEST(I,ITER)-SPB) /SPB) /TRIALS 
SIESTS(I)  =  SIESTS(I)  +  ( (SIEST(I.ITER)-SPB) /SPB)**2/TRIALS 

444  CONTINUE 
ENDIF 

IF  (OFILT)  THEN 
DO  445  1=1,10 

SFESTM(I)  =  SFESTM(I)  +  ((SFEST(I,ITER)-SPB) /SPB) /TRIALS 
SFESTS(I)  =  SFESTS(I)  +  ( (SFEST(I,ITER)-SPB)/SPB) **2/TRIALS 

445  CONTINUE 
ENDIF 

DO  450  I=0,WNUM 

WMEANA(I)  =  WMEANA(I)  +  WMEAN(I,  ITER) /TRIALS 
WVARA(I)  =  WVARA(I)  +  WVAR( I, ITER)  /TRIALS 
450       CONTINUE 
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ENDIF 

IF(ZC)  THEN 
DO  455  1=1,1 

ZESTM(I)  =  ZESTM(I)  +  ( (ZEST(I,ITER)-SPB)  /SPB)  /TRIALS 
ZESTS(I)  =  ZESTS(I)  +  ( (ZEST(I,ITER)-SPB) /SPB) ••2/TRIALS 
455       CONTINUE 
ENDIF 

1500    CONTINUE 

C 

C  END  OF  SIMULATION  LOOP 

C 

(;»••••»••••••»•••»•»•»»•«»*••••»»«»••»«**•«»•••»»*•»»»»»»•«««»«.««,«,,,» 

C 

C  OBTAIN  VARIANCE  OF  STATISTICS 

C 

IF  (ADF)  THEN 
IF(OFILT)  THEN 

SBWV  =  SBWS  -  SBWM«»2 
IFfSBWV.GE.O.)  SBWV  =  SBWV*«0.5 
NSBWS  =  NSBWS**0.5 
ENDIF 

IF  (SMM)  THEN 
DO  490  1=1,10 

SIESTV(I)  =  SIESTS(I)  -  SIESTM(I)»*2 
IF(SIESTVtt).GE.O.O)  SIESTV(I)  =  SIESTV(I)«0.5 
490         CONTINUE 
ENDIF 

IF(OFILT)  THEN 
DO  495  1=1.10 

SFESTV(I)  =  SFESTS(I)  -  SFESTM(I)»*2 
IF(SFESTVCI).GE.O.O)  SFESTV(I)  =  SFESTV(I)*»0.5 
495         CONTINUE 
ENDIF 
ENDIF 

IF  (ZC)  THEN 
DO  500  1=1,1 

ZESTV(I)  =  ZESTS(I)  -  ZESTM(I)»»2 
IF(ZESTVU).GE.O.O)  ZESTV(I)  =  ZESTV(I)**0.5 
500       CONTINUE 

ENDIF 
C 

Ca<«M<»MttMMMM««M«Mt«tM«»«M«t»««»»tM«Mat««M««taMMMt,M 

c 

C     PLOTTING  AND  MODIFICATION 
C 

CALL  PTITLE(TITLE,PRE,POSTA,FW,PREB¥,POSTBW,SNRDB,SNR,R,LIM) 
IF    (TOIAS)    THEN 

WRITE(STITLE,14)    'ADF-DELAT=    '.DELAY,',   N=',WNUM,',   MU=',MU ,',    BIAS- 
ELSE 
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WRITE(STITLE,14) 
ENDIF 


'ADF-DELAY=  '.DELAY,',  N=',WNUM,',  Mn=',MU 


2000    TYPE  *,  *  ' 

TYPE  22 
22      FORMATC  Do  yon  want  to:',/ 


i 
z 
i 
i 
i 


1.  Plot',/ 

2.  Change  parameters',/ 

3.  Enter  user  controlled  ADF  routine', 

4.  Write  results*,/ 

5.  Quit',/ 

6.  Write  output  to  disk') 


TYPE  11,  '/Enter  selection: 
READ  (•,•)  CHOICE 


WRIT  =  .FALSE. 
IF(CHOICE.EQ.6)  WRIT  = 


.TRUE. 


GOTO  (5000,6000,7000,8000,9000,5000)  CHOICE 
C 

C 

C  PLOTTING  AND  WRITING  TO  DISK 

C 

5000 

23  'ORMATC   Do  you  want    to  plot(write) : ' ,/ 

DLD  output',/ 

Frequency  components  of  FSK  signal',/ 

Phase  of  FSK  envelope',/ 

Autocorrelation  estimates',/ 

Weight  mean  and  variance  data',/ 

Intermediate  ANC  error',/ 

Intermediate  ADF  output',/ 

Secondary  ANC  input',/ 

Secondary  ANC  error',/ 

Secondary  ADF  output',/ 

Or  do  you  want  to  return  to  main  menu') 
TYPE  11,  '/Enter  selection:    * 
READ(*,*)  CHOICE 

GOTO(5100, 5200, 5300, 5400, 5500, 5600, 5700, 5750, 5800, 5900, 2000)  CHOICE 
C 

c 

C  DLD  OUTPUT 

C 

5100    DO  5110  1=1, N 

5110      LEFT(I)  =  REAL(SU)) 


TYPE   •, 

»    , 

TYPE  23 

FORMAT ( ' 

Do  yo 

z                  , 

'        1. 

z 

*        2. 

■                    i 

3. 

x                   > 

4. 

z                   , 

5. 

z                   , 

6. 

z                   , 

'        7. 

*                   , 

8. 

z                   , 

9. 

z                   , 

'      10. 

z                    , 

*      11. 

WTNUM  =  N 
BOTH  =  .FALSE. 


Ill 


IF(WRIT)    GOTO  5998 

LTI  =    'REAL  PART  OF  DLD  OUTPUT' 

LUN  =    '    ' 

RTI  =    'NONE' 

HT  =  'SAMPLE  NUMBER' 

HU  =  '  ' 

STITLE  -  '  ' 

XSCALE  -  1 

PNUM  =  N 

GOTO  5999 
C 

c 

C  FREQUENCY  COMPONENTS  OF  FSK  SIGNAL 
C 

5200    IF(.NOT.FSK)  THEN 
TYPE  •,  '  ' 
TYPE  *,  'NOT  ENABLED' 
GOTO  2000 
ENDIF 

DO  5210  1=1, N 
5210      LEFT(I)  =  FSIG(I) 

TTNUM  ■  N 
BOTH  »  .FALSE. 
IF(WRIT)  GOTO  5998 

LTI  =  'FREQUENCY  COMPONENTS  BEFORE  DLD' 

LUN  =  '  ' 

RTI  =  'NONE* 

HT  -  'FREQUENCY' 

HU  =  'HERTZ' 

STITLE  -  *  ' 

XSCALE  =  FO 

PNUM  =  N/2 

GOTO  5999 
C 

£:••••••*»«»»*••*»•«»»»»*»»••»• 

c 

C  PHASE  OF  FSK  ENVELOPE 
C 

5300    IF(.NOT.FSK)  THEN 
TYPE  •,  '  ' 
TYPE  *,  'NOT  ENABLED' 
GOTO  2000 
ENDIF 

DO  5310  1=1, N 
5310      LEFT(I)  =  PFSK(I) 
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WTNUM  =  N 

BOTH  =    .FALSE. 

IF (WRIT)    GOTO  5998 

LTI  = 'PHASE  OF  FSK  ENVELOPE' 

LUN  =  'PI  RADIANS' 

RTI  ■  'NONE' 

HT  =  'SAMPLE  NUMBER' 

HU  =  '  ' 

STITLE  =  '  ' 

XSCALE  =  1 

PNDM  =  N 

GOTO  5999 
C 

C 

C  AUTOCORRELATION  ESTIMATES 

C 

5400    IF(. NOT. AUTO)  THEN 

TYPE  •,  '  ' 

TYPE  »,  '  AUTOCORRELATION  ESTIMATE  NOT  ENABLED' 

GOTO  2000 
ENDIF 

TYPE  •,  '  * 

TYPE  9, '/Plot  how  many  points?  (max=',N, ')   * 

READ  (•,*)  PNUM 

DO  5410  1=1, PNUM 
LEFT(I)  =  SAUTO(I) 
RIGHT(I)  =  NAUTO(I) 
5410    CONTINUE 

WTNUM  =  PNUM 
BOTH  =  .TRUE. 
IF(WRIT)  GOTO  5998 

LTI  -'SIGNAL  AUTOCORRELATION' 

LUN  =  '  ' 

RTI  =  'NOISE  AUTOCORRELATION* 

RUN  =  '  • 

HT  =  'LAG  COMPONENT' 

HU  =  '  • 

STITLE  =  •  ' 

XSCALE  =  1 

GOTO  5999 
C 

c 

C  WEIGHT  MEAN  AND  VARIANCE  DATA 
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C 

C  PLOT  OR  WRITE  DATA  OR  DFT  OF  DATA 

C 

5500    TYPE  •,  '  ' 

TYPE  11, '/Do  yon  want  weight  sequence  or  DFT  of  the  seqnence?  (W/D) 
READ  (*,12)  C 
FFT  =  .FALSE. 

IFCC.EQ.'D*  .OR.  C.EQ.'d*)  THEN 
FFT  =  .TRUE. 

TYPE  11, '/Do  yon  want  to  window  the  seqnence?  (Y/N)   ' 
READ(*,12)  C 
IF(C.EQ.'Y'  .OR.  C.EQ.'y')  THEN 

WINDOW  =  .TRUE. 
ELSE 

WINDOW  =  .FALSE. 
ENDIF 

IF(.NOT.WRIT)THEN 
TYPE  •,  '  ' 

TYPE  11, '/Plot  np  to  what  frequency?   ' 
READ(*,*)  FREQP 
ENDIF 
ENDIF 

TYPE  *,  '  ' 
C 

C  PLOT  OR  WRITE  AVERAGE  DATA  OR  DATA  OF  ONE  TRIAL 
C 

TYPE  15 
15      FORMAT  ('  Plot  (write)  either:',/ 

1  ,'0.   Average  of  mean  and  variance',/ 

x         , '      #.   Mean  and  variance  of  one  run') 

TYPE  16  .'/Enter  selection  as  either  0  or  run  number  (np  to  ' 
X  TRIALS,'):   ' 

READ(*,*)  CHOICE 

IF(CHOICE.EQ.O)  THEN 
IF(FFT)  THEN 

TYPE*, '  ' 

TYPE*,* DFT  NOT  VALID  FOR  AVERAGE' 

GOTO  5500 
ENDIF 
DO  5510  J=0,WNUM 

LEFT(J+1)  =  WMEANA(J) 

RIGHT(J+1)    =  WVARA(J) 
5510  CONTINUE 

WTNUM  =  WNUM  +  1 

BOTH  =   .TRUE. 

IF(WRIT)    GOTO  5998 
C 

C     AVERAGE  DATA 
C 

LTT  =    'AVERAGED  MEAN' 
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urn  =  •  • 

RTI  =  'AVERAGED  VARIANCE' 

RUN  =  '  ' 

WRITE(HT,13)  'WEIGHT  MEAN  AND  VARIANCE  AVERAGE,  '.TRIALS,'  TRIALS' 

HU  =  *  ' 

PNUM  =  INTO  +  1 

XSCALE  =  1 

GOTO  5999 

ELSEIF(CHOICE.GT.O  .AND.  CHOICE. LE. TRIALS)  THEN 
C 

C  DATA  OF  INDIVIDUAL  TRIAL 
C 

IF(.NOT.FFT)  THEN 
DO  5520  J=0,WNUM 

LEFT(J+1)  =  WMEAN( J, CHOICE) 
RIGHT(J+1)  =  WVAR( J, CHOICE) 
5520        CONTINUE 

WTNUM  =  WNUM+1 

BOTH  =  .TRUE. 

IF (WRIT)  GOTO  5998 


LTI  = 

'WEIGHT  MEAN' 

LUN  = 

f   t 

RTI  = 

'WEIGHT  VARIANCE' 

RUN  = 

f    ( 

HT  =  * 

MEAN  AND  VARIANCE 

OF 

WEIGHT  VECTORS 

HU  =  * 

i 

PNUM  ■ 

WNUM  +  1 

XSCALE 

=  1 

GOTO  5999 

ELSE 

c 

c 

DFT 

OF 

DATA 

C 

DO  5530  J=0,WNUM 

AUTOC(J+l)  =  CMFLX(WMEAN(WNUM-J,  CHOICE),  0.) 
IF(WINDOW)  AUTOC(J+l)  =  AUTOC(J+l)  *  FLOAT(J)/WNUM 
5530        CONTINUE 

DO  5535  J=WNUM+1,2«WNUM 

AUTOC(J+l)  =  CMPLX(WMEAN(J-WNUM,CHOICE),0.) 

IF(WINDOW)   AUTOC(J+l)    =  AUTOCCT+l)    *    (1.0  -  FLOAT(J-WNUM)/WNUM) 
5535  CONTINUE 

DO  5540  J=2»WNUM+1,8*WNUM-1 
5540  AUTOC(J+l)   =  CMPLX(0.,0.) 

CALL  DFT(AUTOC,8*WNUM,0) 
MAXF  =  0. 

DO  5550  J=1,4*WNUM 
5550  MAXF  =  MAX(CABS(AUTOC(J) ) ,MAXF) 

DO  5560  J=1,4*WNUM 
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55«0  LEFT(J)    =  20*ALOG10(CABS(AUTOC(J))/MAXF> 

DO  5570  J=0,WNUM 

AUTOC(J+l)  =  CMFLX(WVAR(WNUM-J,  CHOICE),  0.) 
IF(WINDOW)  AUT0C(J+1)  =  AUT0C(J+1)  *  FLOAT(J) /WNUM 
5570        CONTINUE 

DO  5575  J=WNUM+1,2*WNUM 

AUTOC(J+l)  =  CMPLX(WVAR(J-WNUM,CHOICE),0.) 

IF(WINDOW)    AUTOC(J+l)    =  AUTOC(J+l)    *    (1.0  -  FLOAT(J-WNUM)/WNUM) 
5575  CONTINUE 

DO  5580  J=2*WNUM+1,8*WNUM-1 
5580  AUTOC(J+l)    =  CMPLX(0.,0.) 

CALL  DFT(AUTOC,8»WNUM,0) 
HAXF  =  0. 

DO  5590  J=1,4*WNUM 
5590  1HAXF  =  MAX(CABS(AUTOC(J) )  ,MAXF) 

DO  5595  J=1,4«WNUM 
5595  RIGHT(J)    =  20*ALOG10(CABS(AUTOC(J) ) /MAXF) 

TOJUM  ■  WNUM*4 
BOTH  =    .TRUE. 
IF(WRIT)    GOTO  5998 

LTI  =    'DFT  OF  WEIGHT  VECTOR  MEAN' 

LUN  =    '    • 

RTI  =  'DFT  OF  WEIGHT  VECTOR  VARIANCE' 

RUN  =  '  ■ 

HT  =  'FREQUENCY' 

HU  =  'HERTZ' 

XSCALE  =  FS/(WNUM*8) 

PNUM  =  FREQP/XSCALE 

GOTO  5999 

ENDIF 
ELSE 

TIPE  •,  *  ' 

TTFE  *,  '  CHOICE  NOT  WITHIN  VALID  RANGE.  TRY  AGAIN' 

GOTO  5500 

ENDIF 
C 

C 

C  INTERMEDIATE  ANC  ERROR 

C 

5600    IF(.NOT.OFILT)  THEN 

TYPE*, 'INTERMEDIATE  FILTERING  NOT  ENABLED' 

GOTO  5000 
ENDIF 
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DO  5610  1=1, N 
5610      LEFT(I)  =  ERORI(I) 

WTNUM  =  N 

BOTH  -  .FALSE. 

IF  (WRIT)  GOTO  5998 

LIT  =  'INTERMEDIATE  ANC  ERROR' 

LUN  =  '  ' 

RTI  =  'NONE' 

HT  =  'SAMPLE  NUMBER' 

HU  =  '  * 

STITLE  =  '  ' 

XSCALE  =  1 

PNUM  =  N 

GOTO  5999 
C 

c 

C  INTERMEDIATE  ADF  OUTPUT 

C 

5700    IF(.NOT.OFILT)  THEN 

TYPE*, 'INTERMEDIATE  FILTERING  NOT  ENABLED' 

GOTO  5000 
ENDIF 

DO  5710  1=1, N 
5710      LEFT(I)  =  YKI) 

FTNUM  =  N 
BOTH  =  .FALSE. 
IF(WRIT)  GOTO  5998 

LTI  -  'INTERMEDIATE  ADF  OUTPUT' 

LUN  -  '  ' 

RTI  =  'NONE' 

HT  -  'SAMPLE  NUMBER' 

HU  =  •  ' 

STITLE  =  '  ' 

XSCALE  =  1 

PNUM  =  N 

GOTO  5999 
C 

C 

C  SECONDARY  ANC  INPUT 

C 

5750    IF(.NOT.OFILT)  THEN 

TYPE*, 'INTERMEDIATE  FILTERING  NOT  ENABLED' 

GOTO  5000 
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ENDIF 

DO  5760  1=1 ,N 
5760      LEFT(I)  =  MDLD(I) 

FTNUM  =  N 
BOTH  ■  .FALSE. 
IF(IRIT)  GOTO  5998 

LTI  -  'INTERMEDIATE  ANC  INPUT' 

LBN  =  '  ' 

HTI  =  'NONE' 

HT  =  'SAMPLE  NUMBER' 

HU  =  *  ' 

STITLE  =  '  ' 

XSCALE  =  1 

PNUM  =  N 

GOTO  5999 
C 

£»»«*«»»«»»»«»••«••»»•**•«*»•» 

c 

C   SECONDARY  ANC  ERROR 

C 

5800    DO  5810  1=1, N 

5810      LEFT(I)  =  EROR(I) 

ymroM  =  n 

BOTH  =  .FALSE. 

IF  (WRIT)  GOTO  5998 

LTI  =  'ANC  ERROR' 

LUN  =  '  ' 

RTI  =  'NONE* 

HT  =  'SAMPLE  NUMBER' 

HU  =  '  ' 

STITLE  =  '  ' 

XSCALE  =  1 

PNUM  =  N 

GOTO  5999 
C 

C 

C  SECONDARY  ADF  OUTPUT 

C 

5900    DO  5910  1=1, N 

5910      LEFT(I)  =  Y(I) 

WTNUM  =  N 
BOTH  =  .FALSE. 
IF(WRIT)  GOTO  5998 
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LTI  »  'ADF  OUTPUT' 

LDN  =  '  ' 

RTI  =  'NONE' 

HT  =  'SAMPLE  NUMBER' 

HU  =  '  * 

STTTLE  =  '  ' 

XSCALE  =  1 

PNUM  «  N 

GOTO  5999 


5998    TYPE  •,  '  ' 

CALL  SGOPENdO,  'WRITE',  'Filename?   ' ,  'NONAME' ,  'REAL'  ,WTNUM) 
CALL  SGTRANdO, 'WRITE', 'REAL', LEFT.WTNUM) 


IF(BOTH)   THEN 

CALL  SGOPENdO, 'WRITE* 
CALL  SGTRANdO, 'WRITE' 

ENDIF 

GOTO  2000 


'Filename?   ', 'NONAME' , 'REAL' ,WTNUM) 
'REAL'  .RIGHT, WTNUM) 


5999    CALL  SFLOT(LEFT,RIGHT,PNUM,NN,LTI,LnN,RTI,RUN,HT,HU,TITLE, 
X  STITLE, XSCALE) 

READ(*,12)  C 
GOTO  2000 


£*•••*•«********•*•*••••••****•••***•••••••*••*••••*•••*•••••••**•*••*** 

C 

C  MODIFIT  PARAMETERS 

C 

6000    TYPE  *,  '  * 


17 


TTPE  17 

FORMAT 

' '    Do  you 

I 

,'             1. 

I 

2. 

X 

3. 

I 

'             4. 

X 

'             5. 

I 

'             6. 

x 

'             7. 

X 

8. 

X 

9. 

I 

10. 

I 

'           11. 

X 

12. 

I 

13. 

I 

'           14. 

x 

*          15. 

X 

*          16. 

Do  700  want  to  change:',/ 
Bit  sequence',/ 
FSK  frequency  shift',/ 
Signal  to  noise  ratio',/ 
Pref liter  bandwidth',/ 
Limiter  enable',/ 
Postfilter  bandwidth',/ 
ADF  delay',/ 
Number  of  weights',/ 
Gradient  constant',/ 
ADF  weight  bias  enable',/ 

ANC  pseudo-correlation  estimate  enabling',/ 
ANC  max/min  estimate  enabling',/ 
ANC  zero  crossing  estimate  enabling',/ 
ANC  output  filter  enabling',/ 
Autocorrelation  enabling',/ 
FSK  signal  information  enabling',/ 
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17.  Simulation  ADF  enabling',/ 

18.  Zero  crossing  enabling',/ 

19.  Number  of  trials  per  run',/ 

20.  Or  do  yon  want  to  make  a  run',/ 

21.  Or  do  yon  want  to  return  to  the  main  menu') 
TYPE  11,  '/Enter  seleetion:   ' 

HEAD(»,»)  CHOICE 


GOTO  (6010,6020,6030,6040,6050,6060,6070,6080,6090,6100,6110, 

6120,6130,6140,6150,6160,6170,6180,6190,1000,2000)  CHOICE 


6010 


CALL  PSEQDENCE(SEQ,BITS) 
GOTO  6000 


6020    CALL  PFSKSHIFT(FW) 
GOTO  6000 

6030    CALL  PNOISE(SNRDB,SNR) 
GOTO  6000 


6040 


CALL  PPREFILTER(PREBW,PRE) 
GOTO  6000 


6050 


CALL  PLIM(LIM) 
GOTO  6000 


6060 


CALL  PPOSTFILTER(POSTBW,POSTA,ZC,POSTBZ) 
GOTO  6000 


6070 


CALL  PADFDELAT(DELAY) 
GOTO  6000 


6080 


CALL  PWEIGHTS(WNUM) 
GOTO  6000 


6090 


CALL  PCONSTANT(MTJ) 
GOTO  6000 


6100 


CALL  PBIAS(TOIAS) 
GOTO  6000 


6110 


CALL  PCORE(SCORR) 
GOTO  6000 


6120    CALL  PHM(SMN) 
GOTO  6000 

6130  CALL  PANCZC(SANCZC) 

GOTO  6000 


6140 


CALL  PFILTHOFILT) 
GOTO  6000 
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6150    CALL  PAUTO(AUTO) 
GOTO  6000 

6160    CALL  PFSKINFO(FSK) 
GOTO  6000 

6170    CALL  PADF(ADF) 
GOTO  6000 

6180  CALL  PZCROSS(ZC) 

GOTO  6000 

6190  CALL  PTRIALS(TRIALS) 

GOTO  6000 
C 

C 

C  ENTER  USER  CONTROLLED  ADAPTIVE  FILTERTNG  ROUTINE 

C 

7000    DO  7100  1=1, N 

MDLD(I)  =  REAL(S(D) 
7100    CONTINUE 

CALL  SADF(MDLD,DLDNN,N,TS, TITLE) 

GOTO  2000 
C 

£•*•»»»•••*••••••••«»•••••»•»•••*••••••*••«•**••«»•••»••**•*•*»••••»•••• 

c 

C  WRITE  PARAMETERS,  STATISTICS,  AND  ESTIMATES 

C 

8000    TYPE  •, *  * 

TYPE  24 
24      FORMAT  ('  Do  yon  want  to  write:',/, 

x  '      1.   Simulation  parameters',/, 

x         '     2.  Estimate  error  statistics',/, 
x  '     3.   Band  rate  estimates',/, 

x         '4.  Everything') 
TYPE  11,  '/Enter  selection:   ' 
READ(*,»)  CHOICE 
C 

C  NUMBER  OF  ESTIMATES  TO  WRITE 
C 

IF(CHOICE.EQ.3  .OR.  CHOICE. EQ. 4)  THEN 
TYPE  •,'  ' 

TYPE  25,'  There  are  '.TRIALS,'  estimates' 

TYPE  11, '/Enter  first  estimate  and  number  of  estimates:   ' 
READ(*,»)  C1.C2 
ENDIF 

GOTO  (8100,8200,8300,8100)  CHOICE 
C 

£•*•••••••*««««***•• ********** 
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C 

C  PARAMETERS 
C 

8100    WRITE'*, 26)  TRIALS, SNRDB.N.FS, BITS, R,FC,FW,2*PREBW 

26      FORMATC   ',/./, 

X         '  SIMULATION  PARAMETERS:',/, 
x       /,'      Simulation  Parameters:',/, 
x  Number  of  trials  per  run  =   ',13,/, 

x         '         Signal-to-noise  ratio  =   '.F4.1,'  dB',/, 
I       /,'      Sampling  Parameters:',/, 
x  Number  of  samples  =   ',14,/, 

x  Sampling  Frequency  ■   ',F6.0,'  Hz',/, 

X       /,'     Bit  Sequence  Parameters:',/, 
x         '         Number  of  bits  =   ',12,/, 
x         '         Baud  rate  =   ',F5.0,'  Hz',/, 
x       /,'     FSK  Modulation  Parameters: ',/, 
x         '         Carrier  frequency  =   ',F8.0,'  Hz',/, 
x  Frequency  shift  =   ',F5.0,'  Hz',/, 

x       /,'     Receiver  Parameters: ',/, 
x         '         Pre-filter  bandwidth  =   ',F5.0,'  Hz') 
IF(.NOT.POSTA)  THEN 

WRITE'*,*)  '        Post  filtering  not  enabled  for  ANC 
IF(ZC)  WRITE(*,31)  POSTBZ 

31  FORMATC         Post-filter  bandwidth  for  ZC  =   '.F5.0,'  Hz') 
ELSE 

WRITE(*,32)  POSTBW 

32  FORMATC         Post-filter  bandwidth  =   ',F5.0,'Hz') 
ENDIF 

IF(LIM)  THEN 

WRITE(*,*)  *        Limiter  enabled' 
ELSE 

WRITE(*,*)  '        Limiter  not  enabled' 
ENDIF 

WRITE(*,27)  DELAT.MtJ.WNUM 
27      FORMATC  ',/, 

X  Adaptive  Noise  Cancelling  Parameters:',/, 

X         '         Reference  input  delay  =   ',12,/, 
x         '         Gradient  constant  =   ',E8.1,/, 
x         *         Number  of  weights  =   ',13) 

IF(WBIAS)  THEN 

WRITE(*,*)  '        Weighting  bias  included* 
ELSE 

WRITE'*,*)  •        Weighting  bias  not  included' 
ENDIF 

WRITE(*,»)  '  * 

WRITE(*,»)  •  • 

IF  (ADF  .OR.  ZC)  WRITE  (*,*)  'ESTIMATORS  ENABLED:' 

WRITE(*,*)  '  ' 

IF(OFILT)  WRITE(*,*)  •     ANC  Intermediate  Filtering  enabled' 

IF(SMM)  WRITEC*,*)  •     Min/Max  estimators  enabled' 
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IF(SANCZC)  WRITEC*,*)  '     ANC-ZC  estimators  enabled1 

IF(SCORR)  WRITEC*,*)  *     Pseudo-Correlation  estimators  enabled' 

IF(ZC)  WRITEC*,*)  '     Zero  Crossing  estimator  enabled' 

IFCCHOICE.EQ.l)  GOTO  2000 
C 

C 

C  WRITE  STATISTICS 

C 

WRITEC*, 28) 
8200    WRITEC*,*)  '  ' 
WRITEC*,*)  *  ' 

WRITEC*,*)  'NORMALIZED  BAUD  RATE  ESTIMATE  ERROR  STATISTICS:' 
WRITEC*,*)  '  ' 
C 

C   INTERMEDIATE  ESTIMATE  STATISTICS 
C 

IF  (ADF)  THEN 
WRITEC*,*)  '  * 

WRITEC*, •)  '  MMI    MM2I   AMI   MLTOI   MZCI   MAI    VAI 
X  '  EXCI' 

WRITEC*,*)  '  ' 

WRITEC*, •)  'Normalized  band  rate  mean  error:* 
WRITE'*, 19)  (SIESTM(I),I=1,8) 

WRITEC*,*)  'Normalized  band  rate  error  variance:' 
WRITE'*, 19)  (SIESTV(I),I=1,8) 

WRITE(*,*)  'Normalized  band  rate  mean  squared  error:' 
WRITE'*, 19)  (SIESTS(I),I=1,8) 
ENDIF 


c 

c 
c 

SECONDARY  ESTIMTATE  STATISTICS 

IF(OFTLT)  THEN 

WRITE'*, •)  '  ' 

WRITE'*,*)  '  MMF    MM2F   AMF 

X 

'  XZCF   TZCF   EZCF 
WRITE(*,*)  '  • 
WRITE'*,*)  'Normalized  baud  rate 

MLTOF   MZCF   MAF    VAF 


WRITE(*,19)  (SFESTM(I),I=1,10) 

WRITEC*,*)  'Normalized  baud  rate  error  variance:' 
WRITE(*,19)  (SFESTV(I),I=1,10) 

WRITEC*,*)  'Normalized  baud  rate  mean  squared  error:' 
WRITE(*,19)  (SFESTS(I),I=1,10) 
C 

C  BANDWIDTH  STATISTICS 
C 

IF(OFILT)  THEN 
WRITEC,*)  •  ' 

WRITEC*,*)  'BANDWIDTH  ESTIMATE:  ' 
WRITE(*,*)  '  ' 
WRITEC*, 30)  '  Bandwidth  mean  =  ',SBWM, 
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I  '    Bandwidth  standard  deviation  =  ' , SBWV 

WRITEt*, 33)  '  Normalized  bandwidth  mean  =  '  ,NSBWM, 
X  '     Standard  deviation  ■  '.NSBWS 

ENDIF 
ENDIF 
C 

C  TRADITIONAL  ZERO  CROSSING  STATISTICS 
C 

IF  (ZC)  THEN 
WRITEt*,*)  '  ' 

WRITEt*,*)  'ZERO  CROSSING  ESTIMATE  ERROR  STATISTICS:* 
WRITEt*,*)  'Normalized  band  rate  mean  error:' 
WRITE(*,19)  fZESTM(I), 1=1.1) 

WRITEt*,*)  'Normalized  band  rate  error  variance:' 
WRITEt*, 19)  (ZESTV(I). 1-1,1) 

WRITEf*,*)  'Normalized  band  rate  mean  squared  error:' 
WRITEt*, 19)  (ZESTS(I), 1=1,1) 
ENDIF 
WRITEt*, *>  '  * 

IFCCH0ICE.EQ.2)  GOTO  2000 
C 

£••*••••••*••••••*•••••••**•*• 

c 

C  WRITE  ESTIMATES 
C 

WRITEt*, 28) 
8300    WRITEt*,*)  '  ' 
WRITEt*,*)  '  ' 

WRITEt*,*)  'BIT  LENGTH  ESTIMATES  (IN  NUMBER  OF  SAMPLES):' 
WRITEt*,*)  '  ' 
C 

C  INTERMEDIATE  ESTIMATES 
C 

IF  (ADF)  THEN 
WRITEt*,*)  *  ' 

WRITEt*,*)    '      MMI  MM2I        AMI        MLTOI        MZCI        MAI  VAI 

X  '   EZCIT' 

WRITEt*,*)    '    ' 
DO  8310   II=C1,C1+C2-1 

WRITEt*, 29)    (SIEST(I,II),I=1,8) 
IF(MOD(II,10).EO.O)    WRITEt*,*)    '    ' 
8310  CONTINUE 


ENDIF 

c 

c 

SECONDARY  ESTIMATES 

C 

IFtOFILT)  THEN 
WRITEt*,*)  *  ' 

WRITEt*,*)  '   MMF 

MM2F 

AMF 

X 

'  XZCF 
WRITEt*,*)  '  ' 

TZCF 

EZCF 

MLTOF        MZCF       MAF  VAF 
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DO  8340  II=C1,C1+C2-1 

WRTTE(*,29)  (SFESTd, II), 1-1,10) 
IF(HOD(II,10).EQ.O)  WRITE(»,»)  '  ' 
8340      CONTINUE 

ENDIF 
C 

C  TRADISTION  ZERO  CROSSING  ESTIMATES 
C 

IF(ZC)  THEN 
WRITEt*,*)  *  ' 

WRITEC,*)  'Zero  Crossing  Estimates:' 
WRITE(«,»)  '  ' 
DO  8360  II-C1.C1+C2-1 

WRITE(»,29)  (ZESTd, II), 1=1,1) 
IF(MOD(II,10).EQ.0)  VRITE(»,«)  '  ' 
8360      CONTINUE 

WRITE(*,»)  '  ' 
ENDIF 

GOTO  2000 


C 

£•*•••*»*•****•*••******•***•*«**••••••••••**»*••••**•••••••**••***•*••• 

C 

C  FORMAT  STATEMENTS  AND  MAJOR  PROGRAM  ENDING 

C 

9000    CONTINUE 

9  FORMAT(A28,I4,A5) 

10  F0RMAT(A80) 

11  FORMAT(A) 

12  FORMAT  (Al) 

13  FORMAT(A34,I3,A7) 

14  FORMAT(A10,I2,A3,I3,A5,F12.10,A6) 
16  FORMAT(A50.I3,A4) 

18  FORMAT  (A2  9, 14,  A3) 

19  FORMATC  ',10(F5.2,2D) 

20  FORMATC  '  ,10(2X,I3,2X) ) 
25  F0RMAT(A11,I3,A11) 

28  FORMATC1*) 

29  FORMATC  '  ,10(F6.1,1X) ) 

30  FORMAT(A18,F7.0,A35,F7.0) 
33  FORMAT(A27,F5.2,A30,F5.2) 

STOP 
END 

C 

c 

C  PARAMTER  INPUT  SUBROUTINES 
C 

SUBROUTINE  PRMTRS 
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REAL  POSTBW.POSTBZ.SNRDB 

INTEGER  SEQ(64), BITS, BIT,N,N1,ND, TRIALS 

CHARACTER*  1  C 

LOGICAL  AUTO,ADF,FSK,SNR,FRE,POSTA,POSTZ,ZC,ZC4,LIM 
C 

<;•••«»•****••••••**•*«•••••*»• 

c 

C   SEQUENCE 
C 

ENTRY  PSEQDENCE(SEQ,BITS) 

TTPE  •,  '  ' 

TYPE  11,' f         Test  sequence  or  pseudo-random  sequence  (T/S)? 

READ  (*,12)  C 

THEN 


100 


IF    (C.EQ.'T'    .OR.    C.EQ. 

't') 

DO  100   1=  1, BITS/2 

SEQ(2*I)    =  1 

SEQ( 2*1-1)    =  0 

CONTINUE 

ELSEIF   (C.EQ.'S'    .OR.    C. 

,m. 

DO  150   1=0,32,32 

SEQU+1)    =  1 

SEQ(I+2)    =  0 

SEQU+3)    -  0 

SEQ(I+4)    -  1 

SEQ(I+5)    -  1 

SEQU+6)    =  0 

SEQ(I+7)    ■  1 

SEQ(I+8)    =  0 

SEQ(I+9)    =  0 

SEQ(I+10)    =  1 

SECKI+11)    =  0 

SEQU+12)    =  0 

SEQ(I+13)    -  0 

SEQU+14)    »  0 

SEQ(I+15)    -  1 

SEQ(I+16)    =  0 

SEQ(I+17)    -  1 

SEQ(I+18)    =  0 

SEQU+19)    =  1 

SEQ(I+20)    =  1 

SEQ(I+21)    =  1 

SEQU+22)    =  0 

SEQ(I+23)    -  1 

SEQ(I+24)    =  1 

SEQ(I+25)    =  0 

SEQU+26)    =  0 

SECKI+27)    =  0 

SEQ(I+28)    =  1 

SEQ(I+29)    =  1 

s')  THEN 
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SEQ(I+30)  =  1 
SEQU+31)  -  1 
SEQ.U+32)  =  1 
150       CONTINUE 
ELSE 

DO  200  1=1, BITS 
200         SEQ(I)  =  0 
ENDIF 

RETURN 
C 

c 

C  NOISE 
C 

ENTRY  PNOISE(SNRDB,SNR) 

TYPE  *,  '  ' 

TYPE  11,  '/    Do  yon  want  to  add  noise  (Y/N) ?    ' 

READ(«,12)  C 

IF  (C.EO.'Y'  .OR.  C.EQ. 'y')  THEN 

SNR  "  .TRUE. 

TYPE  11,  '/    Signal  to  noise  ratio  (dB)?    ' 

READ  (*,*)  SNRDB 
ELSE 

SNR  =  .FALSE. 
ENDIF 

RETURN 
C 
£•*•••*•••••**•••*•••••••••••• 

C 

C  PREFILTER 

C 

ENTRY  PPREFILTER(PREBW3,PRE) 

TYPE  *,  '  ' 
C       TYPE  11,  '/    Do  yon  want  to  perform  pre-f iltering? 
C       READ  (*,12)  C 
C  =  'Y' 
IF  (C.EQ.'Y'  .OR.  C.EQ.'y')  THEN 

PRE  =  .TRUE. 

TYPE  11,  '/    Pre-filter  bandwidth?    ' 

READ  (*,•)  RFBW 

PREBW3  =  RFBW/ 2. 
ELSE 

PRE  =  .FALSE. 

PREBW3  =  0.0 
ENDIF 

RETDRN 
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C 
(;•••«••••*•*••••»*«»•••••••••* 

C 

C     POST  FILTER 

C 

ENTRY  PPOSTFILTER(POSTBW,POSTA,ZC4, POSTBZ) 

TYPE  •,  '  ■ 

TYPE  11,  '/    Do  yon  want  to  perform  post-filtering?   ' 

READ  (»,12)  C 

IF  (C.EQ.'Y'  .OR.  C.EQ. 'y')  THEN 

POSTA  =  .TRUE. 

TYPE  11,  '/    Post-filter  bandwidth? 

HEAD  (*,*)  POSTBW 

POSTBZ  =  POSTBW 
ELSE 

POSTA  =    .FALSE. 

IF(ZC4)    THEN 

TYPE  11,  '/    Post-filter  bandwidth  for  ZC  algorithm?    ' 
HEAD  (*,•)  POSTBZ 

ENDIF 
ENDIF 

RETDRN 
C 

C 

C  LIMITER 

C 

ENTRY  PLIM(LIM) 

TYPE  »,  '  * 

TYPE  11,  '/    Do  yon  want  to  enable  the  limiter?   ' 

READ(«,12)  C 

IF(C.EQ.'Y'  .OR.  C.EQ. 'y')  THEN 

LIM  =  .TRUE. 
ELSE 

LIM  =  .FALSE. 
ENDIF 

RETDRN 
C 

C 

C  DLD  DELAY 

C 

ENTRY  PDELAY(ND) 

TYPE  *,  *  * 

TYPE  11,  '/    Delay  factor  of  DLD  (an  integer)?    ' 

READ(*,»)  ND 
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RETURN 
C 

C 

C  FSK  FREQUENCY  SHIFT 

C 

ENTRY  PFSKSHIFT(FW) 

TYPE  *,  '  ' 

TYPE  11  ,'/    FSK  frequency  shift? 

READ  (•,•)  FW 

RETURN 
C 
c.«. ..**.,»,*,*.,.**.. .„,,.., 

C 

C  NUMBER  OF  TRIALS 

C 

ENTRY  PTRIALS  (TRIALS) 


TYPE  *,  '  ' 

TYPE  11,  '/Number  of  trials  per  run?    ' 

READ  (*,•)  TRIALS 

RETURN 
C 

c 

C  TRADITIONAL  ZERO  CROSSING  ENABLE 
C 

ENTRY  PZCROSS(ZC) 

TYPE  *,  '  ' 

TYPE  11,  '/Traditional  zero  crossing  algorithm  enabled?   (Y/N) 

READC12)  CANS 

IF  (CANS.EQ.'Y'  .OR.  CANS.EQ.'y')  THEN 

ZC  =  .TRUE. 
ELSE 

ZC  =  .FALSE. 
ENDIF 

RETURN 
C 

C 

C  ADAPTIVE  NOISE  CANCELLING  ENABLE 

C 

ENTRY  PADF(ADF) 

TYPE  *,  '  * 
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TYPE  11,  '/Adaptive  noise  cancelling  enabled?   (T/N)   ' 

READ(*,12)  CANS 

IF  (CANS.EQ.  'T»  .OR.  CANS.EQ.'y')  THEN 

ADF  =  .TRUE. 
ELSE 

ADF  =  .FALSE. 
ENDIF 

RETURN 
C 

C****»f  ********************* 

c 

C  AUTOCORRELATION  ENABLE 
C 

ENTRY  PAUTO(AUTO) 

TYPE  •,  '  ' 

TYPE  11,  '/Autocorrelation  enabled?   (Y/N)   ' 

READ(«,12)  CANS 

IF  (CANS.EQ.'Y'  .OR.  CANS.EQ.'y')  THEN 

AUTO  =  .TRUE. 
ELSE 

AUTO  =  .FALSE. 
ENDIF 

RETURN 
C 

c***************************** 

C 

C  FSK  SIGNAL  INFORMATION  ENABLE 

C 

ENTRY  PFSKINFO(FSK) 

TYPE  *,  '  ' 

TYPE  11,  '/Phase  and  frequency  of  FSK  signal  enabled?  (Y/N)   ' 

READ(*,12)  CANS 

IF  (CANS.EQ.'Y'  .OR.  CANS.EQ.'y')  THEN 

FSK  =  .TRUE. 
ELSE 

FSK=  .FALSE. 
ENDIF 

RETURN 

11  FORMAT(A) 

12  FORMAT  (Al) 

END 
C 
C****************************************************»*****tt*********** 

c 

C  NOISE  GENERATION  SUBROUTINE 
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SUBROUTINE  SNOISE(CNOISE,N,  A.TS.PREBW,  SNRDB, LSNR) 

COMPLEX  CNOISE(2048) 

REAL  A.TS.PREBW, SNRDB,  VAR,  MEAN, SNR,  RNOISE, INOISE 

rNTEGER  N 

LOGICAL  LSNR 

IF(. NOT. LSNR)  THEN 

DO  300  1=1, N 
300         CNOISE(I)  =  CMPLX(O.,0.) 
ELSE 

SNR  =  10.**(SNRDB/10.0) 

VAR  =  (A*»2/2)/(TS*PREBW»SNR) 

MEAN  -  0.0 

DO  400  1=1,  N 

RNOISE  -  GAUSS (MEAN, VAR) 
INOISE  =  GAUSS(MEAN.VAR) 
CNOISE(I)  =  CMPLX (RNOISE, INOISE) 
400       CONTINUE 
ENDIF 

RETURN 

END 
C 
<;*****»»*••**•*•••»*«*•••*«»*«••••*«•***»«»**•«««»*«*»«»»»»*,««««««„»», 

c 

C  PLOT  TITLE  FORMATION  SUBROUTINE 
C 

SUBROUTINE  PTITLE(TITLE,PRE,POSTA,FW,PREBW,POSTBW,  SNRDB,  SNR.R.LIM) 

IMPLICIT  NONE 

REAL  FW,PREBW,POSTBW,R,  SNRDB 
CHARACTER  TITLE*80,PRETI*13,POSTTI*14 
CHARACTER  STI*14,FTI*12,SNRTI*10,LTI*9 
LOGICAL  SNR,PRE,POSTA,LIM 

WRITE    (PRETI.19)    ',    PRE-BW=* ,IFIX(PREBW»2. ) 

19  FORMAT(A,I4) 

WRITE(POSTTI,20)    ',    POST-BW=',    IFIX(POSTBW) 

20  FORMAT  (A, 14) 

TOITE(STI,21)  'BAUD  RATE=  * , IFIX (R) 

21  FORMAT(A10,I4) 

WRITE(FTI,23)  »,  SHIFT=' ,IFIX(FW) 

23  FORMAT(A,I4) 

WRITE  (SNRTI.24)  ',  SNRDB=* ,  IFIX  (SNRDB) 

24  FORMAT(A,I2) 

IF(.NOT.SNR)  SNRTI  =  • ,  NO  NOISE' 
IF(LIM)  THEN 
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LTI  =    ', 

LIMITER 

LTI  =    * 

1 

IF  (PRE)  THEN 
IF(POSTA)  THEN 

TITLE  =  STI//FTI//PRETI//POSTTI//SNRTI//LTI 
ELSE 

TITLE  =  STI//FTI//PRETT//',  NO  POST' //SNRTI/ /LTI 
ENDIF 
ELSE 

IF(POSTA)  THEN 

TITLE  =  STI//FTI//',NO  PRE '//POSTTI// SNRTI/ /LTI 
ELSE 

TITLE  =  STI//FTI//',NO  PRE,  NO  POST* //SNRTI/ /LTI 
ENDIF 
ENDIF 

RETURN 
END 
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c 

C  SADF 

C 

C  VAX-11  FORTRAN  SOURCE  FILENAME:  SADF1.FOR 

C 

C       DEPARTMENT  OF  ELECTRICAL  ENGINEERING    KANSAS  STATE  UNIVERSITY 

C 

C       REVISION        DATE  PROGRAMMER(S) 

C       


C        1.0  MAT  6,  1985  MARC  D.  BRACK 

C 

C********************************************************************** 

c 

C       CALLING  SEQUENCE 

C 

C  CALL   SADFU.XNN, LENGTH, TS, TITLE) 

C 

C  PURPOSE 

C 

C  Implements  the  adaptive  noise  canceller  in  a  nser 

C  controlled  environment 

C 

C       ROUTINE  (S)  ACCESSED  OR  CALLED  BY  THIS  ROUTINE 

C 

C  SPLOT 

C  BUTER 

C  DFT 

C 

C       ROUTINE (S)  ACCESSED  OR  CALLED  BY  EMBEDDED  ROUTINE(S) 

C 

C  NONE 

C 

C  ARGUMENT(S)    REQUIRED  FROM  THE  CALLING  ROUTINE 

C 

C  X       Ontpnt  of  the  delay-line  discriminator  (DLD) 

C 

C  XNN     Output  of  the  DLD  when  the  input  was  noise  free 

C 

C  LENGTH  Number  of  samples  in  DLD  output 

C 

C  TS      Sampling  interval 

C 

C  TITLE   Major  title  for  plotting  purposes 

C 

C       ARGUMENT(S)  SUPPLIED  TO  THE  CALLING  ROUTINE 

C 

C  NONE 

C 

C 

SUBROUTINE  SADF(X,XNN,LENGTH,TS. TITLE) 
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IMPLICIT  NONE 

INTEGER  N, I. NUM1, LENGTH, DELAY,N2,NN,PNUM,P2,TN,TDELAY,TMIJI 

INTEGER  PFIRST, CHOICE, CS1.CS2, J, N3 

REAL  MD,W(2049, 0:256)  ,WMEAN(0:256)  ,WAR(0:256)  ,WMEANS(0:256) 

REAL  EROR(2048) ,LEFT(2048) ,RIGHT(2048) ,1(2048) ,X(2048) 

REAL  FEROR(2048) ,FT(2048) ,BW,MAXF,FX(2048) 

REAL  WDFT(4096)  ,TMU,TS,FREQP,XNN(2048)  .XSCALE 

COMPLEX  CW(2048),CW2(2048),BUTER 

CHARACTER*80  NAME1,NAME2,BNK,NAME,LT,RT,HT,TITLE,STITLE,HU 

CHARACTER  DELAYC»4 , NC*4 , MUC*11 , WAIT'l , CANS*1 , CANS2  «1 

LOGICAL  TOLAS, FILT 

C 

C  INPtTT  PARAMETERS  FOR  ROUTINE 

C 

CALL  PADFDELAY  (DELAY) 

CALL  PWEIGHTS(N) 

CALL  PCONSTANT(MU) 

CALL  PBIAS(TOIAS) 

CALL  PFILTER(FILT.BW) 
C 
C***********************************<* •••••••••*•••••*••••••••••••••*•• 

C 

C  ADAPTIVE  NOISE  CANCELLER  SECTION,  LMS  ALGORITHM 

C 

1000  CONTINUE 

DO  100   1=0, N 
WMEAN(I)    =  0. 
WMEANS(I)    =  0. 
WVAR(I)    =  0. 
100  CONTINUE 

DO  120  J=1,LENGTH+1 
DO  110   1=0, N 
110  W(J,I)    =  0.0 

120  CONTINUE 

DO  160  J=l, LENGTH 
IF    (TOIAS)    THEN 

Y(J)    =  W(J,0) 
ELSE 

Y(J)    =  0.0 
ENDIF 

DO  130  1=1, N 
130  IF(J-I+l.GT.O)   Y(J)    =  Y(J)    +  X(J-I+1) •W(J,I) 

IF(J+DELAY.LE.LENGTH)    THEN 

EROR(J)    =  XU+DELAY)    -  Y(J) 
ELSE 
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EROR(J)  =  -Y(J) 
ENDIF 

W(J+1,0)  =  W(J,0)  +  2.»MU»EROR(J) 

DO  140  1=1, N 
140         IF(J-I+l.GT.O)  W(J+1,I)  =  W(J,I)  +  2.»MU»EROR(J)*X(J-I+l) 

IF(J.GT.N)  THEN 
DO  150  1=0, N 

WMEAN(I)  =  WMEAN(I)  +  W(J+1, 1) /LENGTH 
WMEANS(I)  =  WMEANS(I)  +  W(J+1,I)»»2/LENGTH 
150         CONTINUE 

ENDIF 
160     CONTINUE 

DO  170  1=0, N 
170       WVAR(I)  =  WMEANS(I)  -  WMEAN(I)*»2 

C 

C  FILTER  THE  ERROR  AND  THE  FILTER  OUTPUT 

C 

IF  (FILT)  THEN 
DO  180  1=0, N 
180         Cff(I+l)  =  CMPLX(WMEAN(I),0.) 

DO  190  I=N+1,4»N-1 
190         CWU+1)  =  CMPLX(0.,0.) 
C 

C  DETERMINE  3  dB  BANDWIDTH  OF  SIGNAL 
C 

CALL  DFT(CW,4»N,0) 

MAXF  =  0 
DO  200  I=l,4«N/2 
200         IF(CABS(CW(I)).GT.MAXF)  MAXF  =  CABS(CWd)) 
MAXF  =  MAXF**2 

DO  210  I=l,4*N/2 

IF(CABS(CW(4*N/2+l-I))*»2.GT.MAXF/2)  THEN 
BW  =  (4*N/2+2-I)/(4«N/2*TS) 
GOTO  220 
ENDIF 
210       CONTINUE 
220       CONTINUE 

TYPE  •,  'BANDWIDTH  =   ',BW 
ENDIF 
C 

C 

C  PLOTTING  AND  PARAMETER  CHANGE 


c 

3000  WRITE(»,3) 

3  FORMATC  Do  yon  want  to  plot:',/ 

x  , '     1 .  DLD  ontpnt ' , / 

x  ,'    2.  Ontpnt  error'./ 

x  »'    3.  Filter  estimate',/ 

i  , *    4.  Weights  of  filter  for  a  specific  vector',/ 

x  ,'    5.  Frequency  response  of  filter',/ 

x  ,'6.  Mean  and  variance  of  weights',/ 

x  ,'7.  History  of  a  specific  weight*,/ 

x  ,'  Or  do  yon  want  to:',/ 

x  , '    8.  Change  the  input  delay',/ 

x  ,'    9.  Change  the  gradient  constant',/ 

x  ,'   10.  Change  the  number  of  weights',/ 

x  ,'   11.  Change  the  weighting  bias'./ 

x  ,'12.  Write  the  results  to  disk',/ 

x  ,'   13.  Change  the  automatic  filter  enable',/ 

x  ,'14.  Perform  a  run',/ 

x  ,'15.  Return  to  the  calling  routine') 

TYPE  1,  '/Enter  selection:    ' 

HEAD  (*,•)  CHOICE 

IF  (CHOICE. LE.  7)  THEN 

1  FORHAT(A) 

2  FORHAT(Al) 
WRITE(»,»)  *  • 

TYPE  1,  '/First  point  to  plot  =   ' 

READ  (•,•)  PFIRST 
C         TYPE  1,  '/interval  between  points  =   ' 
C         READ  (•,•)  NN 
ENDIF 

NN  =  1 

GOTO  (3100,3150,3200,3300,3400,3500,3600,4000, 
X        4100,4200,4300,4400,4500,1000,9000),  CHOICE 
C 

C****»»***«*»*********«*****»»»»*****»««***»»«*««»»»«**e*»***««»»«****» 
C 

C  PLOTTING 
C 

C  PLOT  DLD  OUTPUT 
C 

3100    TYPE  1,  '/Number  of  points  to  plot?   ' 
READ  (•,»)  PNUM 
IF  (FILT)  THEN 

TYPE  1,  '/Plot  DLD  output  or  filtered  DLD  output?  (D/F)   ' 
READ  (»,2)  CANS 
ELSE 

CANS  =  'D' 
ENDIF 
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C     PERFORM  INTERMEDIATE  FILTERING 
C 

IF   (CANS.EQ.'F'    .OR.    CANS.EQ. '  f ')    THEN 
DO  290  J=l. LENGTH 
290  CT(J)    =  CMPLX(X(J).0.) 

CALL  DFT(CW,LENGTH,0) 
DO  300  J-l,  LENGTH/  2 

CW(J)    -  CW(J)*BTJTER(4,BW,FL0AT(J-l)/(TS*LENGTfl)) 
CW(LENGTH+1-J)   =   CW(LENGTH+1-J)  * 
X  BUTER(4,BW,FL0AT(-J)/(TS»LENGTH)) 

300  CONTINUE 

CW( LENGTH/ 2+1)    =  CMFLX(0.,0.) 
CALL  DFT(CW,  LENGTH,  1) 

DO  310  J-l, LENGTH 
310  FXCT)      =  REAL(CW(J)) 

ENDIF 

DO  3110   I=PFIRST, LENGTH 

IF  (I+DELAT.LT.  LENGTH)    THEN 

RIGHT (I-PFIRST+1)    =  XNN(I+DELAT) 
ELSE 

RIGHT (I-PFIRST+1)    =0.0 
ENDIF 
IF    (CANS.EQ.'F'    .OR.    CANS.EQ.  'f')   THEN 

LEFT  (I-PFIRST+1)    =  FX(I) 
ELSE 

LEFT( I-PFIRST+1)    =  X(I) 
ENDIF 
3110  CONTINUE 

RT  =    'NOISESLESS  SIGNAL' 

IF   (CANS.EQ.'F'. OR. CANS.EQ. 'f')   THEN 
LT  =    'FILTERED  OLD  OUTPUT' 

ELSE 

LT  =  *DLD  OUTPUT' 

ENDIF 

HT  =  'SAMPLE  NUMBER* 

HU  =  '  • 

XSCALE  ■  1 
C       PNUM  =  LENGTH  -  PFIRST  +  1 

GOTO  3900 
C 

C  PLOT  ERROR  OF  ANC 
C 
3150    TYPE  1,  '/Number  of  points  to  plot?   * 

READ  (*,•)  PNUM 

IF  (FILT)  THEN 

TYPE  1,  '/Plot  error  or  filtered  error?  (E/F)   ' 
READ  (*,2>  CANS 

ELSE 
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CANS  =    'E' 
ENDIF 
C 

C     PERFORM  INTERMEDIATE  FILTERING 
C 

IF   (CANS.EQ.'F'    .OR.    CANS.EQ.'f')   THEN 
DO  260  J-l, LENGTH 
260  CW(J)    =  CMPLX(EROR(J),0.) 

CALL  DFT(CW,LENGTH,0) 
DO  270  J-l. LENGTH/ 2 

CT(J)    =   CW(J)*BUTER(4,BW,FLOAT(J-l)/(TS»LENGTH)) 
CW(LENGTH+1-J)   =   CT(LENGTH+1-J)« 
X  BUTER(4,BW,FLOAT(-J)/(TS«LENGTH)) 

270  CONTINUE 

Cf (LENGTH/2+1)    =  CMPLX(0.,0.) 
CALL  DFT (CW, LENGTH,  1) 

DO  280  J-l, LENGTH 
280  FEROR(J)      =  REAL(Ct(J)) 

ENDIF 

DO  3160  I-PFIRST, LENGTH 

IF(I+DELAT.LT.LENGTH)    THEN 

RIGHT  (I-PFIRST+1)    =  XNNU+DELAY) 
ELSE 

RIGHT (I-PFIRST+1)    =0.0 
ENDIF 
IF    (CANS.EQ.'F'    .OR.    CANS.EQ.'f)   THEN 

LEFT (I-PFIRST+1)    =  FEROR(I) 
ELSE 

LEFTU-PFIRST+1)  =  EROR(I) 
ENDIF 
3160    CONTINUE 

RT  =  'NOISESLESS  SIGNAL* 

IF  (CANS.EQ.'F'.OR.CANS.EQ.'f)  THEN 
LT  =  'FILTERED  ERROR' 

ELSE 

LT  =  'ERROR' 

ENDIF 

HT  -  'SAMPLE  NUMBER' 

HU  =  *  ' 

XSCALE  =  1 
C       PNOM  =  LENGTH  -  PFIRST  +  1 

GOTO  3900 
C 

C  PLOT  ADF  OUTPUT 
C 
3200    TYPE  1,  '/Number  of  points  to  plot?   ' 

READ  (•,•)  PNUM 

IF  (FILT)  THEN 
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TYPE  1,  '/Plot  filter  output  or  filtered  filtered  output?  (0/F) 
READ  (*,2)  CANS 
ELSE 

CANS  =  *0' 
ENDIF 
C 

C  PERFORM  INTERMEDIATE  FILTERING 
C 

IF    (CANS.EQ.  'F'    .OR.    CANS.EQ. 'f')    THEN 
DO  230  J-l. LENGTH 
230  CT(J)    =  CMPLX(Y(J),0.) 

CALL  DFT(C¥, LENGTH, 0) 
DO  240  J-l, LENGTH/2 

CW(J)    =   Clf(J)*BUTER(4,BW,FL0AT(J-l)/(TS»LENGTH)) 
CWLENGTH+1-J)   =  Cf  (LENGTH+1-J)* 
I  BUTER(4,BW,FL0AT(-J)/(TS*LENGTH)) 

240  CONTINUE 

CIT(LENGTH/2+l)    =  CMPLX(0.,0.) 
CALL  DFT(CTf,LENGTH,l) 

DO  250  J-l, LENGTH 
250         FT(J)   -  REAL(CWU)) 
ENDIF 

DO  3210  I- PFIRST,  LENGTH 

IF  (I+DELAT.LT.  LENGTH)    THEN 

RIGHT  (I-PFIRST+1)    =  XNNU+DELAY) 
ELSE 

RIGHT(  I-PFIRST+1)    =0.0 
ENDIF 
IF   (CANS.EQ. 'F'    .OR.    CANS.EQ. 'f')   THEN 

LEFT (I-PFIRST+1)    =  FT (I) 
ELSE 

LEFT (I-PFIRST+1)  =  T(I) 
ENDIF 
3210    CONTINUE 

RT  -  'NOISELESS  SIGNAL' 

IF  (CANS.EQ. 'F'  .OR.  CANS.EQ. 'f')  THEN 
LT  =  'FILTERED  FILTER  OUTPUT' 

ELSE 

LT  =  'FILTER  OUTPUT' 

ENDIF 

HT  =  'SAMPLE  NUMBER' 

HU  =  '  ' 

XSCALE  =  1 
C  PNUM  =  LENGTH  -  PFIRST  +  1 

GOTO  3900 
C 

C  PLOT  A  SPECIFIC  WEIGHT  VECTOR 
C 
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3300    IF(PFIRST.GT.N)  THEN 

WRITEC*,*)  'ERROR:   THERE  ARE  NOT  THAT  MANY  WEIGHTS* 
GOTO  3000 
ENDIF 

WRTTE(*,*)  '  ' 

TYPE  1,  '/Weights  for  what  vector  number?   (a,b)   * 
READ  (•,*)  CS1.CS2 
IF  (CS2.GE.1)  THEN 
DO  3310  1=0, N 

LEFTU-PFIRST+2)  =  W(CS1,I-PFIRST+1) 
RIGHT (I-PFIRST+2)  =  W(CS2,I-PFIRST+1) 
3310      CONTINUE 

WRITE(RT,8)  'VECTOR  NUMBER  ',CS2 

8  FORMAT(A17,I4) 
ELSE 

DO  3320  1=0, N 
3320        LEFT (I-PFIRST+2)  ■  W(CS1,I-PFIRST+1) 
RT  =  'NONE* 
ENDIF 

WRITE(LT,8)  'VECTOR  NUMBER  ',CS1 

HT  =  'WEIGHT  NUMBER' 

HU  =  '  ' 

NN  =  1 

XSCALE  =  1 

PNUM  =  N+l 

GOTO  3900 

C 

C  PLOT  FREQUENCY  RESPONSE  OF  ADF 

C 

3400    N2  =  1 

3410    N2  =  N2*2 

IF(N2.LT.N+1)  GOTO  3410 

WRITE*:*,*)  '  ' 

WRITE(*,6)  N+l,'  points,  need  to  zero  pad  to  at  least  *,N2,'  points' 
6       F0RMAT(I4,A,I4,A) 

TYPE  1,  '/Number  of  points  to  zero  pad  to  =   ' 
READ  (*,*)  N3 

N2  =  1 
3415    N2  =  N2*2 

IFCN2.LT.N3)  GOTO  3415 

WRITE(*,»)  'Zero  padded  to  *,N2,'  points' 

WRITE  (*,*)  '  ' 
WRITE  (*,9) 

9  FORMAT  (*  Frequency  response  desired  for:',/, 

X         '    0.   Mean  and  variance  of  weights',/, 
X         '    #.   Specific  vector  of  weights') 
TYPE  1,  '/Enter  0  or  vector  number:   ' 
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READ    (*,*)CS1 

IF    (CS1.NE.0)    THEN 
DO  3420   1=0, N 
3420  CT(I+1)    =  CMPLX(W(CS1,I),0.0) 

IF(N2.GT.N+1)    THEN 
DO  3430  I=N+1,N2-1 
3430  CW(I+1)    =  CMPLX(0.,0.) 

ENDIF 
ELSE 

DO  3440   1=0, N 

CWU+1)    =  CMPLX(WMEAN(I),0.0) 
CW2(I+1)    =  CMPLX(WVAR(I),0.0) 


3440      CONTINUE 


IF(N2.GT.N+1)  THEN 
DO  3450  I=N+1,N2-1 


CW(I+1)  =  CMFLX(0.,0.) 
CW2(I+1)  =  CMPLX(0.,0,) 
3450        CONTINUE 
ENDIF 
ENDIF 

CALL  DFT(C»,N2,0) 
CALL  DFT(CW2,N2,0) 

DO  3460  I=l,N2/2 

LEFT(I)  =  CABS(CWd)) 
RIGHT(I)  =  CABS(CW2(D) 
3460    CONTINUE 

IF(PFIRST.GT.N2/2)  THEN 

WRITE(»,»)  'THERE  ARE  ONLY  *,N2/2,'  RESPONSES' 

GOTO  3000 
ENDIF 

IF  (CS1.EQ.0)  THEN 

LT  =  'RESPONSE  OF  MEAN' 

RT  =  'RESPONSE  OF  VARIANCE' 
ELSE 

WRITE(LT.IO)  'RESPONSE  OF  VECTOR  ',CS1 
10        FORMAT (A19, 14) 

RT  =  'NONE' 
ENDIF 

HT  =  'FREQUENCY  RESPONCE' 
HU  =  'HERTZ' 
XSCALE  =  l./(N2»TS) 
WRITE(*,»)  •  ' 

TYPE  l,*/Plot  up  to  what  frequency?  (4800  maximum)   ' 
READ  (•,•)  FREQP 
PNUM  =  FREQP/XSCALE 
GOTO  3900 
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C  PLOT  MEAN  AND  VARIANCE  OF  WEIGHTS 

C 

3500    DO  3510  1=0, N 

LEFT(I+1)  =  WMEAN(I) 
RIGHTU+1)  =  WVAR(I) 
3510    CONTINUE 

LT  =  'MEAN  OF  WEIGHTS' 
RT  =  'VARIANCE  OF  WEIGHTS' 
HT  =  'WEIGHT  NUMBER' 
HU  =  '  * 
XSCALE  =  1. 
PNUM  =  N+l 
GOTO  3900 
C 

C  PLOT  ONE  WEIGHT  THROUGHOUT  THE  UPDATE  PROCEDURE 
C 

3600    TTPE  I,  '/History  of  what  number  weights?   (enter  a,b)   ' 
READ  (*,*)  CS1.CS2 
IF(CS2.GE.O)  THEN 
DO  3610  J-l, LENGTH 
LEFT(J)  =  W(J,CS1) 
RIGHT(J)  -  W(J,CS2) 
3610      CONTINUE 

WRITE(RT,7)  'CROSS-SECTION  NUMBER  ',CS2 
7         FORMAT (A21, 14) 
ELSE 

DO  3620  J=l, LENGTH 
3620        LEFT(J)  =  W(J,CS1) 
RT  ■  'NONE' 
ENDIF 

WRITE(LT,7)  'CROSS-SECTION  NUMBER  \CS1 

HT  =  'WEIGHT  CROSS-SECTION' 

HU  =  '  ' 

PNUM  =  J 

XSCALE  =  TS 

GOTO  3900 

3900    CONTINUE 

C 

C     FORM  SUB-TITLE  AND  CALL  PLOTTING  ROUTINE 

C 

BNK  =    '    ' 

IF   (WBIAS)    THEN 

WRITE(STITLE,5)  'ADF-DELAY=' .DELAY, ' ,  N=',N,',  MU=',MU,',  BIAS' 
5         FOHMAT(A10,I2,A4,I4,A5,F12.10,A6) 
ELSE 

WRITE(STITLE,5)  *ADF-DELAY=' .DELAY, ' ,  N»',N,',  MU=',MU 
ENDIF 

CALL   SPLOT(LEFT, RIGHT, PNUM, NN.LT, BNK, RT,  BNK, HT, 
X  HU, TITLE, STITLE, XSCALE) 

RBAD(»,2)    WAIT 
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GOTO  3000 
C 

C 

C  PARAMETER  CHANGE 
C 

4000    CALL  PADFDELAT (DELAY) 
GOTO  3000 

4100    CALL  PCONSTANT(MU) 
GOTO  3000 

4200    CALL  PWEIGHTS(N) 
GOTO  3000 

4300    CALL  PBIAS(WBIAS) 
GOTO  3000 

4400    WRITE<*,*)  'NOT  AVAILABLE  AT  THIS  TIME' 
GOTO  3000 

4500    CALL  PFILTER(FILT.BW) 
GOTO  3000 

9000    CONTINUE 

RETURN 

END 
C 
C 
C 
C**********************»f  ******************************************** 

c 

C  THE  FOLLOWING  SUBROUTINES  ARE  QUERIES  FOR  THE  ADAPTIVE  NOISE 

C  CANCELLER  PARAMETERS.   THEY  ARE  USED  IN  BOTH  THE  USER 

C  CONTROLLED  ANC  SUBROUTINE  AND  THE  MAIN  SIMULATION  ROUTINE. 

C 

SUBROUTINE  PARAMETERS 

REAL  MU.BW 

INTEGER  DELAY, N 

CHARACTER  CANS»1 

LOGICAL  WBIAS.FILT.OFILT,  SZC,  SWV,  SMM,  SCORR,  SANCZC 

1  FORMAT  (A) 

C 

C  NUMBER  OF  WEIGHTS  OF  WEIGHT  VECTOR 


C 


ENTRY  PWEIGHTS(N) 

WRITE(*,»)  •  • 

TYPE  1,  '/    Number  of  weights  = 

READ  (•,•)  N 


143 


RETURN 
C 

C  DELAY  OF  INPUT  TO  ADAPTIVE  FILTER 
C 

ENTRY  PADFDELAY  (DELAY) 

WRITE(«,»)    '    * 

TYPE  1,    '/  Reference    input   delay  =      ' 

READ    (*,*)    DELAY 

RETURN 
C 

C     RATE  OF  CONVERGENCE  FACTOR 
C 

ENTRY  PCONSTANT(MU) 

WRITE(«,»)    '    ' 

TYPE  1,  '/    Gradient  proportionality  constant  =   ' 

READ  (*,*)  MU 

RETURN 
C 

C  WEIGHT  BIAS  ENABLING 
C 

ENTRY  PBIAS(WBIAS) 

WRITE(*,»)    '    ' 

TYPE  1.  '/    Do  yon  want  to  include  bias  for  weights?   (Y/N)   ' 

READ  (»,2)  CANS 
2       FORMAT  (Al) 

IF  (CANS.EQ.'Y'.OR.CANS.EQ.'y')  THEN 
WBIAS  =  .TRUE. 

ELSE 

WBIAS  =  .FALSE. 

ENDIF 

RETURN 
C 

C  BANDWIDTH  ESTIMATOR  FOR  USER  CONTROLLED  ANC 
C 

ENTRY  PFILTER(FILT.BW) 

WRITE(*,»)  •  • 

TYPE  1,  '/    Do  yon  want  to  make  a  bandwidth  estimate  (Y/N)   ' 

READ  (*,2)  CANS 

IF  (CANS.EQ.'Y'.OR.CANS.EQ.'y')  THEN 
FILT  =  .TRUE. 

ELSE 

FILT  =  .FALSE. 

ENDIF 

RETURN 
C 

C  INTERMEDIATE  FILTERING  ENABLE  FOR  USER  CONTROLLED  ANC 
C 

ENTRY  PFILTKOFILT) 

WRITE(«,*)  *  • 

TYPE  1,  •/    Do  yon  want  to  enable  the  intermediate  filter?   (Y/N) 

READ  (»,2)  CANS 

IF  (CANS.EQ.'Y'.OR.CANS.EQ.'y')  THEN 
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OFILT  =  .TRUE. 

ELSE 

OFILT  =  .FALSE. 

ENDIF 

RETURN 
C 

C  ENABLE  ESTIMATES  BASED  ON  THE  WEIGHT  VECTOR 
C 

ENTRY  PHH(SHH) 

WRITE(*,»)  '  ' 

TTPE  1,  '/    Do  yon  want  to  enable  the  Mai/Min  estimators?   (T/N) 

READ  (*,2)  CANS 

IF  (CANS.EQ.'Y'.OR.CANS.EQ. 'y')  THEN 
SMM  =  .TRUE. 

ELSE 

SMM  -  .FALSE. 

ENDIF 

RETURN 
C 

C  ENABLE  ESTIMATES  BASED  ON  SIMILARITIES  BETWEEN  WEIGHT  VECTOR  AND 
C  AUTOCORRELATION  ESTIAMTE 
C 

ENTRY  PCORR(SCORR) 

WRITE(*.»)  '  • 

TYPE  3,'f         Do  yon  want  to  enable  the*, 
X     '/  psendo-corr.  estimators?   (Y/N)   ' 
3        FORMAT (A, A) 

READ  (*,2)  CANS 

IF  (CANS.EQ.'Y'.OR.CANS.EQ. 'y')  THEN 
SCORR  =  .TRUE. 

ELSE 

SCORR  ■  .FALSE. 

ENDIF 

RETURN 
C 

C  ENABLE  ESTIMATES  BASED  ON  ZERO  CROSSING  ANALYSIS 
C 

ENTRY  PANCZC(SANCZC) 

WRITE(*,»)  *  • 

TYPE  1,  •/    Do  yon  want  to  enable  the  ANC-ZC  estimators?   (Y/N) 

READ  (»,2)  CANS 

IF  (CANS.EQ.'Y'.OR.CANS.EQ. 'y')  THEN 
SANCZC  =  .TRUE. 

ELSE 

SANCZC  =  .FALSE. 

ENDIF 

RETURN 

END 
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C 

C  SSADF 

C 

C  VAX-11  FORTRAN  SOURCE  FILENAME:  SADF3 

C 

C       DEPARTMENT  OF  ELECTRICAL  ENGINEERING    KANSAS  STATE  UNIVERSITY 

C 

C       REVISION        DATE  PROGRAMMER  (S) 

C       

C        0.0  MAT  6,  1985  MARC  D.  BRACK 

C 

£•»«•«•»»»«•»•»»••••*»»«»*••»•»••••»»«»•••»••»•»«»»»«»••«»»»»»*«»*»»«»» 

c 

C       CALLING  SEQUENCE 

C 

C  CALL  SSADF(X,LENGTH,MU,N,DELAT,TOIAS,WMEAN,WVAR, 

C  FTLT,SZEST.TS,MMEST,BW,SMM, 

C  ERORI,YI.EROR,T,SCORR,SANCZC> 

C 

C  PURPOSE 

C 

C  This  subroutine  implements  the  Adaptive  Noise 

C  Canceller  Band  Rate  Estimator.   The  routine  is 

C  to  be  nsed  inconj unction  with  SBAUD. 

C 

C  ROUTINE(S)    ACCESSED  OR  CALLED  BY  THIS  ROUTINE 

C 

C  DFT 

C  BUTER 

C  ZCROSS 

C  BITLENGTB 

C 

C  ROUTINE(S)   ACCESSED  OR  CALLED  BY  EMBEDDED  ROUTINE(S) 

C 

C  NONE 

C 

C       ARGUMENT(S)  REQUIRED  FROM  THE  CALLING  ROUTINE 

C 

C  X       Array  of  input  sequence  (real) 

C 

C  LENGTH  Number  of  elements  of  X  (integer) 

C 

C  MU      Gradient  constant  for  LMS  algorithm  (real) 

C 

C  N       Number  of  weights  in  weight  vector  (integer) 

C 

C  DELAY   Number  of  samples  for  delay  (integer) 

C 

C  WBIAS   Flag  for  weight  bias,  TRUE  if  bias 

C  is  to  be  included  (logical) 

C 
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C  FILT 

C 

C 

C  TS 

C 

C  SMM 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C  HHEST 

C 

C 

C  BW 

C 

C  X 

c 
c 

C  EEORI 

C 

C 

C  TI 

C 

C 

C  EROR 

C 

C 

C  T 

C 

C 

C  WMEAN 

C 

C  WAR 

C 

c 


Flag  for  intermediate  filtering,  TRUE  if 
filtering  to  be  done  (logical) 

Sampling  interval  of  input  sequence  (real) 

Flag  for  bitlength  analysis,  TRUE  if 
analysis  to  be  done  (logical) 


SCORR   Flag  for  psendo-correlation  analysis,  TRUE 
if  analysis  to  be  done  (logical) 

SANCZC  Flag  for  zero  crossing   analysis,  TRUE  if 
analysis  to  be  done  (logical) 

ARGUHENT(S)  SUPPLIED  TO  THE  CALLING  ROUTINE 

SZEST   Array  for  estimates  made  with  secondary 
adaptive  noise  canceller  (real) 


Array  for  estimates  made  with  intermediate 
adaptive  noise  canceller  (real) 

Bandwidth  estimate  of  demodulator  output  (real) 

Array  of  intermediately  filtered  demodulator 
outpute  (real) 

Array  of  error  sequence  from  intermediate 
adaptive  noise  canceller  (real) 

Array  of  adaptive  filter  output  sequence  from 
intermediate  adaptive  noise  canceller  (real) 

Array  of  error  sequence  from  secondary 
adaptive  noise  canceller  (real) 

Array  of  adaptive  filter  output  suquence  from 
secondary  adaptive  noise  canceller  (real) 

Array  of  mean  of  filter  weights  (real) 

Array  of  variance  of  filter  weights  (real) 


SUBROUTINE  SSADF(X,LENGTH,MU,N,DELAY,WBIAS,WMEAN,WVAR, 
X  FILT,  SZEST, TS,MMEST,B¥,  SMM, 

X  ERORI,TI,EROR,T,  SCORR,  SANCZC) 

IMPLICIT  NONE 

INTEGER  N, I, J, LENGTH, DELAY,NUM,NUM2. MINIM, MINJV 

INTEGER  HFLAG(10),HISTA(2,10,257),MAX,ACCUM  . 
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REAL  MU,W(0:256)  ,WMEAN(0:N)  ,WVAR(0:N)  ,WMEANS(0:256)  ,MAXF 

REAL  X(LENGTH),SZEST(10),BW,MMEST(10),MINF,R(512) 

REAL  1(2048), EROR(2048)  ,TS,SBITLU0)  ,TEMPMU,L( 512) 

REAL  TI(2048),ERORI(2O48) 

COMPLEX  CW(2048),BUTER 

LOGICAL   SMM,FILT,WBIAS,TWO,INTERM,SCORR,SANCZC 

CHARACTER* 80  B 

INTERM  =    .FALSE. 
TEMPMII  =  MU 

IF(FILT)  MU  =  0.000000001 
C 

c 

C  ADAPTIVE  NOISE  CANCELLING  SECTION 

C 

1000    CONTINUE 

DO  100  1=0, N 
WMEAN(I)  =  0. 
WMEANS(I)  =  0. 
W(I)  =  0.0 
100     CONTINUE 

C 

C  LMS  IMPLEMENTATION  OF  ANC 

C 

DO  160  J=l, LENGTH 
IF  (TO I AS)  THEN 

Y(J)  =  W(0) 
ELSE 

T(J)  =  0.0 
ENDIF 

DO  130  1=1, N 
130  IFU-I+1.GT.0)   T(J)    =  T(J)    +  KJ-I+Df  (I) 

IF(J+DELAT.LE. LENGTH)  THEN 

EROR(J)  =  XCJ+DELAY)  -  T(J) 
ELSE 

EROR(J)  =  -T(J) 
ENDIF 

W(0)  =  W(0)  +  2.»MU»EROR(J) 

DO  140  1=1.  N 
140         IF(J-I+1.GT.0)  W(I)  =  W(I)  +  2.*MU»ER0R(J)«X(J-I+1) 
C 

C  END  OF  LMS  STRUCTURE 
C 

C  CALCULATE  MEAN  AND  MEAN  SQUARE  OF  WEIGHTS 
C 


148 


IF(J.GT.N)  THEN 
DO  150  1=0, N 

WMEAN(I)  -  WMEAN(I)  +  W(I)  /(LENG1H-N) 
WMEANS(I)  =  WMEANS(I)  +  W(I)**2/(LENGTH-N) 
150         CONTINUE 
ENDIF 

160     CONTINUE 

C 

C  END  OF  LMS  ALGORITHM  IMPLEMENTATION  LOOP 

C 

DO  170  1=0, N 
170       WVAR(I)  =  WMEANSd)  -  WMEAN(I)**2 
C 

C  END  OF  ANC  IMPLEMENTATION 
C 

£••»»•»»»•••*•»«*»»•»••»•••«•«»«••••«*»»»*»••»•••••»»•»•»•••*«««»•••«*» 

C 

C  OBTAIN  PSEUDO-CORRELATION  ESTIMATES 

C 

IF(SCORR)  THEN 
DO  5620  J=0,N 

CW(J+1)  =  CMPLX(WMEAN(N-J),0.)  •  FLOAT(J)/N 
5620      CONTINUE 

DO  5625  J=N+1,2*N 

Cff(J+l)    =   CMPLX(1IMEAN(J-N),0.)    •    (1.0  -  FLOAT(J-N)/N) 
5625  CONTINUE 

DO  5630  J=2»N+1,4*N-1 
5630  CW(J+1)    =   CMPLKO.,0.) 

CALL  DFT(CW,4«N,0) 
MINF  =  10000000 
DO  5640  J=1,2*N 

L(J)    =  20*ALOG10(CABS(CW(J))) 
IF(CABS(CW(J)).LT.MINF)    THEN 
MINIM  -  J 
MINF  =  CABS(CW(J)) 
ENDIF 
5640  CONTINUE 

DO  5660  J=0,N 

Cf(J+l)    =   CMPLX(WAR(N-J),0.)    •  FLOAT(J)/N 
5660  CONTINUE 

DO  5665  J=N+1,2*N 

CTU+1)  =  CMPLX(WVAR(J-N),0.)  •  (1.  -  FLOAT(J-N) /N) 
5665      CONTINUE 

DO  5670  J=2»N+1,4*N-1 
5670        CW(J+1)  =  CMPLX(0.,0.) 
CALL  DFT(C1T,4»N,0) 
MINF  =  10000000 
DO  5680  J=1,2*N 

R(J)  =  20»ALOG10(CABS(CW(J))) 
IF(CABS(CW(J)).LT.MINF)  THEN 
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MINJV  =  J 
MINF  =  CABS(CWU)) 
ENDIF 
5680      CONTINUE 

ENDIF 
C 

C 

C  OBTAIN  ESTIMATES  BASED  ON  MINIMUM  AND  MAXIMUM  OF  MEAN  AND  VARIANCE 

C  AND  FROM  ZERO  CROSSING  ANALYSIS 

C 

TWO  =  .TRUE. 

IF(SMM  .AND.  .NOT.INTERM)  THEN 

CALL  BITLENGTH(WMEAN,WVAR,N,DELAT,MMEST,TWO) 
MMEST(6)  =  FL0AT(2»(4*N))/MINJM 
MMEST(7)  =  FL0AT(2*(4»N)) /MINJV 
CALL  ZCROSS(EROR,LENGTH,MMEST(10).l) 
ENDIF 

IF(FILT  .AND.  INTERM)  THEN 

CALL  BITLENGTH(WMEAN,WVAR,N,DELAY,SZEST,TWO) 
SZEST(6)  =  FLOAT(2«(4*N))/MINJM 
SZESTC7)  =  FLOAT(2»(4»N)) /MINJV 
CALL  ZCROSS(X,LENGTH,MMEST(8),2) 
CALL  ZCROSS(I,LENGTH.SZEST(8),2) 
CALL  ZCROSS(EROR,LENGTH,SZEST(10),1) 
ENDIF 
C 

C 

C  IMPLEMENT  INTERMEDIATE  FILTERING.   THE  FILTER  BANDWIDTH  IS  BASED 

C  ON  THE  3  dB  BANDWIDTH  FOUND  BY  TAKING  THE  DFT  OF  THE  AVERAGE 

C  FILTER  WEIGHTS. 

C 

IF  (FILT  .AND.  .NOT.INTERM)  THEN 
INTERM  =  .TRUE. 
DO  210  1=0, N 
210         CW(I+1)  =  CMPLX(WMEAN(I).0.) 

DO  220  I=N+1,4*N-1 
220         CWU+1)  =  CMFLX(0.,0.) 
C 

C  DETERMINE  3  dB  BANDWIDTH  OF  SIGNAL 
C 

CALL  DFT(CW,4«N,0) 

MAXF  -  0 
DO  230  I=1.4»N/2 
230         IF(CABS(CW(I)).GT.MAXF)  MAXF  =  CABS(CWU)) 
MAXF  =  MAXF**2 

DO  240  I=l,4*N/2 
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IF(CABS(CW(4*N/2+l-I))*»2.GT.MAXF/2)    THEN 
BW  =    (4*N/2+2-I)/(4*N/2*TS) 
GOTO  250 
ENDIF 
240  CONTINUE 

250  CONTINUE 

C 

C     PERFORM  FILTERING 
C 

DO  320  J-l, LENGTH 
320  CW(J)    =   CMPLX«(J),0.0) 

CALL  DFT (CW, LENGTH, 0) 

DO  330  J-l, LENGTH/2 

CW(J)    ■  CW(J)»BUTER(4,BW,FL0AT(J-1)/(TS*LENGTH)) 

CW(LENGTH+1-J)  =  CW(LENGTH+1-J)  *BUTER(4,BW, FLOAT (-J)  /(TS*LENGTH) ) 
330       CONTINUE 

CW( LENGTH/ 2+1)  =  CMPLX(0.,0.) 

CALL  DFT(CW,LENGTH,1) 

DO  340  J-l, LENGTH 
ERORKJ)  =  EROR(J) 
YI(J)  =  Y(J) 
340         X(J)  =  REAL(CTKJ)) 
C 

C  ENTER  ANC  FOR  SPECIFIED  MU  AFTER  INTERMEDIATE  FILTERING  IS  COMPLETE 
C 

MU  =  TEMPMU 
GOTO  1000 

ENDIF 

RETURN 
END 
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C 

C  BITLENGTH 

C 

C  VAX-11  FORTRAN  SOURCE  FILENAME:  BIT.FOR 

C 

C       DEPARTMENT  OF  ELECTRICAL  ENGINEERING    KANSAS  STATE  UNIVERSITY 

C 

C       REVISION        DATE  PROGRAMMER(S) 

C 

C        0.0  MAT  6.    1985  MARC  D.  BRACK 

C 

c 

C       CALLING  SEQUENCE 

C 

C  CALL  BITLENGTH(WMEAN,WVAR,N,DELAT,BITL,TWO) 

C 

C       PURPOSE 

C  This  routine  makes  estiamtes  of  the  bad  rate  by 

C  estimating  the  bit  interval  in  number  of  sampling 

C  intervals.   The  estimate  is  taken  from  the  weight 

C  vector  information. 

C 

C       ROUTINE  (S)  ACCESSED  OR  CALLED  BY  THIS  ROUTINE 

C 

C  NONE 

C 

C       ROUTINE(S)  ACCESSED  OR  CALLED  BY  EMBEDDED  ROUTINE(S) 

C 

C  NONE 

C 

C       ARGUMENT(S)  REQUIRED  FROM  THE  CALLING  ROUTINE 

C 

C  WMEAN   Array  of  length  N  of  mean  of  the 

C  weight  vector  (real) 

C 

C  WAR    Array  of  length  N  of  variance  of  the 

C  weight  vector  (real) 

C 

C  N       Number  of  weight  of  the  weight  vector  (integer) 

C 

C  DELAY   Delay  of  the  primary  input  of  the  ANC  (integer) 

C 

C  TWO     TRUE  if  the  estimates  involving  the 

C  variance  are  to  be  made  also  (logical) 

C 

C       ARGUMENT(S)  SUPPLIED  TO  THE  CALLING  ROUTINE 

C 

C  BITL   Array  containing  estimates  (real) 

C 
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SUBROUTINE  BITLENGTH(WMEAN, WVAR.N, DELAY,  BITL, TWO) 

IMPLICIT  NONE 

INTEGER  MEXT(256),FRMIN,ZC 

INTEGER  RI,LI,DELAY,N,I,Q1,Q2,Q3,Q4,Q5,Q6,TEMP,FRMAX,L,R,J 

REAL  WMEAN(0:N),WVAR(0:N),BITL(10) 

LOGICAL  MINMAX,SIDED1,SIDED2,TW0,NSET,NSET1,NSET2,KEY 

DO  100  1=1,10 
100       BITL(I)  =  1  -  DELAY 
C 

C  FIND  ALL  RELATIVE  MINIMA  OF  THE  WEIGHT  MEAN  FILE  AND  ALL  RELATIVE 
C  MAXIMA  OF  THE  WEIGHT  VARIANCE  FILE 
C 

NSET  =  .TRUE. 
NSET1  =  .TRUE. 
NSET2  =  .TRUE. 
DO  200  1=1, N-l 

IF(WMEAN(I).LT.WMEAN(I-1)  .AND.  WMEAN(I)  .LT.WMEAN(I+1) )  THEN 
MEXT(I)  =  1 
IF (NSET)  THEN 
NSET  =  .FALSE. 
FRMIN  =  I 
ENDIF 
ELSEIF(WMEAN(I).GT.WMEAN(I-1)  .AND.  WMEAN(I)  .GT.WMEAN(I+1) )  THEN 
MErr(I)  =  -1 
IF(NSETl)  THEN 
NSET1  =  .FALSE. 
FRMAX  =  I 
ENDIF 
ELSE 

MEXT(I)  -  0 
ENDIF 
IF(.N0T.NSET1  .AND.  NSET2)  THEN 

IF(WMEAN(I)*WMEAN(I+1).LE.0.0)  THEN 
ZC  =  I 

NSET2  =  .FALSE. 
ENDIF 
ENDIF 
200     CONTINUE 

C 

C  ONE  ESTIMATE  IS  THE  FIRST  TIME  A  RELATIVE  MAXIMUM  AND  A  RELATIVE 

C  MINIMUM  OCCUR  AT  THE  SAME  POINT.   IF  THIS  DOES  NOT  HAPPEN  DEFAULT 

C  TO  FIRST  RELATIVE  MINIMUM  PAST  THE  FIRST  RELATIVE  MAXIMUM 

C 

IF  (TWO)  THEN 

DO  300  I=FRMAX,N 

IF(MEXTU).EQ.l)  THEN 
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CALL  FMAX(WVAR,I,0,N,KEY) 
IF  (KEY)  THEN 
BITL(l)  =  I 
GOTO  310 
ENDIF 
ENDIF 
300       CONTINIIE 

BITL(l)  =  FRMIN 
310       CONTINUE 
C 

C  A  SECOND  ESTIMATE  IS  IF  A  VARIANCE  MASIMUM  IS  WITHIN  TWO  SAMPLES  OF 
C  THE  MEAN  MINIMUM 
C 

DO  320  I=FRMAX,N 

IF(MErT(I).EQ.l)  THEN 

CALL  FMAX(WVAR,I,2,N,KEY) 
IF(KEY)  THEN 
BITL(2)  =  I 
GOTO  330 
ENDIF 
ENDIF 
320       CONTINUE 

BITL(2)  =  FRMIN 
330       CONTINUE 

ENDIF 
C 

C  THE  THIRD  ESTIMATE  IS  THE  ABSOLUTE  MINIMUM  OF  THE  MEAN 
C 

BITL(3)  =  FRMAX 

DO  400  I=FHMAX+1,N 
400       IF(WMEAN(I).LT.WMEAN(BITL(3)))  BITL(3)  =  I 
C 

C  THE  FOURTH  ESTIMATE  IS  THE  FIRST  RELATIVE  MINIMUM  LESS  THAN  ZERO 
C  IF  NOT  FOUND,  DEFAULT  TO  THE  FIRST  RELATIVE  MINIMUM 
C 

DO  500  I=FRMAX,N 

IF(MEXT(I).EQ.l  .AND.  WMEAN(I) .LT.O)  THEN 
BITL(4)  =  I 
GOTO  510 
ENDIF 
500     CONTINUE 

BITL(4)  =  FRMIN 
510     CONTINUE 
C 

C  THE  FIFTH  ESTIMATE  IS  THE  FIRST  ZERO  CROSSING  AFTER  THE  FIRST  RELATIVE 
C  MAXIMUM  OF  THE  MEAN 
C 

BITL(5)  =  ZC 
C 

C  ADD  THE  DELAY  DIFFERENCE 
C 
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DO   800   1=1,10 
800  BITL(I)    =  BITL(I)    +  DELAY  -  1 

RETURN 

END 
C 
C 

C 

c 

C  THIS  SUBROUTINE  IS  TO  DETERMINE  IF  A  MAXIMUM  OF  THE  VARIANCE 

C  IS  CLOSE  TO  THE  MINIMUM  OF  OF  THE  WEIGHT  MEAN. 

C 

C  The  minimum  of  the  mean  is  located  at  weight  number  I.   If  the 

C  maximum  occurs  within  0  samples  of  I  KEY  is  returned  with  the 

C  value  of  TRUE,  otherwise  it  is  returned  with  a  value  of  FALSE. 

C 

SUBROUTINE  FMAX(X,I,Q,N,KEY) 

IMPLICIT  NONE 
INTEGER  I,Q,N,J 
REAL  X(0:N) 
LOGICAL  KEY 

KEY  =  .FALSE. 

DO  100  J=I-Q,I+Q 

IF(J-l.GE.O  .AND.  J+l.LE.N-1)  THEN 

IF(X(J).GE.X(J-1)  .AND.  X(J)  .GE.XU+1) )  KEY  =  .TRUE. 
ENDIF 
100     CONTINUE 

\ 
RETURN 
END 
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C 

C  ZCROSS 

C 

C  VAX-11   FORTRAN  SODRCE  FILENAME:  ZC.FOR 

C 

C       DEPARTMENT  OF  ELECTRICAL  ENGINEERING    KANSAS  STATE  UNIVERSITY 

C 

C       REVISION        DATE  PROGRAMMER(S) 

C       

C        0.0  MAT  6,    1985  MARC  D.  BRACK 

C 

C 

C       CALLING  SEQUENCE 

C 

C  CALL  ZCROSS(SIG,N,ZEST,NEST) 

C 

C       PURPOSE 

C 

C  This  routine  forms  an  estimate  of  the  band  rate  using 

C  zero  crossing  analysis. 

C 

C       ROOTINE(S)  ACCESSED  OR  CALLED  BY  THIS  ROUTINE 

C 

C  NONE 

C 

C       ROUTINE(S)  ACCESSED  OR  CALLED  BY  EMBEDDED  ROUTINE(S) 

C 

C  NONE 

C 

C       ARGUMENT(S)  REQUIRED  FROM  THE  CALLING  ROUTINE 

C 

C  SIG     Array  of  length  N  of  input  sequence  (real) 

C 

C  N       Number  of  samples  of  input  (integer) 

C 

C  NEST    Number  of  estimates  to  makeu  the  different 

C  estimates  are  based  on  different  criteria. 

C  NEST  can  be  either:  1,  2,  or  3 .   Generally, 

C  all  estimates  give  same  result. 

C 

C       ARGUMENT(S)  SUPPLIED  TO  THE  CALLING  ROUTINE 

C 

C  ZEST    Array  of  length  NEST  of  estimates  (real) 

C 

c 

SUBROUTINE  ZCROSS(SIG,N,ZEST,NEST) 

IMPLICIT  NONE 

INTEGER  I,N,DIFF(2048),HIST(2048),J,K,SUCCES,.FAIL,NUM,NUM2.EST 
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INTEGER  MAXW.NEST 

REAL  SIG(N)  ,ZEST(NEST)  ,DIFFA(2,2048)  ,AMIN,TEMP,NORM,AMINT 

REAL  RATIO(10),TOL(10),LMAXW 

DO  100  1=1, N 
DIFF(I)  =  0 
HIST(I)  =  0 
100     CONTINUE 

1  =  1 

J  =  0 

K  =  1 
C 

C  DETERMINE  THE  LOCATIONS  OF  THE  ZERO  CROSSINGS  AND  RECORD  THE 
C  DIFFERENCES  IN  SAMPLES  BETWEEN  THE  CROSSINGS 
C 

200     DO  210  WHILE(SIG(I)»SIG(I+l).GT.O) 
1=1+1 

IF(I.GT.N-l)  GOTO  220 
210     CONTINUE 

215     IF  (SIG(I+1).EQ.0)  THEN 

1=1+1 

IFU.GT.N-1)  GOTO  220 

GOTO  215 
ELSE 

1=1+1 
ENDIF 

IF  (I.GT.N-1)  GOTO  220 

J=J+1 

DIFF(J)  =  I  -  K 

K  =  I 

GOTO  200 

220     CONTINUE 

C 

C  MAKE  A  HISTOGRAM  OF  THE  DIFFERENCES|i  THE  INDEX  OF  THE  HIST  ARRAY  IS 

C  THE  DIFFERENCE  BETWEEN  CROSSINGS,  THE  VALUE  OF  HIST  FOR  A  CERTAIN 

C  INDEX  IS  THE  NUMBER  OF  TIMES  A  DIFFERENCE  WHICH  EQUALS  TO  THE  INDEX 

C  OCCURRED 

C 

DO  300  1=1, J 
300       HIST(DIFFU))  =  HIST(DIFFd) )  +  1 
C 

C  MAKE  IMPULSE  FUNCTION  ON  REAL  LINE|i  THE  LOCATION  OF  THE  IMPULSES  ARE 
C  THE  AVERAGE  OF  A  GROUP  IN  THE  HISTOGRAM  WHERE  A  GROUP  IS  A  COLLECITON 
C  OF  ENTRIES  NOT  HAVING  TWO  OR  MORE  ADJACENT  EMPTY  BINS.   THE  WEIGHTS 
C  OF  THE  IMPULSES  ARE  THE  NUMBER  OF  HISTOGRAM  ENTRIES  THAT  CONTRIBUTED 
C  TO  THE  AVERAGE  OF  THE  RESPECTIVE  IMPULSE.   THE  ARRAY  DIFFA  CONTAINS 
C  THE  LOCATION  OF  THE  IMPULSES  IN  THE  FIRST  ROW  AND  THE  WEIGHT  IN  THE 
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C  SECOND  ROW 
C 

J=l 

NUM2  =  0 
DO  320  1=1, N 
C 

C  TO  DETECT  A  GROUP,  FIRST  FIND  A  NON-EMPTY  BIN  THAT  HAS  TWO  EMPTY  BINS 
C  IMMEDIATELY  PRECEEDING  IT  (CARE  MOST  BE  TAKEN  IN  CHECKING  THE  FIRST 
C  TWO  BINS) .   THEN  GATHER  IN  BIN  ENTRIES  FROM  THE  FIRST  BIN  OF  THE  GROUP 
C  UNTIL  TWO  ADJACENT  EMPTY  BINS  ARE  FOUND 
C 

IF(I.LE.2)  THEN 

IF(HIST(l).NE.O    .AND.    I.EQ.l)   THEN 
NUM  =  HIST(I) 
DIFFA(l.J)    =  HIST(I)*I 
K=l 
310  IF(.NOT.  (HISTU+KKEQ.O    .AND.    HISTU+K+1)  .EQ.O) )   THEN 

NUM  =  NUM  +  HIST(I+K) 

DIFFA(l.J)    =  DIFFA(l.J)    +  HIST(I+K)*(I+K) 
K=K+1 
GOTO  310 
ENDIF 

DIFFA(1,J)    =  DIFFA(1,J)/NUM 
DIFFA(2,J)    =  NUM 
NUM2   =  NUM2   +  NUM 
J=J+1 
ENDIF 

IF(HIST(2).NE.O    .AND.    HIST(l) .EQ.0    .AND.    I.EQ.2)THEN 
NUM  =  HIST(I) 
DIFFA(l.J)    =  HIST(I)*I 
K=l 
312  IF(. NOT.  <HIST(I+K). EQ.O    .AND.    HIST(I+K+1)  .EQ.O) )    THEN 

NUM  =  NUM  +  HIST(I+K) 

DIFFA(l.J)    =  DIFFA(l.J)    +  HIST(I+K)*(I+K) 
K=K+1 
GOTO  312 
ENDIF 

DIFFA(l.J)    ■  DIFFA(1,J)/NUM 
DIFFA(2,J)   -  NUM 
NUM2   =  NUM2   +  NUM 
J=J+1 
ENDIF 
ELSE 

IF(HIST(I).NE.O    .AND.    HIST(I-l) .EQ.O    .AND.    HIST(I-2) .EQ.O)    THEN 
,  NUM  =  HIST(I) 
DIFFA(l.J)    =  HIST(I)*I 
K=l 
315  IF(. NOT.  (HISTU+K). EQ.O    .AND.    HIST(I+K+1)  .EQ.O) )   THEN 

NUM  =  NUM  +  HIST(I+K) 

DIFFA(l.J)    =  DIFFAd.J)    +  HIST(I+K)  »(I+K) 
K=K+1 
GOTO  315 
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ENDIF 

DIFFA(l.J)  =  DIFFAd.D/NUM 
DIFFA(2,J)  =  NOM 
NUM2  =  NOM2  +  NDM 
J=J+1 
ENDIF 
ENDIF 
320     CONTINUE 

J=J-1 
AMIN  =  N 
DO  330  1=1, J 
330       AMIN  =  MIN(AMIN,DIFFA(1,I)) 

TOL(l)  =  0.05 

TOL(2)  =  0.20 

TOL(3)  =  0.10 

RATIO(l)  =  2 

RATIOU)  =  2 

RATIOC3)  =  2 

DO  390  EST=1,NEST 
C 

C  DETERMINE  THE  NUMBER  OF  SUCCESS  AND  FAILURES  OF  MEETING  THE  TOLERENCE 
C   CRITERION  WHERE  WHERE  A  SUCCES  IS  WHEN  THE  LOCATION  OF  THE  IMPULSE 
C  CONSIDERED  IS  AN  INTEGRAL  DIVISOR  OF  THE  A  LOCATION  TO  THE  RIGHT. 
C  EACH  IMPULSE  TO  THE  RIGHT  THAT  SATASFIES  THE  CRITERION  CONTRIBUTES  ITS 
C  WEIGHT  TO  THE  TOTAL  NUMBER  OF  SUCCESS 
C 

335       SUCCES  =  0 
FAIL  =  0 
DO  340  1=1, J 

IF(DIFFA(1,I).GE.AMIN)  THEN 
TEMP  =  DIFFA( 1,1) /AMIN 

NORM  =  (TEMP  -  ANrNT(TEMP))/ANINT(TEMP) 
IF(ABS(NORM).LE.TOL(EST))  THEN 
SUCCES  =  SUCCES  +  DIFFA(2,I) 
ELSE 

FAIL  =  FAIL  +  DIFFA(2,I) 
ENDIF 

IF(DIFFA(1,I).EQ.AMIN)  NUM  =  DIFFA(2,I) 
ENDIF 
340       CONTINUE 
C 

C  CHECK  THE  RATIO  CRITERION.   THERE  MUST  BE  A  CERTAIN  RATIO  BETWEEN  THE 
C  NUMBER  OF  SUCCESS  AND  NUMBER  OF  FAILURES.   IN  ADDITION  THE  NUMBER  OF 
C  SUCCESS  MUST  BE  GREATER  THAN  ONE-FOURTH  OF  THE  TOTAL  NUMBER  OF 
C   CROSSINGS  AND  THE  NUMBER  OF  ENTRIES  CONTRIBUTING  TO  THE  IMPULSE 
C  BEING  CONSIDERED  MUST  BE  GREATER  THAN  ONE-TENTH  OF  THE  TOTAL  NUMBER 
C  OF  CROSSINGS 
C 

IF(SUCCES.GE.RATTO(EST)»FAIL  .AND.  SUCCES. GE.NUM2/ 4  .AND. 
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X        NnM.GT.NUM2/10)  THEN 
ZEST(EST)  =  AMIN 
ELSE 

AMINT  =  AMIN 
AMIN  =  N 
DO  360  1=1, J 
360  IF(MIN(DIFFA(1.I),AMIN).GT.AMTNT)  AMIN  =  MIN(DIFFA(1,I)  ,AMIN) 

IF(AMIN.LT.N)  THEN 

GOTO  335 
ELSE 
C 

C   IF  NO  IMPOLSE  MEETS  THE  CRITERIA,  CHOOSE  THE  MAXIMUM  WEUGHT  IMPULSE 
C  AS  THE  DEFAULT  ESTIMATE 
C 

MAXW  =  0 

DO  370  1=1, J 

MAXW  =  MAX(DIFFA(2,I),FLOAT(MAXW)) 
370  CONTINUE 

DO  380  1=1, J 

IF  (MAXW.EQ.DIFFAU.D)  THEN 
LMAXW  =  DIFFAU.I) 
GOTO  385 
ENDIF 
380  CONTINUE 

385  ZEST(EST)  =  LMAXW 

ENDIF 
ENDIF 
390     CONTINUE 

RETURN 
END 
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C 

C  DLD 

C 

C  VAX-11  FORTRAN  SOURCE  FILENAME:  DLD. FOR 

C 

C       DEPARTMENT  OF  ELECTRICAL  ENGINEERING    KANSAS  STATE  UNIVERSITY 


160 


c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


DLD 

VAX-11  FORTRAN  SOURCE  FILENAME: 
DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
REVISION        DATE 


DLD.  FOR 

KANSAS  STATE  UNIVERSITY 

PROGRAMMER(S) 


0.0 
0.1 


JAN  23,  1983 
MAY  6,  1985 


DONALD  M.  HUMMEL S 
MARC  D.  BRACE 


CALLING  SEQUENCE 

CALL  DLD  (TD,FC,FLAG,S,N,TF) 

PURPOSE 

This  subroutine  models  a  delay  line  discriminator. 
No  pre-  or  post  filtering  is  done.  The  resulting 
waveform  is  returned  in  the  input  array. 

ROUTINE(S)  ACCESSED  OR  CALLED  BY  THIS  ROUTINE 

DFT 
ROUTINE(S)  ACCESSED  OR  CALLED  IN  EMBEDDED  ROUTINES 

NONE 

ARGUMENT(S)  REQUIRED  FROM  THE  CALLING  ROUTINE 

TD=Time  delay  of  discriminator 

FOCenter  frequency 

FLAG=Integer  representing  type  of  discriminator 

0  —  Cos  type 

1  —  Sin  type 

S=Array  containing  input  signal 

N=#  of  points  in  input  array  (a  power  of  2) 

TF=Input  period  (=1/ fundamental  freq) 

ARGUMENT (S)  SUPPLIED  TO  THE  CALLING  ROUTINE 

S-Array  containing  output  signal 


1«1 


COMPLEX  S(2048),S2(2048), FACTOR 

INTEGER  FLAG 
C 

C       PERFORM  THE  TIME  DELAY 
C 

DO  19  I  -  1,N 

S2(I)    =   S(I) 

19  CONTINUE 

CALL  DFT(S2,N,0) 

DO  20   I  =  l.N/2 

ANGLE  =  -(2»3.141592»FLOAT(I-l)/TF»TD) 
FACTOR  =  CMFLX(COS(ANGLE),SIN(ANGLE)) 
S2(I)    =  S2(I)    •  FACTOR 
ANGLE  -  -(2»3.141592*FLOAT(-I)/TF»TD) 
FACTOR  =  CMPLX(COS(ANGLE),SIN(ANGLE)) 
S2(N+1-I)    =  S2(N+1-I)    •  FACTOR 

20  CONTINUE 

CALL  DFT(S2,N,1) 
C 

C     FIND  THE  INPUT  FOR  THE  POST-FILTER 
C 

ANGLE  =  3.14159*(2*FC»TD  +    .5*FLAG) 

FACTOR  =   CMPLX(COS(ANGLE),SIN(ANGLE)) 

DO  40   1=1, N 

S(I)    =  S(I)*CONJG(S2(I))*FACTOR 
S(I)    =   CMPLX(4.»REAL(S(I)),0.) 
40  CONTINUE 

RETURN 
END 
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Abstract 


The  traditional  method  of  FSK  signal  band  rate  estimation  is 
to  nse  a  delay-line  discriminator  to  demodulate  the  signal  and 
then  perform  a  zero- crossing  analysis  to  obtain  the  band  rate. 
This  stndy  proposes  an  alternate  method  which  implements  an 
adaptive  noise  canceller.  The  weight  vector  of  the  adaptive 
noise  canceller,  npon  convergence,  can  be  used  directly  to  deter- 
mine the  band  rate.  However,  the  input  signal  is  random  making 
convergence  difficult  to  obtain  or  sustain.  The  adaptive  noise 
canceller  can  be  used  indirectly  to  improve  the  performance  of 
the  zero-crossing  analysis.  The  way  in  which  it  does  this  is  by 
making  a  rough  but  consistent  estimate  of  the  bandwidth  of  the 
message.  This  parameter  is  not  known  due  to  the  surveillence 
situation  the  receiver  is  intended  to  be  used  in.  Then  the 
delay-line  discriminator  output  can  be  filtered  to  produce  a 
cleaner  input  for  the  zero-crossing  analysis  which  will  improve 
perf ormace . 

The  performance  of  the  traditional  z e r o- c r o s s i ng  analysis 
and  the  direct  and  indirect  implementation  of  the  adaptive  noise 
canceller  was  judged  using  a  Monte  Carlo  simulation.  The  results 
showed  that  both  adaptive  noise  cancelling  methods  performed 
better  than  the  zero-crossing  analysis  alone. 


